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RESUMEN 



Se presenta el estudio de dos edificios esbeltos de 21 pisos con base en de marcos ductiles 
de acero estructural, disenados para distorsiones de entrepiso cercanas al limite maximo A=1.2% 
estipulado por el Reglamento del Distrito Federal vigente en Mexico (RCDF-93). Los edificios no 
cumplen con la relation de esbeltez limite H/L< 2.5 estipulada en el apartado seis de las Normas 
Tecnicas Complementarias por Sismo vigentes (NTCS‘95), por lo que el factor reductivo de 
fuerzas sismicas con fines de diseno (Q') utilizado para el diseno fue afectado por un factor 
reductivo de 0.8 (Q -0.8Q), segun se especifica en la section 4.1 de las NTCS‘95. Las estructuras 
se disenaron para la zona III definida por las NTCS‘95, es decir, la zona de lago o de terreno 
blando. Se realizaron analisis dinamicos no lineales de marcos ante la action de acelerogramas 
registrados y artificiales asociados al espectro de diseno. Se comparan los niveles y distribution 
de fluencia estructural y las demandas de ductilidad maximas de cada nivel con respecto al factor 
de comportamiento sismico (Q) utilizado en el diseno. Se comparan, ademas, las distorsiones 
maximas de diseno con las distorsiones dinamicas maximas obtenidas de los analisis no lineales. 



Del estudio se desprende que aunque las demandas de ductilidad de entrepiso maximas 
se encuentran dentro de niveles razonables y que la fluencia de elementos estructurales 
corresponde al mecanismo ultimo de falla de viga debibcolumna fuerte, una interpretation 
erronea de las NTCS'95 con respecto a como deben revisarse las restricciones de desplazamientos 
laterales para fines de diseno puede conducir a distorsiones dinamicas maximas de entrepiso del 
orden de 2.4%, que sobrepasan notablemente el limite maximo de distorsion de entrepiso 
permitido por el reglamento. A estos niveles de distorsion se pueden esperar danos considerables 
en elementos no estructurales, asi como problemas de funcionamiento y operation del edificio 
asociado a la conducta humana durante sismos. 



SUMMARY 



The study of two 21 story, slender, special moment resisting steel buildings is presented. 
Both buildings were designed according to the seismic provisions of RCDF-93 code of Mexico City 
for a maximum story drift ratio close to the limiting drift A=1.2% established by the code. As both 
buildings do not satisfy the limiting slenderness ratio H/L< 2.5 established by RCDF-93 for 
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regular structures, then, the Q 'factor allowed by the code to reduce seismic forces for design 
purposes was affected by a 0.8 reduction factor, as established by RCDF-93 (Q-0.8Q'). Both 
buildings were designed for the lakebed region of Mexico City (soft soil conditions). Nonlinear 
dynamic analyses of representative frame models of subject structures were conducted. Several 
recorded and simulated accelerograms associated to the design spectra of RCDF-93 for the 
lakebed region were used for the nonlinear dynamic analyses. Response modification factors (Q) 
used in the design of subject buildings were compared against story ductility computed from the 
nonlinear dynamic analyses. Story drift ratios associated to the original design according to 
RCDF-93 were compared with peak dynamic story drift angles computed from nonlinear dynamic 
analyses. Structural yielding was studied and associated to hysteretic, deformation and strength 
demands. 

The study shows that although peak story ductility demands are within reasonable 
bounds, and the yielding mapping suggests the potential formation of a weak-beam, strong- 
column failure mechanism, peak dynamic drift angles are much larger than anticipated in the 
design (up to 2.4%). For drift angles higher than 2%, one should expect considerable damage in 
non- structural elements, important problems in the normal operation of the building, as well as 
potential panic in people, particularly people living in the lake-bed zone of Mexico City, where 
strong earthquakes are felt for a time interval of two or three seconds. Most residents of Mexico 
City that felt the M s =8.1 September 19, 1985 earthquake are particularly sensitive and still 
afraid of moderate and strong earthquakes felt in the city. 

1. INTRODUCTION 



Las Normas Tecnicas Complementarias para Diseno por Sismo (NTCS-95) del 
Reglamento del Distrito Federal vigente (RCDF-93) establece, en su apartado seis, las once 
condiciones que debe cumplir una estructura para ser considerada como regular y, de esta 
manera, utilizar directamente el factor reductivo de fuerzas sismicas con fines de diseno (Q 1 ) en 
funcion de los factores de comportamiento sismico (Q) del apartado cinco. En caso de que no se 
cumpla una, varias o todas las condiciones de regularidad del apartado seis, Q' debera 
multiplicarse por un factor reductivo de 0.8, segun se especifica en la section 4.1 de las NTCS-95. 
Esta disposition reglamentaria esta basada exclusivamente en la intuition y no se cuenta con 
estudios que la validen. Ademas, las normas de referenda no indican como proceder en caso de 
que se violen simultaneamente varias condiciones de regularidad, en cuyo caso el castigo de 0.8 
del factor reductivo puede no ser suficiente. 

A1 respecto, los comentarios a las NTCS‘87, que aplican a las NTCS‘95 puesto que no han 
sufrido modificaciones en este criterio, dicen lo siguiente- “La reduction de Q’ establecida para 
edificios que no cumplen con requisites de regularidad refleja el resultado de la experiencia y de 
la intuition. En efecto, se espera mejor comportamiento ante sismo de los edificios que satisfacen 
dichos requisites! ademas, un examen de los efectos de temblores intensos recientes confirma en 
terminos generates la menor vulnerabilidad de tales inmuebles. Sin embargo, no existen analisis 
ni pruebas de laboratorio que fundamenten cuantitativamente la reduction de Q’ en un 20% para 
los edificios menos regulares, asi que esta disposition no necesariamente implica que dichas 
estructuras desarrollan menores factores de ductilidad” (Rosenblueth y Gomez 1991). 

Por lo tanto, es evidente que ni siquiera puede asegurase que el factor Q' este relacionado 
con la esbeltez, o con alguna de las otras condiciones de regularidad, sino que la reduction de Q’ 
es simplemente una medida preventiva arbitraria. En tal caso, es importante reflexionar sobre 
las siguientes cuestiones (entre muchas otras)- ante la action de sismos, <-,que debilidades 
presenta una estructura por el hecho de ser esbelta? ^.Podra reducirse la vulnerabilidad ante 
sismos de una estructura esbelta con solo modificar el valor propuesto del factor Q’? De hecho, 
^podra realmente disminuirse la vulnerabilidad de estructuras irregulares reduciendo el valor de 
Q’ cuando no se satisface ninguna condition de regularidad? ^Seria mas prudente tomar medidas 
preventivas diferentes en el diseno de estructuras irregulares dependiendo de la condition de 
irregularidad que no se satisfaga, como lo hace, por ejemplo, en el reglamento UBC (Uniform 
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Building Code)? <^Es adecuado el factor de 0.80 para estructuras que violan la condition de 
esbeltez u otras de las condiciones de regularidad? La contestation a estas interrogantes requiere 
de muchos estudios detallados, de los cuales el estudio que aqui se presenta es solamente un 
modesto esfuerzo que intenta dar respuesta a algunas de ellas. 

Es un hecho que los efectos de las irregularidades pueden afectar el comportamiento de la 
estructura de diferentes maneras, y la influencia de estos no ha sido cuantificada suficientemente 
en estudios previos a nivel mundial realizados en las decadas de los anos setenta, ochenta y 
principios de los noventa, sobre todo en el intervalo de comportamiento inelastico, aunque en el 
ultimo lustro se han realizando algunos estudios para condiciones de irregularidad tales como 
piso debil, estructuras esbeltas, estructuras alargadas, irregularidades importantes en planta y/o 
elevation, asi como irregularidades en la distribution de masas, rigideces y resistencia (por 
ejemplo, Valmundsson y Nau 1997, Ruiz et al. 1997, Tena et al. 1996, Tena 1997, Tena-Colunga y 
del Valle 1998, Tena 2000, Tena 2001 y Tena, 2002). El estudio comparative de las coincidencias 
y diferencias principales entre las disposiciones de diseno sismico para estructuras irregulares de 
los reglamentos mexicanos y el reglamento norteamericano UBC se ha presentado recientemente 
en la literatura (Tena-Colunga y del Valle 1998, Tena-Colunga 1999). 

Con respecto al diseno de estructuras irregulares, cabe senalar que las Normas Tecnicas 
Complementarias para Diseno por Sismo para el Distrito Federal de 2004 (NTCS‘2004 2004), que 
estan por entrar en vigor, contemplan algunos cambios en el diseno de estructuras irregulares. 
Entre los mas importantes esta la definition de estructuras fuertemente irregulares, que 
esencialmente son aquellas con un piso suave o una excentricidad torsional muy pronunciada. 
Para estas estructuras, Q' debe multiplicarse por un factor de 0.7. Otra modification es que si 
una estructura no satisface una condition de regularidad, Q' se multiplicand por 0.9 y por 0.8 
cuando no cumpla con dos o mas de las condiciones de regularidad del apartado 6 de las NTCS- 
2004, que son esencialmente los mismos que se definen en las NTCS‘95 y NTCS‘87. 

En el presente trabajo se resume brevemente los estudios realizados en estructuras 
esbeltas de 21 pisos a base de marcos ductiles de acero estructural, disenados para distorsiones 
de entrepiso cercanas al limite maximo A=1.2% estipulado por el RCDF-93 que no cumplen con la 
relation de esbeltez limite H/L< 2.5 estipulada en el apartado seis de las NTCS‘95, y que se 
presentan en detalle en trabajos mas extensos (Cheja 1997, Tena et al. 1996 y Tena 1997), donde 
ademas se presenta el estudio realizado para edificios de ocho pisos con relaciones de esbeltez 
similares. El diseno de cada estructura se hizo utilizando el factor de comportamiento sismico 
adecuado conforme a las NTCS‘95, como se presenta brevemente en las siguientes secciones. Este 
estudio es tambien representative de la propuesta de NTCS-2004, ya que de acuerdo con ellas, el 
diseno sismico por condition de irregularidad seria el mismo que se establece en las NTCS‘95, 
ademas que los criterios especificos para el diseno de marcos ductiles de acero estructural son 
esencialmente los mismos, sobre todo los que rigen el diseno de vigas y columnas. 



2. ESTRUCTURAS EN ESTUDIO 



En este estudio se eligio como estructura esbelta prototipo al proyecto arquitectonico de 
una de las torres de unos condominios de lujo ubicados en la calle Ruben Dario en la colonia 
Polanco, en una de las zonas residenciales mas exclusivas de la ciudad de Mexico. Se considero 
que este tipo de proyecto es representative de estructuras esbeltas comunes que pueden 
construirse en los proximos anos en distintas zonas de la ciudad de Mexico, incluyendo la zona del 
lago. 



El conjunto condominal cuenta con un terreno con dimensiones de 64x57 m. y un area de 
3648 m 2 , y esta compuesto por tres torres (A, B y C) de 21 pisos y 66 m. de altura sobre la 
superficie del terreno, y de tres sotanos comunes donde se alojan las areas de estacionamiento y 
servicios (SI, S2 y S3). Por su parte, las torres A y B son iguales entre si, mientras que la torre C, 
aunque presenta la mismas altura que las dos anteriores, es un edificio mas regular tanto en 
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planta como en elevation. Los detalles del proyecto se presentan en Cheja (1997); sin embargo, a 
continuation se resumen las caracteristicas generales de la torre A, cuya geometria fue 
seleccionada para el presente estudio. 

La torre A esta conformada por un mezanine en el primer piso, y diez departamentos de 
dos niveles cada uno en los pisos subsecuentes. El mezanine tiene una altura de 4.00 m., la planta 
baja de cada departamento es de 3.40 m. y la planta alta de 2.80 m. Los elevadores y las escaleras 
de servicio solo tienen entrada en la planta baja de cada condominio, y una escalera interior 
ubicada entre los ejes C, D, 8 y 9 (figura l) permite el acceso a la planta alta. Todas las losas 
tienen un saliente entre los ejes 7 y 8 de forma sensiblemente triangular y, excepto en la azotea, 
otro entre los ejes 8 y 9 donde se aloja una de las jardineras. La jardinera ubicada entre los ejes C 
y D solo aparece en la planta baja de cada departamento, ya que en dicha area existe doble altura 
y, por lo tanto, en la losa de entrepiso es espacio vacio. El area entre los ejes D, E, 6 y 7 esta 
destinada para las escaleras de servicio. Los claros principales tienen longitudes entre 3.80 m. el 
mas pequeno y 6.00 m. el mas grande. Las dimensiones gruesas de la planta son de 15.9 x 20.0 m, 
siendo esta simetrica respecto a un eje central paralelo al lado mayor de 20.0 m. 




Se eligio la geometria de la torre “A” como prototipo por ser adecuada para el estudio de 
estructuras esbeltas, ya que no cumple con la condition de regularidad por esbeltez (condition 2 
NTCS-95). La condition de esbeltez no se cumple en ambas direcciones en planta (figura l), ya 
que para la direction X, la relation de esbeltez es 66/15.9=4.15, mientras que para la direction Y 
dicha relation resulta igual a 66/20.0=3.30; ambas mayores al limite de 2.5. Ademas, la torre A 
tampoco cumple con las condiciones de regularidad 6, 7, y 11 de las NTCS-95, segun se demuestra 
con detalle en Cheja (1997). 
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2.1 Diseno sismico de los edificios IR6 e IR6B 



Se realizo el diseno sismico de dos modelos, denominados IR6 e IR6B, con base en el 
proyecto arquitectonico de la torre “A”, con la salvedad que en el presente estudio se considera 
que la estructura tiene condition de empotramiento perfecto en su base, lo que obedece a que se 
pretende estudiar exclusivamente la validez de las disposiciones debidas a la irregularidad de la 
estructura, tales como la esbeltez! despreciando otras acciones que pueden ser significativas como 
la interaction suekrestructura. Cabe senalar que el que no se consideren la interaction suekr 
estructura y otros efectos no implica de ninguna manera que no se consideren importantes, mas 
bien, el razonamiento asociado a esta decision es que el incluir otros efectos (como la interaction) 
no permitiria discernir claramente en que medida afectan las condiciones de irregularidad 
(particularmente la esbeltez), por si solas, al diseno y a la respuesta dinamica no lineal de 
estructuras disenadas conforme a las NTCS‘95, que es realmente uno de los objetivos centrales 
que se persiguen con los estudios que se estan conduciendo sobre este tema. 

Los modelos IR6 e IR6B se disenaron considerando que pertenecen al grupo “B” segun el 
RCDF‘93, y como tales fueron considerados los factores de carga para el diseno. Los valores del 
espectro de diseno utilizados para el analisis sismico son los correspondientes a la zona III (figura 
2), es decir, la zona de mayores ordenadas espectrales de diseno, y en la cual resulta mas 
interesante evaluar a las NTCS‘95. La estructuracion considerada son marcos ductiles a 
momento, por lo que el factor de comportamiento sismico empleado fue Q=4. Debido a que no 
cumplen con las condiciones de regularidad antes destacadas, el factor reductivo de fuerzas 
sismicas con fines de diseno (Q') se multiplied por el factor reductivo de 0.8 especificado en la 
section 4.1 de las NTCS‘95, y dadas las caracteristicas dinamicas de los modelos (T<Tb), en este 
caso Q'=0.8Q=4x0.8=3.2. Para las restricciones de desplazamientos laterales que indica el 
articulo 209 del RCDF‘93, se acepto que para las estructuras en cuestion los elementos no 
estructurales estan totalmente desligados de la estructura principal, por tanto, se considero como 
limite de la distorsion de entrepiso de diseno el valor de A=1.2%. 




Figura 2. Espectros de diseno para estructuras del grupo B, zona III del RCDF-93 



La diferencia basica entre el diseno de IR6 e IR6B estriba exactamente en como se revisan 
las restricciones de desplazamientos laterales para fines de diseno (Ad), ya que para algunos la 
redaction actual de las NTCS‘95 no es totalmente clara al respecto. De acuerdo con la section 4.1 
de las NTCS‘95, “las deformaciones laterales se calcularan multiplicando por Q las causadas por 
las fuerzas sismicas reducidas cuando se emplee el metodo estatico de analisis que se detalla en 
la section 8 de las presentes normas 0 el analisis modal de la section 9”. Si llamamos Aq- a las 
deformaciones laterales causadas por las fuerzas sismicas reducidas por Q’, entonces, de acuerdo 
con las NTCS‘95, A D = QAq , que es el criterio con que se diseno IR6B. En cambio, IR6 se diseno 

considerando que A D =Q'A q -, ya que existe quien interpreta que el castigo por el efecto de 
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irregularidad es reducir Q por 0.8, por lo que consideran logico revisar las deformaciones 
suponiendo que Q'=Q, que no es lo que se estipula en las NTCS‘95 y, de hecho, las NTCS‘2004. 
Por tanto, en este caso concrete, se considero que A D = 3.2A Q . para IR6 y A D = 4A Q - para IR6B. De 

hecho, de acuerdo con las NTCS‘95, las estructuras irregulares deben disenarse para resistir 
fuerzas sismicas reducidas mayores que las estructuras regulares y, si se observa con detalle, el 
criterio de revision para deformaciones laterales, aunque pareciera ser el mismo que para las 
estructuras regulares es, de hecho, mas estricto. Este criterio se mantiene en la propuesta de 
NTCS-2004. 

En la estructuracion con base en marcos de acero ductiles se utilizaron secciones cajon 
para las columnas, pues sus propiedades son mas convenientes para dar estabilidad ante carga 
lateral y evitar problemas de pandeo y, ademas se usan frecuentemente en la practica mexicana. 
Se seleccionaron secciones I (perfiles tipo W) para las vigas principales y secundarias, verificando 
que cumplieran los criterios para secciones compactas correspondientes. Como sistema de piso se 
propuso una losa de lamina de acero acanalada y concreto reforzado con malla electrosoldada de 5 
cm. de espesor, soportada por trabes secundarias. El diseno estructural de IR6 e IR6B se llevo a 
cabo con apoyo del programa ETABS. Se realizaron las siguientes hipotesis en el modelado de las 
estructuras. El modulo de elasticidad del acero se tomo como E = 2.039x10® Kg./cm^. Se 
utilizaron las propiedades geometricas nominales de los perfiles W en las vigas y las propiedades 
geometricas efectivas en las columnas de section cajon. La zona rigida que se considero tanto 
para las columnas como para las trabes es del 50% de la longitud total de apoyo entre los 
elementos, segun se recomienda para este tipo de edificios. Se realizo un diseno conforme a un 
analisis modal espectral elastico tridimensional conforme a las NTCS‘95 y al RCDF-93. En el 
analisis estructural fueron tornados en cuenta los efectos secundarios (P-A) que se inducen a la 
estructura por los desplazamientos laterales debidos al sismo, los cuales deben de incluirse en el 
diseno de estructuras esbeltas. 

El diseno preliminar de los elementos estructurales principales de acero se hizo conforme 
al reglamento norteamericano “Load and Resistance Factor Design” (LRFD-93), utilizando el 
modulo STEELER del ETABS. Posteriormente, se verified que las columnas y trabes mas 
esforzadas de dicho diseno, cumplieran con las especificaciones de las Normas Tecnicas 
Complementarias Para Diseno y Construction de Estructuras Metalicas (NTCEM-95, 1995); como 
se ilustra en detalle en Cheja (1997). Los detalles de los disenos se presentan en Cheja (1997), 
Tena et al. (1996) y Tena (1997) y se resumen en las tablas 1 y 2, donde se presentan las secciones 
finales de diseno para columnas y trabes, respectivamente, y en las figuras 3 y 4, donde se 
presentan los niveles de esfuerzo de trabes y columnas asociados a los disenos de IR6 e IR6B, 
respectivamente. 



Tabla 1 . Secciones finales de diseno de columnas 



Modelo IR6 


Modelo IR6B 


Niveles 


Columnas 


Dimensiones (cm.) 


Niveles 


Columnas 


Dimensiones (cm.) 


1-2 


C1-C5, C8, C9, 
C12, C13, C16 a 
C20 


70x60x2.54 


1-4 


C1-C20 


70x70x2.86 


1-2 


C6, C7, CIO, Cll, 
C14, C15 


50x50x2.54 


3-10 


C1-C20 


50x50x2.54 


5-10 


C1-C20 


60x60x2.38 


11-16 


C1-C20 


45x40x2.54 


11-16 


C1-C20 


50x50x1.91 


17-21 


C1-C20 


40x40x1.27 


17-21 


C1-C20 


40x40x1.59 



99 



Comportamiento sismico de edificios esbeltos con base en marcos ductiles de acero 
estructural disenados conforme al reglamento de la Ciudad de Mexico 



Tabla 2. Secciones finales de diseno para trabes 



Niaioveles 


Ceaw2556 


Modelo IR6B 


Trabes 


Perfil 


Trabes 


Perfil 


1-9 


T-l, T-3 a T-15 


W24xl03 


T-l, T-3 a T-15 


W24xl31 


T-2 


W24xll7 


T-2 


W24xl46 


T-16 a T-31 


W24x84 


T16 a T-31 


W24xl03 


10-12 


T-l, T-3 a T-15 


W24xl03 


T-l, T-3 a T-15 


W24xl31 


T-2 


W24xll7 


T-2 


W24xl46 


T-16 a T-31 


W24x76 


T-16 a T-31 


W24x94 


13 


T-l, T-3 a T-15 


W24x84 


T-l, T-3 a T-15 


W24xll7 


T-2 


W24xl03 


T-2 


W24xl31 


T-16 a T-31 


W24x68 


T-16 a T-31 


W24x94 


14-16 


T-l, T-3 a T-15 


W24x76 


T-l, T-3 a T-15 


W24xll7 


T-2 


W24x84 


T-2 


W24xll7 


T-16 a T-31 


W24x68 


T-16 a T-31 


W24x94 


17 


T-l, T-3 a T-15 


W24x62 


T-l a T-15 


W24x94 


T-2, T-16 a T-31 


W24x68 


T-16 a T-31 


W24x94 


18-21 


T-l, T-3 a T-31 


W24x62 


T-l a T-15 


W24x94 


T-2 


W24x68 


T-16 a T-31 


W24x76 



En las figuras 3 y 4 se presentan las distorsiones de entrepiso de diseno calculadas para 
las dos direcciones principales de los modelos IR6 e IR6B respectivamente, y se comparan con los 
limites establecidos por el RCDF-93. Tambien se ilustran los niveles de esfuerzo del diseno final 
de trabes y columnas. Puede observarse que en los niveles 6 a 14 se presentan distorsiones en la 
direccion X muy cercanas al limite de 1.2%; y que tanto las columnas de los primeros niveles, 
como las trabes paralelas al eje Y de los niveles intermedios, se encuentran sometidas a un alto 
nivel de esfuerzos. Se observa que la direccion critica ambos modelos es la direccion X. Es 
evidente tambien que muchos de los elementos se encuentran sobredisenados, en gran medida, 
como consecuencia de tomar secciones tipo, ya que como se aclaro anteriormente, se pretende que 
el diseno sea lo mas simple, realista y construible posible, con distorsiones cercanas al limite 
maximo permitido por el RCDF-93, pero no necesariamente optimo en cuanto a peso y/o nivel de 
esfuerzos. En este respecto, parece que el diseno de IR6B es ligeramente mas optimo en cuanto al 
nivel de esfuerzos de los elementos estructurales que el de IR6, conforme a las acciones 
estipuladas por el Reglamento. 




Figura 3. Niveles de esfuerzo y distorsiones de diseno, modelo IR6 



Arturo Tena Colunga y Carlos Cheja Mochon 
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Figura 4. Niveles de esfuerzo y distorsiones de diseno, modelo IR6B 



2.2 Caracteristicas dinamicas elasticas de los modelos IR6 e IR6B 



Los resultados del analisis modal del diseno final de los modelos IR6 e IR6B se resumen 
escuetamente en la tabla 3, donde se indica tambien el peso total de cada estructura. Ambos 
modelos son flexibles, pero un poco “rlgidos” (valga la aparente contradiction) para ser 
estructuras de marcos de acero de 21 niveles y 66 m de altura, debido a que se considero que los 
muros de bloque de concreto que se ubican en el area del cubo de elevadores participan en la 
rigidez lateral initial de las estructuras. Las estructuras son mas flexibles en la direction E _ W, en 
la cual son mas esbeltas. Los modos de traslacion estan moderadamente acoplados, debido a las 
importantes excentricidades de masa que existen por el uso de areas de los entrepisos y a la 
ubicacion de los 



Tabla 3. Caracteristicas dinamicas de los modelos IR6 e IR6B 



MODO 


Direction 


Modelo IR6 
(Wi6=5542.7 Ton) 


Modelo IR6B 
Wi6B=5704.5 Ton) 


T (s) 


Masa modal (%) 


T (s) 


Masa modal (%) 


N-S 


E-W 


9 


N-S 


E-W 


0 


1. Primer modo 

traslacion 


E-W 


1.923 


4.64 


63.02 


6.45 


1.646 


3.59 


63.10 


5.79 


2. Primer modo 

traslacion 


N-S 


1.882 


70.51 


4.51 


0.07 


1.615 


70.08 


3.50 


0.05 


3. Primer modo torsion 


Rotation 


1.352 


0.11 


5.64 


69.45 


1.197 


0.07 


5.22 


69.04 


4. Segundo modo 
traslacion 


E-W 


0.660 


2.20 


10.49 


0.83 


0.566 


4.03 


9.47 


0.83 


5. Segundo modo 
traslacion 


N-S 


0.646 


11.33 


2.25 


0.14 


0.559 


9.92 


4.10 


0.30 


6. Segundo modo torsion 


Rotation 


0.479 


0.03 


1.91 


10.88 


0.430 


0.02 


1.89 


11.32 



muros de bloque de concreto (figura l). Los modos de traslacion en direction N _ S presentan menor 
acoplamiento en rotation que los modos de traslacion en direction E _ W. Los modos de torsion 
presentan mayor acoplamiento con traslacion en la direction E-W. El peso total de IR6 es 
Wi6=5,542.7 ton y el de IR6B es Wjgg=5,704.5 ton. El peso de los elementos estructurales 
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principales es aproximadamente 14% el peso total de la estructura para IR6 y 16% para IR6B. Se 
observa que el interpretar con rigor las disposiciones de las NTCS‘95 para revisar las 
restricciones de desplazamientos laterales para fines de diseno (Ad) en el diseno de estructuras 
irregulares guio al diseno de una estructura mas pesada y rigida (IR6B) con respecto al diseno 
resultante con la interpretacion que muchos otros estructuristas le dan a dichas normas (IR6). 
Sera importante evaluar si esta diferencia en la interpretacion puede guiar a diferencias de 
comportamiento notables entre ambos modelos ante la accion de sismos intensos. Este es un 
aspecto sobre el cual reflexiona el presente trabajo. 



3. ANALISIS NO LINEALES 



Los disenos de los modelos IR6 e IR6B, conforme a los lineamientos de las NTCS‘95 y 
NTCEM-95 del RCDF-93, se evaluaron mediante la simulacion de su respuesta ante la accion de 
acelerogramas asociados al sismo de diseno realizando analisis dinamicos no lineales en el 
programa DRAIN-2DX (Prakash et al. 1992). 

Debido a que en la investigacion en curso uno de los objetivos es estudiar si los factores de 
comportamiento sismico para el diseno de estructuras irregulares en la zona III de las NTCS'95 
son indicados, es necesario utilizar acelerogramas representatives cuyas caracteristicas esten 
relacionadas al sismo maximo esperado y contemplado por el presente reglamento y 
acelerogramas que pudieran estar contemplados en reglamentos futures. Por lo tanto, se 
seleccionaron para el estudio de IR6 e IR6B cinco acelerogramas, dos registrados en el valle de 
Mexico durante el sismo del 19 de septiembre de 1985 (M g =8.l) y tres acelerogramas sinteticos 
generados para un sismo de M s =8.2 para tres estaciones de la red acelerografica del valle de 
Mexico. 



Los acelerogramas registrados durante el sismo del 19/09/85 que fueron seleccionados 
corresponden a los sitios Secretaria de Comunicaciones y Transportes en su componente E _ W 
(SCT-EW) y Tlahuac bombas en su componente N-S (TBOM-NS). Las caracteristicas de los 
registros difieren en cuanto a duracion, contenido de frecuencias, periodo dominante y respuestas 
espectrales maximas. Los registros de SCT tienen una duracion de un poco mas de 160s y 
aceleracion maxima del terreno de aproximadamente 0.17g en su componente E _ W, con respuesta 
espectral maxima de l.Og para un periodo dominante de 2.0s. El registro de Tlahuac bombas en 
direccion N-S tiene una duracion de 150s con espectro de respuesta con periodos dominantes 
cercanos a los 2.0s, con aceleraciones maximas del terreno de aproximadamente 0.14g y 
aceleraciones espectrales maximas de 0.75g. 

Los registros artificiales para estaciones de la red acelerografica de la ciudad de Mexico se 
generaron para un sismo postulado de M g =8.2, basado en las funciones de transferencia promedio 
calculadas para cada estacion con base en la informacion capturada por la red acelerografica a 
partir de 1988, conforme a los procedimientos expuestos en Perez Rocha (1998). Se seleccionaron 
tres estaciones cuyos niveles de respuesta espectral de pseudoaceleraciones resultaran similares 
o mayores al acelerograma registrado en SCT-EW durante el sismo del 19 de septiembre de 1985 
y que tuvieran periodos dominantes de respuesta mayores a 1.3s y menores a 2.2s dadas las 
caracteristicas de las estructuras a estudiar. Se selecciono, por tanto, la estacion 05 (Las Cibeles), 
ubicada en la colonia Condesa, con periodo dominante de 2.15s y aceleracion maxima espectral de 
0.86g; la estacion 56, ubicada en la calle de Cordoba, en la colonia Roma, con periodo dominante 
de 2.20s y aceleracion maxima espectral de 1.44g (figura 5), y la estacion 84 (Culhuacan), ubicada 
en San Francisco Culhuacan, al sur de la ciudad de Mexico, con periodo dominante de 1.38s y 
aceleracion maxima espectral de 1.03g. Por fines de espacio no se muestran todos los registros y 
sus caracteristicas espectrales, pero pueden consultarse en otros trabajos (por ejemplo, Cheja 
1997, Tena et al. 1996, Tena 1997). 
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Los analisis dinamicos no lineales se condujeron exclusivamente en los marcos mas 
esforzados en cada direccion (marco 10 en la direction X, marco B en la direction Y), es decir, se 
obtuvieron marcos dinamicamente equivalentes con el procedimiento explicado en detalle en 
Cheja (1997). Se procedio de esta manera, ya que en estudios paralelos que se condujeron en 
edificios esbeltos de acero de ocho niveles cuya geometria se definio de manera que fuera analoga 
a la de IR6 e IR6B (Cheja 1997, Tena et a/1996), se observaron respuestas dinamicas muy 
similares en analisis conducidos en marcos equivalentes con respectos a modelos que consideran 
la interaction entre los marcos que resisten a las fuerzas sismicas en una direccion dada, por lo 
que se decidio, para una primera evaluation, economizar en el esfuerzo y recursos de computo 
necesarios para la evaluation de la respuesta dinamica de IR6 e IR6B. 




En los analisis no lineales se considero un nivel de amortiguamiento viscoso equivalente 
^=3%, ya que el nivel de amortiguamiento viscoso equivalente medido experimentalmente en 
algunas estructuras de acero varia entre C,=2% y ^=4%. Se considero que los elementos 
estructurales son capaces de desarrollar sus capacidades plasticas teoricas y que se presenta un 
comportamiento bilineal del acero estructural, donde inicialmente se considero una rigidez 
posterior a la fluencia del 2% de la rigidez elastica, que es la minima observada en elementos de 
acero estructural en pruebas experimentales y que ha sido reportada indirectamente en la 
literatura (por ejemplo, Mamaghani et al. 1995). En estudios subsecuentes se han considerado 
tambien pendientes posteriores a la fluencia de 5% y 10%, para evaluar su influencia en la 
respuesta dinamica no lineal, resultados que tambien se discuten en este trabajo. 

Se estudiaron los siguientes parametros dinamicos para los marcos 10 (direccion X) y B 
(direccion Y) de los modelo IR6 e IR6B: (a) curvas de histeresis de entrepiso, (b) mapeo de la 
fluencia de elementos, en tiempos especificos, en envolventes, y su monitoreo en el dominio del 
tiempo, (c) envolventes de distorsiones de entrepiso maximas (A), (d) envolventes de cortantes de 
entrepiso maximos (V/Wt), (e) demandas maximas de ductilidad de entrepiso (p), (f) distorsion de 
entrepiso asociada a la primera fluencia de cualquier elemento estructural que forma parte del 
entrepiso (Af e ), (g) distorsion de entrepiso maxima asociada al medio ciclo de histeresis de carga 
maxima (A me ), (h) cociente del numero de ciclos donde se presento respuesta inelastica (Cine) entre 
el numero total de ciclos a que fue sujeta la estructura ante un acelerograma dado (Ctot), (i) 
cociente entre la rigidez de entrepiso secante minima inelastica (kimin) 0 “pico a pico” asociada a 
los ciclos histereticos de amplitud maxima y la rigidez elastica del entrepiso (k e i) y, (j) cociente 
entre el promedio de las rigideces de entrepiso secantes (kiprom) asociadas a los medios ciclos 
histereticos inelasticos (kmc) y la rigidez elastica del entrepiso (k e i). Algunos de estos conceptos se 
ilustran en la figura 6. 

La information analizada es muy extensa, ya que se realizaron analisis para los marcos y 
modelos de referenda ante la action de todos los acelerogramas senalados anteriormente, 
considerando pendientes posteriores a la fluencia de 2%, 5% y 10%, y por razones de espacio en 
este trabajo solamente se presentan y discuten algunos resultados del estudio descrito. Las 
mayores respuestas para los marcos 10 y B de los modelos IR6 e IR6B se obtuvieron ante la 
action de los acelerogramas SCT85-EW y S56-EW; sin embargo, solo se presentan los resultados 
de los marcos 10 de los modelos IR6 (figura 8) e IR6B (figura 7) y de los marcos B de IR6 (figura 
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10) e IR6B (figura 9) ante la action del acelerograma S56-EW, ya que, en general, las respuestas 
mas criticas se asocian a este acelerograma. 




De las envolventes de respuesta maximas del marco 10 del modelo IR6B (el que cumple 
con todas las disposiciones de las NTCS'95 para edificios irregulares) se observa que su respuesta 
no lineal es reducida (figura 7), con demandas maximas de ductilidad de 1.3 asociadas a rigideces 
secantes minimas inelasticas (pico a pico) del 73% de su rigidez elastica (k2=2%, nivel 5), a un 
valor minimo observado del promedio de las rigideces secantes inelasticas de medio ciclo del 83% 
de su rigidez elastica (k 2 = 2 %, nivel l) y distorsiones de entrepiso maximas muy por debajo de las 
calculadas en su diseno. Se presento respuesta elastica en numerosos niveles. En aquellos 
entrepisos donde se presento respuesta inelastica, la distorsion de entrepiso asociada a la 
primera fluencia del mismo vario de Af e =0.31% (k2=2%, nivel l) a Af e =0.65% (k2=2%, nivel 8), 
obteniendose en algunos entrepisos respuestas elasticas para distorsiones tan altas como 
A=0.76% (k2=2%, nivel 12), lo que llama poderosamente la atencion. La distorsion de entrepiso 
maxima asociada al medio ciclo de histeresis de carga maxima fue A me =0.97% en los niveles 4 a 6, 
la cual se encuentra por debajo del limite maximo permitido por el RCDF-93 (A=1.2%). 




Figura 7. Envolventes de respuesta maxima, marco 10 del modelo IR6B, acelerograma S56-EW 
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Como se puede observar de la figura 4, el diseno de IR6B en esta direction, a pesar de ser 
mas flexible, esta asociado a niveles de esfuerzos de las vigas inferiores, en promedio, a aquellos 
calculados para la direction Y. La fluencia se presento exclusivamente en vigas, y observando al 
nivel de esfuerzos de diseno asociados al diseno, no se puede explicar el por que existe tanta 
variabilidad en la distorsion asociada a la primera fluencia entre los distintos entrepisos (Af e ), 
aspecto que se analizara con mayor detalle en estudios futuros. Ademas, para este edificio no se 
observan grandes diferencias en las respuestas dinamicas maximas del marco 10 con respecto a 
la hipotesis hecha con respecto a la pendiente posterior a la fluencia de los elementos 
estructurales (2%, 5% o 10%). 

El cortante basal que resiste el marco 10 ante la action del acelerograma S56-EW varfa 
entre el 5.8% y 6.0% del peso total del edificio IR6B, por lo que crudamente se puede estimar que 
la resistencia a cortante basal de IR6B en direction X debe encontrarse entre el 29% (k2=2%) y 
30% (ka=10%) del peso total de la estructura, ya que en esta direction participan cinco marcos de 
caracterfsticas muy similares en cuanto a estructuracion y rigidez. Por tanto, si se toma en 
cuenta que la estructura se diseno para Q'=3.2 y en la meseta del espectro de diseno para zona III 
(c=0.40, figura 2), entonces, el cortante sismico no factorizado en una direction de analisis serfa 

c 

crudamente estimado como V = — IE = 0.125JE , por lo que la sobrerresistencia observada para 

este diseno en direction X, ante la action del acelerograma de referenda, debe variar crudamente 
entre 2.3 y 2.4, ignorando aun otras fuentes posibles de sobrerresistencia como la action 
tridimensional de la estructura y la contribution en la rigidez y resistencia del sistema de piso. 
En estudios futuros se pretende cuantificar de una manera mas precisa la sobrerresistencia que 
se presenta en esta y otras estructuras. 

De las envolventes de respuesta maximas del marco 10 del modelo IR6 (donde se estudia 
una variation en la interpretation de las disposiciones de las NTCS‘95 de como revisar las 
restricciones de desplazamientos laterales para fines de diseno, Ad, para edificios irregulares), 
figura 8, se observa que su respuesta no lineal, aunque reducida, es mas importante que la 
observada en el modelo IR6B. La demanda maxima de ductilidad de entrepiso es de 2.26 (k2=2%, 
nivel 5) asociada a rigideces secantes minimas inelasticas (pico a pico) del 43% de su rigidez 
elastica (k2=2%, nivel 5), a un valor minimo observado del promedio de las rigideces secantes 
inelasticas de medio ciclo del 68% de su rigidez elastica (k 2 = 2 %, nivel 12). Las distorsiones de 
entrepiso maximas sobrepasan a las calculadas en su diseno en los niveles 3 a 9, y al limite 
maximo permitido por el RDF-93 A=1.2%. La distorsion maxima corresponde al sexto entrepiso y 
es A=1.58%, calculada esta distorsion con respecto a la position inicial de equilibrio. Se presento 
respuesta elastica exclusivamente en los niveles 18 y 20. En aquellos niveles donde se presento 
respuesta inelastica, la distorsion de entrepiso en que se presento su primera fluencia vario de 
Af e =0.30% (k2=2%, nivel l) a Af e =0.97% (k2=2%, nivel 16), obteniendose en algunos entrepisos 
respuestas elasticas para distorsiones tan altas como A=1.07% (k2=2%, nivel 18), coincidiendo 
esta gran variation con lo observado en IR6B. La distorsion de entrepiso maxima asociada al 
medio ciclo de histeresis de carga maxima fue A me =2.2% en los niveles 5 y 6, valor que se 
encuentra muy por encima del limite maximo permitido por el RCDF-93 (A=1.2%). 

La fluencia se presento principalmente en vigas, aunque las columnas exteriores de los 
niveles 1, 3 y 4 experimentaron fluencias incipientes (se presentaron tres veces en 2500 pasos de 
integration analizados). Al igual que para IR6B, observando al nivel de esfuerzos asociados al 
diseno, no se puede explicar el por que existe tanta variabilidad en la distorsion asociada a la 
primera fluencia entre los distintos entrepisos (Afe). Al igual que para IR6B, no se observan 
grandes diferencias en las respuestas dinamicas maximas del marco 10 de IR6 con respecto a 
pendientes posteriores a la fluencia de los elementos estructurales de 2%, 5% o 10%. El cortante 
basal que resiste el marco 10 de IR6 ante la action del acelerograma S56-EW varfa entre 5.1% 
(k2=2%) y 5.5% (k2=10%) del peso total del edificio IR6. Siguiendo un razonamiento analogo al 
presentado para IR6B, se puede estimar crudamente que la resistencia a cortante basal de IR6 en 
direction X debe encontrarse entre el 25% y 28% del peso total de la estructura, y como en este 
caso el cortante sismico no factorizado en una direction de analisis serfa crudamente estimado 
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tambien como V = 0.125W , la sobrerresistencia observada para este diseno en direccion X, ante la 
accion del acelerograma de referenda, debe variar crudamente entre 2.0 y 2.25. 
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Figura 8. Envolventes de respuesta maxima, marco 10 del modelo IR6, acelerograma S56-EW 



Se puede concluir que la respuesta dinamica del marco 10 del modelo IR6B, que se diseno 
completamente conforme a las disposiciones de las NTCS‘95 del RCDF-93 para estructuras 
irregulares, es mas favorable en todos los rubros que la del modelo IR6, donde se estudio una 
variation en la interpretation de las disposiciones de las normas para la revision de las 
restricciones de desplazamientos laterales para fines de diseno, que como puede observarse a 
partir de la discusion y de las figuras 7 y 8, conduce a desplazamientos y distorsiones de 
entrepisos mayores a los permitidos por el Reglamento, que si bien no estan asociados a dano 
estructural o demandas inelasticas severas, a estos niveles de distorsion se pueden esperar danos 
considerables en elementos no estructurales, asi como problemas de funcionamiento y operation 
del edificio asociados a la conducta humana durante sismos que deben ser tornados en cuenta por 
el proyectista. 

En la direccion Y se presentan mayores demandas no lineales que en la direccion X ante 
la accion del acelerograma S56'EW, principalmente en los niveles inferiores, a juzgar de 
comparar las respuestas de los marcos B de IR6B e IR6 (figuras 9 y 10) con las de los marcos 10 
discutidas anteriormente (figuras 7 y 8). De las envolventes de respuesta maximas del marco B 
del modelo IR6B se observa que su respuesta no lineal es reducida (figura 9), con demandas 
maximas de ductilidad de 2.1 asociadas a rigideces secantes minimas inelasticas (pico a pico) del 
54% de su rigidez elastica (k2=2%, nivel 5), a un valor minimo observado del promedio de las 
rigideces secantes inelasticas de medio ciclo del 91% de su rigidez elastica (k2=2%, niveles 3 y 4) y 
distorsiones de entrepiso maximas que sobrepasan las calculadas en su diseno en los niveles 1 a 8 
(k2=1.17%, nivel 5), pero que no sobrepasan el limite del limite maximo de distorsion (A=1.2%) 
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permitido por el Reglamento. Se presento respuesta elastica en los niveles superiores 15 a 21. En 
aquellos entrepisos donde se presento respuesta inelastica, la distorsion de entrepiso asociada a 
la primera fluencia del mismo vario de Af e =0.28% (k2=2%, nivel l) a Af e =0.61% (k2=2%, nivel 12), 
obteniendose en algunos entrepisos respuestas elasticas para distorsiones de hasta A=0.43% 
(k2=2%, niveles 15 y 18). La distorsion de entrepiso maxima asociada al medio ciclo de histeresis 
de carga maxima fue A me =0.97% en los niveles 4 a 6, la cual se encuentra por debajo del limite 
maximo permitido por el RCDF-93 (A=1.2%). La distorsion de entrepiso maxima asociada al 
medio ciclo de histeresis de carga maxima fue Ame=1.4% en los niveles 4 y 5, valor que se 
encuentra por encima del limite maximo permitido por el RCDF-93 (A=1.2%). Como se puede 
observar de la figura 4, el diseno de IR6B en esta direccion, es mas rigido y esta asociado a 
niveles de esfuerzos mayores para las vigas, en promedio, a aquellos calculados para la direccion 
X. 
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Figura 9. Envolventes de respuesta maxima, marco B del modelo IR6B, acelerograma S56-EW 



La fluencia se presento exclusivamente en vigas, y observando al nivel de esfuerzos 
asociados al diseno, aqui tampoco se puede explicar el por que existe tanta variabilidad en la 
distorsion asociada a la primera fluencia entre los distintos entrepisos (Af e ). Aunque no se 
observan grandes diferencias en las respuestas dinamicas maximas del marco B con respecto a la 
hipotesis hecha con respecto a la pendiente posterior a la fluencia de los elementos estructurales 
(2%, 5% o 10%), si son un poco mas apreciables que para el marco 10 (figura 7) y, en general, 
tienden a disminuir las demandas de ductilidad y a incrementar las rigideces secantes inelasticas 
y la distorsion de entrepiso asociada a la primera fluencia. 

El cortante basal que resiste el marco B ante la accion del acelerograma S56-EW varia 
entre el 7.5% (k2=2%) y 7.8% (k2=10%) del peso total del edificio IR6B, por lo que crudamente se 
puede estimar que la resistencia a cortante basal de IR6B en direccion Y debe encontrarse entre 
el 30% y 31.2% del peso total de la estructura, ya que en esta direccion participan cuatro marcos 
de caracteristicas muy similares en cuanto a estructuracion y rigidez. Por tanto, la 
sobrerresistencia observada para este diseno en direccion Y, ante la accion del acelerograma de 
referenda, debe variar crudamente entre 2.4 y 2.5. 
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Figura 10. Envolventes de respuesta maxima, marco B del modelo IR6, acelerograma S56-EW 



Al igual que para los marcos 10, de las envolventes de respuesta maximas del marco B del 
modelo IR6 (figura 10) se observa que su respuesta no lineal es mas importante que la observada 
en el modelo IR6B (figura 9). La demanda maxima de ductilidad de entrepiso es de 2.5 (k2=2%, 
nivel 5) asociada a rigideces secantes minimas inelasticas (pico a pico) del 40% de su rigidez 
elastica (k2=2%, nivel 5), a un valor minimo observado del promedio de las rigideces secantes 
inelasticas de medio ciclo del 72% de su rigidez elastica (k2=5% y 10%, nivel 5). Las distorsiones 
de entrepiso maximas sobrepasan a las calculadas en su diseno en los niveles 1 a 14, y al limite 
maximo permitido por el RDF-93 A=1.2% en los niveles 3 a 7. La distorsion maxima corresponde 
al cuarto entrepiso y es A=1.63% (k2=10%) calculada esta distorsion con respecto a la position 
inicial de equilibrio. Se presento respuesta elastica exclusivamente en los niveles 19 a 21. En 
aquellos niveles donde se presento respuesta inelastica, la distorsion de entrepiso asociada a su 
primera fluencia vario de Af e =0.33% (k2=2%, nivel l) a Af e =0.53% (k2=2%, nivel 14), obteniendose 
en algunos entrepisos respuestas elasticas para distorsiones de A=0.43% (k2=2%, nivel 19), 
coincidiendo esta variation con la observada en IR6B para el marco B. El intervalo de variation 
observado en los marcos B es menor al observado en los marcos 10. La distorsion de entrepiso 
maxima asociada al medio ciclo de histeresis de carga maxima fue Ame=2.4% en los niveles 4 y 5 
(k2=10%), valor que se encuentra muy por encima del limite maximo permitido por el RCDF-93 
(A=1.2%). 



La fluencia se presento principalmente en vigas, excepto de una columna exterior del 
nivel 1, que experimento fluencia incipiente en su base (de una a cinco veces en 2500 pasos de 
integration analizados). Se pueden hacer las mismas observaciones que se han hecho para el 
marco B de IR6B con respecto a las respuestas dinamicas maximas, dependiendo de las 
pendientes posteriores a la fluencia de los elementos estructurales de 2%, 5% o 10%. 
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El cortante basal que resiste el marco B de IR6 ante la action del acelerograma S56-EW 
varla entre 6.6% (k2=2%) y 7.4% (k2=10%) del peso total del edificio IR6, por lo que se puede 
estimar crudamente que su resistencia a cortante basal en direction Y debe variar entre el 26% y 
30% del peso de IR6, y que la sobrerresistencia observada para este edificio debe variar 
crudamente entre 2.08 y 2.40. Por tanto, se puede concluir, al igual que marco el marco 10, que la 
respuesta dinamica del marco B del modelo IR6B es mas favorable en todos los rubros que la del 
marco B del modelo IR6, por las mismas razones expuestas anteriormente. 



4. COMENTARIOS FINALES 



Se presento el estudio de estructuras esbeltas de 21 niveles de acero estructural, 
disenadas para distorsiones de entrepiso de diseno cercanas al limite maximo A=1.2% permitido 
por el RCDF-93, para marcos ductiles. El diseno de cada estructura irregular se realizo utilizando 
que el factor reductivo de fuerzas sismicas con fines de diseno (Q') establecido por las NTCS'95. 
Se estudiaron dos modelos, el modelo IR6B que cumple con todas las disposiciones de diseno de 
las NTCS'95 para estructuras irregulares y el modelo IR6, donde se estudio una variation en la 
interpretation de las disposiciones de las NTCS'95 de como deben revisarse las restricciones de 
desplazamientos laterales para fines de diseno (Ad). Se efectuaran analisis dinamicos no lineales 
empleando acelerogramas representatives para la zona del lago de la ciudad de Mexico, cuyas 
caracteristicas esten relacionadas al sismo maximo esperado y contemplado por el reglamento 
vigente (acelerogramas registrados durante el sismo del 19/09/85) y acelerogramas artificiales 
que pudieran estar contemplados en reglamentos futuros pues, de hecho, el sismo de M s =8.2 
proveniente de las costas de Guerrero es uno de los escenarios sismicos que se estudiaron en la 
revision de los espectros de diseno para la propuesta de NTCS'2004, que aun no entran en vigor. 
Los resultados obtenidos para los modelos en estudio permiten concluir, de manera preliminar, lo 
siguiente. 

Conforme a lo observado en el estudio del modelo IR6B, las disposiciones de diseno para 
estructuras irregulares de las NTCS'95 permiten llegar a disenos razonablemente seguros para 
estructuras esbeltas con base en marcos de acero estructural aun cuando estos se disenen para 
altas distorsiones (cercanas al limite maximo A=1.2% estipulado por el RCDF'93) y factores de 
comportamiento sismico altos (Q=4), ya que aun para acelerogramas asociados a sismos 
postulados intensos, la respuesta no lineal es moderada, estable y concentrada en las vigas. Las 
demandas maximas de ductilidad de entrepiso (p) son pequenas e incluso, bastante inferiores al 
factor de comportamiento sismico Q asociado a su diseno. Cabe senalar que aunque p y Q no son 
exactamente el mismo concepto, si existe relation entre ambos parametros, y en este respecto, 
este estudio y otros que conduce el primer autor estan encaminados tambien para establecer 
relaciones entre ambos parametros. Sin embargo, como se observo para el marco B, existe la 
posibilidad que las distorsiones de entrepiso maximas sobrepasen a las calculadas en su diseno, 
aunque en el caso particular de IR6B no se observo que sobrepasaran al limite maximo A=1.2% 
permitido por el RCDF'93. 

Como se indico en secciones anteriores, las disposiciones de diseno por irregularidad de la 
propuesta de NTCS'2004 que competen a los edificios en estudio se mantienen iguales, por lo que 
tambien deben esperarse disenos seguros con la futura norma para este tipo de edificios. 

Conforme a lo observado en el estudio del modelo IR6, la variation estudiada en la 
interpretation de las disposiciones de las NTCS'95 de como deben revisarse las restricciones de 
desplazamientos laterales para fines de diseno (Ad) tambien permiten llegar a disenos 
razonablemente seguros para estructuras esbeltas con base en marcos de acero estructural en 
cuanto a niveles de esfuerzos y demandas no lineales, aun cuando estos se disenen para altas 
distorsiones (cercanas al limite maximo A=1.2% estipulado por el RCDF'93) y factores de 
comportamiento sismico altos (Q=4), ya que aun para acelerogramas asociados a sismos 
postulados intensos, la respuesta no lineal es moderada, estable y concentrada en las vigas. Sin 
embargo, se observa que para este caso, las distorsiones de entrepiso maximas sobrepasan 
notablemente, tanto a las calculadas en su diseno, como al limite maximo de distorsion (A=1.2%) 
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permitido por el RCDF-93. A estos niveles de distorsion se pueden esperar danos considerables en 
elementos no estructurales, asi como problemas de funcionamiento y operation del edificio 
asociados a la conducta humana durante sismos, por lo que debera procurarse ser mas claro y 
explicito en las futuras normas de diseno por sismo y en sus comentarios como deben revisarse 
las restricciones de desplazamientos laterales para fines de diseno en el diseno sismico de 
estructuras irregulares. Las demandas maximas de ductilidad de entrepiso (p) son moderadas, 
ligeramente mayores a las que se presentaron en IR6B e inferiores al factor de comportamiento 
sismico Q asociado a su diseno. 

Con respecto a la pendiente posterior a la fluencia supuesta en las curvas primarias de los 
elementos estructurales (2%, 5% o 10%), no se presentaron grandes diferencias en las respuestas 
dinamicas maximas de los modelos IR6 e IR6B, observandose diferencias mas apreciables para 
los marcos en direction Y que para los marcos en direction X. En general, el aumento en el valor 
de la pendiente posterior a la fluencia tiende a disminuir las demandas de ductilidad, a 
incrementar las rigideces secantes inelasticas, la distorsion de entrepiso asociada a la primera 
fluencia y aumentar la capacidad a cortante de los entrepisos (sobrerresistencia). Habra que 
evaluar en estudios futures la influencia de la pendiente posterior a la fluencia en sistemas 
estructurales que sean sujetos a una respuesta no lineal mas amplia, con el objeto de determinar 
si a estos niveles de excitation se observan o no diferencias mayores en las respuestas dinamicas 
maximas con respecto a las valoradas en IR6 e IR6B. 

Ademas, se deben realizar estudios adicionales donde se evaluen otras fuentes de rigidez 
y sobrerresistencia adicional, como la participation del sistema de piso en la rigidez de las vigas 
considerandolas como vigas de section compuesta, en la respuesta dinamica no lineal y ver de que 
manera afectan la respuesta de estos disenos cuando se sujetan a los acelerogramas empleados en 
este estudio u otros asociados a los espectros de diseno. 

En el presente estudio se proponen algunos parametros que permiten valorar mas 
claramente a la respuesta no lineal de la estructura, ya que en ocasiones los lazos histereticos, el 
mapeo de fluencias, determination de distorsiones maximas de entrepiso, calculo de demandas de 
ductilidad y de rigideces secantes “pico a pico” no son suficientes para comparar la intensidad de 
respuesta no lineal de una estructura cuando se sujeta a movimientos del terreno de distintas 
caracteristicas, o se quiere evaluar el efecto de distintas fuentes de sobrerresistencia. En este 
respecto, el presente estudio ilustra que varios de los parametros dinamicos propuestos resultan 
muy utiles para valorar la importancia y caracteristicas de la respuesta no lineal mas claramente 
que a partir de los lazos de histeresis, y que miden de una manera indirecta la fatiga (de bajo y 
alto ciclaje), entre los que destacan, (l) el cociente entre el numero de ciclos donde se presenta 
respuesta inelastica (Cine) y el numero total de ciclos a que fue sujeta la estructura ante un 
acelerograma dado (Ctot) y, (2) el cociente entre el promedio de las rigideces de entrepiso secantes 
(kiprom) asociadas a los medios ciclos histereticos inelasticos y la rigidez elastica del entrepiso (k e i). 
De hecho, este ultimo parametro resulta mas idoneo para valorar la respuesta inelastica de una 
estructura sujeta a un acelerograma dado que, por ejemplo, el cociente entre la rigidez de 
entrepiso secante minima inelastica (kimin) o “pico a pico” asociada a los ciclos histereticos de 
amplitud maxima y la rigidez elastica del entrepiso (k e i). Este ultimo parametro se utiliza mucho 
en investigaciones experimentales, donde se tiene un mejor control de las amplitudes de 
deformation y de los ciclos a que se sujetan los modelos experimentales, situation que no se da 
cuando una estructura es sujeta a movimientos aleatorios del terreno. 

El primer autor considera que para que los reglamentos evolucionen a criterios de diseno 
por desempeno y por desplazamiento de una manera realista y transparente, se deben conducir 
estudios analiticos encaminados a determinar los niveles de distorsion de entrepiso que deben 
permitirse a edificios, dependiendo del material y sistema estructural utilizado. Resulta 
incongruente tener los mismos limites para todos los sistemas y materiales, independientemente 
de si se trata de estructuras a base de mures de mamposteria o marcos ductiles de acero, como se 
establece en el RDF-93. Existen algunos indices que pueden proponerse a partir de los estudios 
experimentales que se han hecho hasta la fecha en modelos de reducido tamano o escala! de 
hecho, en este tipo de estudios se basan los limites de distorsion que ahora se proponen en el 
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apendice A de las NTCS‘2004, donde ya se diferencian las distorsiones que deben permitirse a 
distintos sistemas estructurales. Sin embargo, muchos de los limites propuestos se basan en 
estudios de sistemas muy simples y, por lo tanto, para aumentar su confiabilidad deben ser 
complementados por modelos analiticos que permitan estudiar un mayor numero de casos y 
variables, utilizando modelos mas complejos que sean representatives de estructuras reales, no 
de modelos simples de laboratorio. 

Por ello, en el presente estudio se determino la distorsion de entrepiso asociada a la 
primera fluencia de cualquier elemento estructural que forma parte del entrepiso (Af e ), ya que 
este indice puede ser muy util para proponer con bases mas firmes, en un futuro, limites de 
distorsion de entrepiso de diseno en los reglamentos que sean congruentes con el material y 
sistema estructural que pretenda utilizarse. Los estudios aqui conducidos muestran una gran 
variabilidad en el valor de la distorsion de entrepiso asociada a la primera fluencia, la que oscilo 
entre Af e =0.28% y Af e =0.97%, obteniendose en algunos entrepisos respuestas elasticas para 
distorsiones tan altas como A=1.07%. A1 parecer, lo anterior puede deberse a un nivel de esfuerzos 
de diseno disparejo en los elementos estructurales en los distintos entrepiso asociado al hecho de 
emplear secciones tipo en el diseno, aspecto que no se ha estudiado con profundidad y que se 
analizara con mayor detalle en estudios futuros. El primer autor ha realizado este mismo tipo de 
analisis en estructuras esbeltas con base en marcos de concreto reforzado, donde se ha observado, 
tanto una variabilidad menor en el valor de Af e , como un valor absoluto menor de Af e con respecto 
al presentado en este estudio para IR6 e IR6B. 

Finalmente, quiere recalcarse que los estudios se realizaron utilizando el programa de 
analisis no lineal DRAIN-2DX, que permite modelar exclusivamente las fallas por flexion y la 
action bidimensional. No se pueden obtener mas que ideas aproximadas de la redistribution de 
esfuerzos tridimensional (como se hizo en los modelos de ocho niveles que no se presentan, pero 
no en IR6 e IR6B), ignorando efectos significativos como es la torsion dinamica no lineal y la 
action multidirectional de las aceleraciones del terreno. De igual manera, no se pueden evaluar 
potenciales fallas por cortante, flexocortante y torsion simultaneamente con las fallas por flexion, 
flexocompresion y flexotension. Se considera importante realizar analisis adicionales en las 
estructuras en estudio tomando en cuenta la action de estos efectos, cuando los recursos de 
computo y de software, permitan realizar este tipo de estudios en modelos tan complejos (aunque 
aparentemente simples) como son los modelos IR6 e IR6B presentados. Para ellos, se requieren 
programas de analisis tridimensionales que modelen con rigurosidad los efectos descritos en el 
intervalo de comportamiento no lineal y en cada paso de integration en el tiempo. Programas 
como el SAP2000 ofrecen un modelado no lineal de baja aproximacion para este tipo de 
estructuras, sobre todo si toma uno en cuenta que la respuesta paso a paso no se calcula con 
rigurosidad en cada paso de integration, sino con un procedimiento “equivalente”, segun se 
documenta en el manual de dicho software. 
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RESUMEN 



Se presents un modelo de generation de acelerogramas artificiales consistente con la 
sismicidad de la region oeste de la Argentina. En su formulation se hace uso del concepto de 
espectro evolutivo, donde la evolution en el campo frecuencial esta materializado por una 
discretization del espectro de Amplitudes de Fourier, definiendo un proceso uniformemente 
modulado en cada franja. 

La modulation en el campo temporal se hace mediante una funcion tipo gamma cuyos 
parametros se determinan en funcion de los momentos temporales del registro, los que a su vez 
son estimados en base a una regresion lineal en funcion de la magnitud y la distancia epicentral; 
de esta manera se concreta la relation con la sismicidad local. 

ABSTRACT 



A model for the artificial generation of accelerograms consistent with the local seismicity 
of the western region of Argentina is presented. It is based upon evolutionary spectrum concepts 
where the evolution in frequency content is taken into account by means of the discretization of 
the Fourier Amplitude Spectra and the definition of a uniform modulated process in each strip. 

The amplitude modulating function is considered to be of gamma type and its parameters 
are determined in function of the temporal moments of the record, which in turn are estimated by 
a multilinear regression analysis in term of magnitude and epicentral distance. In this way the 
relation to local seismicity is considered. 

1. INTRODUCTION 



La Republica Argentina cuenta con una amplia zona adyacente a la Cordillera de los 
Andes con una actividad sismica importante, evidenciada en las cronicas historicas y en los 
ultimos anos, en el catalogo instrumental. Se cuenta con algunos registros obtenidos 
principalmente en la zona cuyana, que es la de mayor actividad sismica, pero ellos son todavia 
insuficientes para las aplicaciones que el calculo estructural exige actualmente. Esta amplia area 
de la region oeste _ central ha sido microzonificada, por lo que el conocimiento que en ella se tiene 
de las condiciones del suelo, tectonica y geomorfologia, es mas acabado que en otras regiones 
argentinas. Se estan realizando tambien grandes esfuerzos para encarar estudios de la 
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sismicidad en la Provincia de Salta, en el norte del pais, que es la segunda zona de importancia 
considerando la peligrosidad sismica, donde no se cuenta con ningun registro de sismos intensos, 
a pesar de que hay evidencias historicas de su ocurrencia. 

Por todo ello es que se considero factible encarar un trabajo de investigation con el 
principal objetivo de definir un modelo de generation de sismos artificiales probables de ocurrir en 
principio, en la region cuyana argentina, para ser aplicados al calculo estructural. Este trabajo 
dio como resultado la realization de dos modelos- uno que reproduce un registro objetivo, con el 
que se realizaron numerosas experiencias numericas a fin de obtener historias de aceleraciones 
concordantes con las de la base de datos disponible, constituida en principio por registros del 
sismo que acontecio en Enero de 1985 en la provincia de Mendoza, cuya sismicidad frecuente se 
debe a movimientos de fallas superficiales cercanas. Basado en este modelo, se propone un 
modelo de generation de registros artificiales, que hace uso de la tecnica de regresion, 
considerando como datos basicos la magnitud y la distancia epicentral que a veces son los unicos 
datos con que puede contarse en algunos sitios, y es el que se explica a continuation. 



2. METODOLOGIA 



Las distintas propuestas y estrategias para encarar la solution de esta problematica, 
pueden englobarse en uno de los dos grupos siguientes ; modelos empiricos y modelos sismologicos. 
Normalmente, los modelos de generation con aplicacion al calculo estructural, que tienen un 
sentido eminentemente ingenieril, siguen la linea de los metodos empiricos. Para ello es necesario 
contar con un buen numero de registros de buena calidad de varios terremotos originados en las 
distintas fuentes sismicas que afecten el sitio en estudio. Se ha sorteado esta dificultad en varios 
paises con pocos registros, usando modelos definidos para otras regiones, y modificando sus 
parametros a condiciones locales, normalmente haciendo regresiones en funcion de magnitud y 
distancia epicentral, que son las magnitudes macrosismicas mas accesibles en todas partes. Los 
modelos sismologicos tienen una aplicacion mas relacionada al ambito de la geofisica y geologia! 
son en general mas complejos y la determination de muchos parametros para definirlos hace 
necesario contar con information y conocimientos sismologicos, que la mayoria de las zonas 
sismicas del tercer mundo no cuentan aun. 

Para este trabajo de investigation se conto con pocos registros, correspondientes a dos 
eventos sismicos acaecidos en las Provincias de Mendoza y San Juan. Se tiene tambien 
information complementaria valiosa, resultante de los estudios de microzonificacion llevados a 
cabo en las areas mas importantes de las mismas. Se penso conjugar en la definition del modelo 
ambas caracteristicas (empiricos y sismologicos); de esta manera se aprovecha la principal 
ventaja de los modelos sismologicos, pero sin perder de vista la aplicacion ingenieril del mismo. 



3. BASES CONCEPTUALES 



El modelo de regresion se basa en el modelo de simulation desarrollado por la autora 
(Orosco, 2003) que hace uso de un acelerograma objetivo a reproducir, por lo que sus bases 
conceptuales se identifican con el mismo. 

A fin de considerar la no estacionariedad de los registros sismicos tanto en el campo del 
tiempo como en el de la frecuencia, el modelo se basa en considerar el movimiento del suelo como 
constituido de varios componentes uniformemente modulados. Esta metodologia fue utilizada 
entre otros, por Der Kiureghian y Crempien (1989^ y Carli y Faravelli (1990^. Tambien en esta 
linea se inscribe la propuesta de Saragoni y Hart (1974) con la diferencia que en el modelo de 
Saragoni se discretiza el movimiento en el dominio del tiempo mientras que en el modelo de 
Crempien la funcion de modulation varia en intervalos determinados en el campo de la 
frecuencia, por lo que presenta continuidad temporal. 
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La aceleracion del suelo debido a un movimiento sismico a{t ) , se representa mediante la 
sumatoria de k procesos uniformemente modulados X ne ^ ( t ) , correspondientes a las k bandas 
consideradas en el campo de las frecuencias. Luego 

a(0 = X *«*(*) (1) 

k 



donde el proceso uniformemente modulado se define como 

XneJ 1 ) = <Pk(t)X e (t) 

por lo que la aceleracion del suelo 

«(0 = Z^Z 

k 

X e ( t ) define en cada banda de frecuencia, un proceso estacionario estadisticamente 

independiente, Gaussiano, de media nula y varianza unitaria, que se expresa mediante la 
conocida expresiom 



( 2 ) 

(3) 



x e (0 = Z 2^S(co)Aco cos (©,.* + ( f >, ) (4) 

i=i 

En la expresion anterior S(ffl) es la funcion de densidad espectral de potencia del proceso 

estacionario en cada franja. Con el objetivo de representar mejor los acelerogramas de la zona 
cuyana argentina, se experimento numericamente con varias funciones de densidad espectral de 
potencia (FDEP), partiendo de una simple funcion constante, utilizada en el modelo de Crempien 
para simular sismos chilenos de subduccion y los del oeste norteamericano (Der Kiureghian y 
Crempien, 1989! Crempien y Araneda, 1993). Cumpliendo con la exigencia de varianza unitaria, 
la FDEP constante se define como: 



S k (a>) = ~ r 

2K + i ~(O k ) 



(5) 



siendo nula fuera de la franja k considerada. 

La funcion espectral adoptada en este trabajo, incorpora en su definition al espectro de 
Boore de la siguiente manera (Orosco, 2003; Queck et al., 1990): 

S 1 (oj) = ^-A:(oj) (6) 

2m 

donde / , denominada intensidad, se define como (Queck et al., 1990) 

I = \\<P(tfdt ( 7 ) 

0 

A a ( CO ) es el espectro de Boore que se expresa como (Boore, 1990) 



ti fl (®) = C5 1 (tu) S 2 (a) S,(co) 



( 8 ) 
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donde C es un factor de escala, S 1 ( CO ) es el espectro de fuente, S 2 ( CO ) tiene en cuenta las 
posibles amplificaciones por efectos locales de sitio y S 3 ( CO ) es el factor de atenuacion. 



En las expresiones (3) y (7) , (p k (/) son funciones de modulacion en funcion del tiempo, 

cuyos parametros se determinan considerando los momentos espectrales del registro, los que se 
estiman para este modelo mediante una regresion multilineal, en funcion de la magnitud y la 
distancia epicentral. Se explica a continuation, detalles del proceso de regresion realizado y la 
base de datos utilizada! en parrafos posteriores se retoma la discusion acerca de la 
determination de la funcion de modulacion. 



4. BASE DE DATOS 



Las tablas 1 y 2 resumen caracteristicas de los sismos utilizados en este estudio. Se los 
dividio en dos grupos, considerando la magnitud. 



Tabla 1 Registros del grupo U1 



Sismo 


M-H 


Registro 


A 


A 

max 


s 


Ia 






(km) 


(cm/seg 2 ) 




(cm 2 /seg 3 ) 


Granada 


M=4 9 


Santa Fe (T) 


40.4 


29.7 


0.67 


0.57 


1984 


H=12 km 


Santa Fe (L) 


40.4 


35.0 


0.67 


0.58 


Beznar (L) 


20.0 


19.80 


0.78 


0.13 


Dursunbey 

1979 


M=5.9 
H=7 km 


Dursunbey 


7.8 


175.2 


0.76 


14.02 






AYLO 


27.2 


98.1 


0.75 


4.17 






AMLO 


28.2 


274.7 


0.57 


36.52 






AMTR 


28.2 


333.5 


0.73 


47.42 






LALO 


50.2 


54.0 


0.61 


12.36 






LATR 


50.2 


59.8 


0.63 


14.28 






LHLO 


33.1 


192.3 


0.76 


17.16 






LHTR 


33.1 


412.0 


0.66 


91.46 






MALO 


20.0 


58.9 


0.65 


1.36 






MATR 


20.0 


186.4 


0.76 


10.03 






POLO 


42.1 


98.1 


0.34 


5.29 






POTR 


42.1 


107.9 


0.34 


5.29 


Mendoza 


M=5.7 


SMLO 


34.3 


78.5 


0.44 


5.78 


1985 


H=12 km 


SMTR 


34.3 


78.5 


0.49 


4.57 






TULO 


55.9 


107.9 


0.43 


9.16 






TUTR 


55.9 


107.9 


0.37 


12.47 






USLO 


77.6 


58.9 


0.40 


3.38 






USTR 


77.6 


107.9 


0.37 


7.13 






ARC 


8.6 


286.3 


0.44 


85.84 






DWY 


16.6 


334.0 


0.62 


94.36 






LWD 


24.0 


291.4 


0.66 


64.31 


Whittier 

Narrows 

1987 




LFS 


31.0 


131.5 


0.52 


13.36 


M=5.9 


PMC 


37.9 


122.1 


0.48 


16.02 


H=14.7 km 


CCY 


51.8 


113.1 


0.61 


11.03 




NWS 


62.0 


85.2 


0.44 


9.10 


Coalinga 














9/5 

Coalinga 

9/7 


M=5.3 


BNT 


12.4 


89.7 


0.74 


3.07 


M=5 


BNT 


12.4 


137.9 


0.58 


12.58 


M=6 


BNT 


12.4 


330.5 


0.57 


90.12 


Coalinga 


22/7 
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Un primer conjunto (llamado Ul) engloba acelerogramas de eventos de magnitud 
comprendida en el intervalo [4.9, 6]; en la tabla 1 se resumen las principales caracteristicas de 
estos eventos. Se detalla la magnitud (M ), profundidad (H) en km, distancia epicentral (A) en 
kilometros, valor pico de la aceleracion (Amax) en cm/seg 2 , factor de Vanmarcke ( S) e Intensidad 
de Arias (Ia) en cm 2 /seg 3 . Un segundo conjunto de datos (al que se denomino U2), ademas de los 
sismos moderados del grupo Ul, incluye registros correspondientes a sismos intensos de mayor 
magnitud, cuyas caracteristicas se muestran en tabla 2. 

5 FORMULAS DE REGRESION 



Evaluando las distintas alternativas y los distintos parametros que pueden ser objeto de 
una regresion, se juzgo conveniente optar por relacionar la energia de los sismos con la distancia 
epicentral y la magnitud. Se cuenta con expresiones de atenuacion de la aceleracion maxima de 
terremotos de corteza superficial para la zona en estudio. Se considero la hipotesis de que la 
energia responde a una ley de comportamiento similar a la de la aceleracion maxima, lo que se 
confirmo en el caso de la base de datos usada. 



Tabla 2 Registros de magnitud mayor a 6 



Sismo 


M-H 


Registro 


A 

(km) 


A 

7 1 max 

(cm/seg 2 ) 


s 


Ia 

(cm 2 /seg 3 ) 


Caucete 1977 


M=7.4 
H= 25 km 


SJLO 

SJTR 


70.0 

70.0 


186.4 

189.3 


0.51 

0.52 


107.35 

133.28 


Loma Prieta 
1989 


M=7 

H=17.6 

km 


Oackland (L) 
Piedmont (T) 
Saratoga (L) 
Saratoga (T) 
Santa Cruz (L) 
Santa Cruz (T) 


95.2 

93.3 
29.8 

29.8 

19.9 
19.9 


306.0 

110.0 
316.2 

494.5 

101.5 
433.1 


0.71 

0.68 

0.62 

0.65 

0.46 

0.45 


81.75 

5.12 

109.40 

146.23 

204.74 

263.73 


Northridge 

1994 


M=6.7 

H=18.4 

km 


Rinaldsi Stn 
(L) 

Rinaldi Stn (T) 

SylmarC.Stn. 

(L) 

SylmarC.Stn 

(T) 

Sylmar Stn E 
(L) 

Sylmar Stn E 
(T) 

Sylmar Hosp 
(L) 

Sylmar Hosp 
(T) 


10.0 

10.0 

12.0 

12.0 

13.0 

13.0 

16.0 
16.0 


804.4 

559.2 
618.0 

696.5 
735.8 
461.1 

863.3 

598.4 


0.81 

0.81 

0.73 

0.74 

0.70 

0.75 

0.73 

0.74 


796.89 

435.87 

425.22 

635.02 

409.93 

319.21 

177.14 

245.44 



5.1 Formula de regresion basada en expresiones de atenuacion existentes 



Se han propuesto varias relaciones de atenuacion para el oeste argentino, que describen la 
variation de la aceleracion pico con la distancia. Bufaliza desarrollo una relation de atenuacion 
usando datos registrados en la zona cuyana argentina, que se expresa como (INPRES, 1995) 



log(a max ) = -0.7837 + 0.353 M - 1.59811og R - 0.000997? 



(9) 
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donde M es la magnitud correspondiente a las ondas superficiales para M > 6 y a las ondas de 
cuerpo para M < 6; R es la distancia hipocentral en km y a ma x es la mayor de las aceleraciones 
pico de las dos componentes horizontales, expresada en fracciones de la aceleracion de la 
gravedad. En el caso de sitios cercanos al piano de falla, la ecuacion (9) deberia ser modificada 
para considerar la saturation del movimiento del terreno, de la siguiente manera 

log(a max ) = -0.7837 + 0.353M - 1 .598 1 log D - 0.00099D ( 10 ) 

donde D se expresa como 



D = tJr 2 + (4.32 io 0073M ) 2 

Se realizo la regresion multilineal mediante la ecuacion 

log( A,. ) = a + bM + c log A + dA 

donde a, b, C y d son los parametros a determinar. Se observa en la ecuacion (12) que son los 
momentos temporales del registro los que se determinan mediante la regresion ya que describen 
de algun modo la energia que contiene y el modo en que esta se distribuye, con lo que la 
information de la genesis, paso y sitio esta contemplada de modo sencillo y general, pero 
suficiente a los fines del calculo estructural. 



( 11 ) 



( 12 ) 



x 10 4 




Figura 1 Relation entre energia y distancia epicentral para sismos grupo U1 
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5.1.1 Analisis de los resultados de la regresion 



En las tablas 3, 4 se resumen los resultados de la regresion para el conjunto de datos 
denominado Ul. En ellas Std son las desviaciones estandar y t es el valor de la prueba de 
Student. Se realizaron ademas las pruebas de significancia de Fischer (Herva et al., 1981) a fin 
de evaluar si al menos uno de los componentes de la formula de regresion utilizada (excepto la 
constante) explican el valor de los parametros que se desean estimar. 



Tabla 3 Resultados de la regresion 1 para grupo Ul 



Parametro 


a 


b 


c 


d 




Std 


Std 


Std 


Std 




t 


T 


t 


t 


A 0 


-9.59 


3.92 


-1.29 


0.013 




3.99 


0.60 


0.94 


0.031 




-2.41 


6.53 


-1.37 


0.42 


A t 


-16.52 


5.05 


-0.53 


0.004 




4.47 


0.67 


1.05 


0.035 




-3.70 


7.50 


-0.50 


0.10 


a 2 


-23.80 


6.24 


0.34 


-0.007 




6.03 


0.91 


4.42 


0.047 




-3.94 


6.86 


0.24 


-0.14 


t f 


-6.92 


1.38 


0.85 


-0.017 




2.35 


0.35 


0.55 


0.018 




-2.94 


3.89 


1.53 


-0.93 


t„ 


-4.28 


0.63 


0.99 


-0.012 




2.37 


0.36 


0.56 


0.018 




-1.80 


1.77 


1.76 


-0.66 



Tabla 4 Resultados de la regresion 1 para el grupo U2 



Parametro 


a 


b 


c 


d 




Std 


Std 


Std 


Std 




t 


T 


t 


t 


An 


0.23 


2.41 


-1.66 


0.010 




2.56 


0.26 


0.74 


0.022 




0.09 


9.17 


-2.24 


0.47 


A , 


-5.39 


3.10 


-0.42 


0.014 




2.78 


0.29 


0.80 


0.024 




-1.94 


10.84 


-0.52 


-0.57 


A 2 


-10.63 


3.73 


0.99 


-0.040 




1.21 


0.37 


1.04 


0.031 




-2.84 


10.14 


0.88 


-1.28 


t f 


-3.43 


0,61 


1.16 


-0.025 


J 


1.24 


0,12 


0.35 


0.011 




-2.84 


4.90 


3.31 


-2.34 


t e 


-4.41 


0.51 


1.36 


-0.027 




1.24 


0.13 


0.36 


0.011 




-3.54 


3.02 


3.79 


-2.52 



Analizando los valores para el grupo Ul de datos, se tiene que el momento /£ 0 tiene un 

valor F de 16.07 que es mucho mayor que F(0.0l; 3, 27)=4.60, donde 0.01 es la densidad de 
probabilidad, 3 es el numero de parametros a determinar (sin considerar la constante) y 27 son 
los grados de libertad (cantidad de datos de la muestra menos numero total de parametros a 
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estimar). Se procedio del mismo modo con el resto de los parametros, obteniendose el valor mas 
bajo de F para la duracion efectiva (6.12). 

Acerca de la prueba de Student una evaluacion general de los parametros a,b,C y d 
considerando t (0.10; 27) =1.70 lleva a la conclusion que el termino con menor significancia en 
todos los casos es d , ya que para el mismo, el valor mas alto obtenido fue de 1 1 1 =0.93. En cuanto 
a la constante C , los valores correspondientes a A ] y X 2 son malos! para los restantes 

parametros, no obstante sus bajos valores, se pueden aceptar. En general, la constante b explica 
bien el modelo propuesto en todos los casos. De este analisis puede inferirse que la distancia 
epicentral no explica muy bien el modelo para esta base de datos. 

Considerando ahora el grupo de datos U2 los valores de F para los momentos espectrales 
se ubican entre 35 y 45 mientras que para la duracion total y efectiva, F es aproximadamente 12; 
considerando el nivel F(0.01! 3, 43) = 4.31, se concluye que en el modelo al menos alguna de sus 
variables es significativa para explicar los valores encontrados. 

En cuanto al valor de las constantes a,b,C y d, se observa que sus valores t 
aumentaron y teniendo en cuenta un mejor nivel de estimacion (t(0.02; 43)= 2.42) permite inferir 
que el termino que multiplica a d explica la duracion total y efectiva pero no asi los momentos 
temporales. En cuanto a la constante C , tal como en el caso anterior no explica los momentos de 
orden 1 y 2 pero si explica, desde el punto de vista estadistico, el valor de la energia, la duracion 
total y la efectiva. El valor de b es muy bueno en cuanto a su significancia en el modelo de 
regresion propuesto. Como conclusion, la magnitud y el logaritmo de la distancia epicentral 
explican muy bien el valor de la energia. Para la duracion efectiva y total todos los parametros 
son significativos, mientras que para los momentos Aj y A 1 , solo la magnitud tiene relevancia. 



5.2 Modelo de regresion propuesto 



En vista del analisis estadistico realizado a los resultados de la regresion hecha en base a 
las formulas de atenuacion existentes, se propone una expresion que determina la energia del 
proceso y la duracion efectiva en funcion de dos parametros mas la constante, dejando de lado la 
distancia epicentral (termino que multiplica la constante d ) por los bajos valores de significancia 
obtenidos en el caso anterior. La nueva formula de regresion tiene la expresiom 

log ;(A 0 ) = a + b log(M) + c log(A) 
log(( 0 ) = a + Mog(M) + clog(A) 

Con el modelo anterior, en la estimacion de los valores de los momentos temporales de 
orden 1 y 2 solo la magnitud tiene relevancia! por ello la formula de regresion para los momentos 
de orden superior se expresa en funcion de A 0 y t e de la siguiente manera- 

l°g(^ 7 )- a + b log(/l 0 ) + c log(( e ) ( = 1,2 ( 14 ) 



donde A 0 y t e son determinados mediante (13). 

Los resultados de la regresion efectuada para los dos grupos de datos se resumen en la 

tabla 4. 
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Tabla 4 Resultados de la regresion con el modelo alternativo 



Parametros 


Grupo Ul 


Grupo U2 




a 


b 


c 


a 


b 


c 




Std 


Std 


Std 


Std 


Std 


Std 




t 


t 


t 


t 


t 


t 


K 


-25.02 


21.25 


-0.91 


-14.44 


15.68 


-1.20 


5.49 


3.19 


0.29 


3.03 


1.54 


0.25 




-4.56 


6.65 


-3.13 


-4.76 


10.18 


-4.73 


A, 


-1.370 


1.125 


0.897 


-0.982 


1.063 


0.953 




0.174 


0.017 


0.037 


0.158 


0.013 


0.040 




-7.89 


65.90 


24.30 


-6.20 


79.46 


24.13 


a 2 


-2.853 


1.259 


1.942 


-1.788 


1.118 


1.975 




0.335 


0.033 


0.071 


0.298 


0.025 


0.074 




-8.51 


38.18 


27.24 


-6.00 


44.40 


26.56 




-5.81 


3.40 


0.64 


-6.42 


3.63 


0.69 




3.28 


1.91 


0.17 


1.61 


0.82 


0.13 




-1.77 


1.78 


3.65 


-3.98 


4.43 


5.11 



5.2.1 Analisis de los resultados del modelo alternativo de regresion 



Para el conjunto de datos Ul, los valores de F para el momento de orden 0 y la duration 
efectiva son de 24.81 y 9.18 respectivamente, lo que considerando un nivel F(0.0l! 2, 28) = 5.45 se 
concluye que las variables del modelo tienen una alta significancia estadistica para explicar los 
parametros mencionados. En cuanto a 1, y 1 2 sus valores son 2,387 y 1,047 respectivamente, lo 
que estadisticamente significa un valor muy bueno de relation entre los parametros involucrados. 

Para el grupo de datos U2 los valores de F se incrementaron con respecto al caso anterior! 
se tiene para los momentos temporales las siguientes cantidades^ 77, 3,489 y 1,363 
respectivamente, mientras que para la duration efectiva es 19 y considerando un nivel F(0.0l! 2, 
42) = 5.18, se concluye que hay una relevante relation entre los parametros estimados y las 
constantes obtenidas por la regresion. 

Analizando a,b y c con la prueba de Student para los datos del grupo Ul, los valores t 
de todas las constantes sobrepasan el nivel t(0.00i; 28)=3.66 excepto el caso de C (que solo esta 
a 0.005 por debajo de este valor) y las constantes a y b de la duration efectiva (que satisfacen el 
nivel t(0.10; 28)=1.70). En cuanto al grupo U2, todas las variables sobrepasan el nivel 
t(0.001!42)=3.55, lo que implica que el modelo es muy bueno desde el punto de vista estadistico. 
Asi, el logaritmo de la energia y de la duration efectiva explican muy bien los momentos 
temporales de orden 1 y 2, lo que puede observarse del analisis de sismos reales. 

Como las regresiones se realizan sobre los logaritmos de la magnitud y distancia 
epicentral, pequenas variaciones en ellos provocan grandes diferencias en los resultados, por ello 
la regresion debe ser realizada con cuidado para obtener resultados confiables. 

Este modelo se denomina de “regresion 2” para diferenciarlo del anterior realizado en 
base a la formula (12) al que se identifica como modelo de “regresion 1”. Se hicieron simulaciones 
con ambas expresiones y en general se observa una mayor estabilidad en las formas de las 
funciones de amplitud determinadas en base a los momentos temporales obtenidos por el modelo 
de regresion 2 con respecto a las obtenidas por la regresion 1, en especial para registros de banda 
angosta que en general corresponde a campo cercano, mientras que las diferencias son 
insignificantes para el caso de eventos de banda ancha. 
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6. PARAMETROS DE CONTROL DEL MODELO 

Los aspectos importantes a considerar para determinar los parametros que controlan el 
proceso de generation de historias con determinadas caracteristicas se explicitan a continuation. 
Se incluyen algunos ejemplos de aplicacion a fin de mostrar sus efectos e influencias en los 
registros obtenidos. 

a) Tipo de registro que se desea obteneti banda angosta o banda ancha. 

b) Distribution de la energia en cada franja frecuencial. 

c) Diseno de la funcion de amplitud. 

d) Diseno de la funcion de densidad espectral de potencia. 

e) Control de la energia en cada franja. 

f) Duration efectiva y total, y aceleracion maxima. 



6.1 Tipo de registro 



Es posible decidir si se desea obtener un acelerograma de banda angosta o de banda 
ancha. Tales caracteristicas son determinadas por el valor del factor de (Vanmarcke, 1980) del 
registro que se define como- 



s= 




JL 

Vi 



donde los A- son los momentos espectrales del registro a(t) definidos por 



(15) 



„ rco, 

A t = £ co' A f (co) i = 0, 1, 2 ( 16 ) 

y A f (C0) es el espectro de amplitudes de Fourier. Para valores de 8 mayores a 0.66 se tienen 

registros de banda angosta y para valores menores a 0.66, registros de banda ancha. La base de 
datos de Mendoza arrojo como conclusion que para registros de banda angosta se obtenian 
mejores resultados en el proceso de simulation considerando una division de tres franjas en el 
campo frecuencial, mientras que para los de banda ancha la division en cinco franjas era la 
apropiada. Si bien esto puede limitar el universo a sismos similares solo a los del evento con que 
se cuenta mayor cantidad de datos, sigue siendo valido el concepto en el que se basa la division 
del espectro- considerar “paquetes de energia”, esto es dividir el espectro de energia del registro 
en tantas franjas como queden determinadas por los valores nulos o casi nulos del espectro de 
energia (considerese que esto es importante en el caso de utilizar el modelo que reproduce las 
caracteristicas de un sismo muestra). 

Es entonces posible “crear” un registro con un numero de paquetes de energia distinto a 
tres y cinco, si asi se desea. Se observo en el proceso de simular algunos sismos distintos a los 
cuyanos que para sismos muy impulsivos la division en dos bandas podria funcionar bien. 
Mientras que en sismos de banda muy ancha, con bajos valores del factor de Vanmarcke (menores 
de 0.45 por ejemplo), a veces es preciso recurrir a mas de 5 divisiones en el campo de la 
frecuencia. El sismo de Caucete (San Juan) es un ejemplo. 
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En base a los datos con que se conto para este estudio, para distancia epicentrales 
mayores a 40 km se recomienda considerar historias de banda ancha, para distancias epicentrales 
menores de 30 km, se aconseja considerar banda estrecha, y para los situados entre 30 km y 40 
km se adopta uno u otro criterio de acuerdo al que mejor se adapte a los resultados que se desean 
obtener. Esto es solo una sugerencia, ya que se observo en la base de datos algunas variaciones a 
la tendencia general, pues las caracteristicas de los movimientos en el sitio de estudio dependen 
grandemente de las condiciones locales de suelo. 



6.2 Distribution de la energia en cada banda de frecuencia 



De todos los analisis numericos realizados se calculo el promedio del monto de la energia 
contenido en cada banda frecuencial, tanto para registros de banda angosta como de banda ancha. 
Esto no impide decidir “a priori” los porcentajes de la energia total con que se dotara a cada banda 
de frecuencia considerada! asi puede tenerse el caso de que se requiera probar el comportamiento 
dinamico no lineal de una estructura sujeta a una historia de aceleraciones en un rango de 
frecuencia determinado (cercano al periodo propio, por ejemplo), y al no contar con registros reales 
con suficiente energia en tal rango, es preciso simular historias en el tiempo que respondan a esa 
caracteristica deseada. 

Para sismos de banda estrecha (para los que se requiere dividir el espectro en tres 
franjas), los valores medios de la fraction de energia total obtenidos para cada banda son: 0.77, 
0.20 y 0.05 con desviaciones estandar igual a 0.13, 0.13 y 0.04 respectivamente. Para sismos de 
banda ancha (con division en cinco franjas), se tiene que el promedio de las fracciones de la 
energia total en cada franja son 0.02, 0.24, 0.52, 0.20 y 0.03 respectivamente, siendo sus 
desviaciones estandar 0.02, 0.18, 0.28, 0.29 y 0.03. Se hace notar nuevamente que estos valores 
son solo indicatives. 



6.3 Diseno de la funcion de amplitud 



A los fines de modular las amplitudes de los acelerogramas simulados, varios 
investigadores hicieron distintas propuestas de funciones de modulation a lo largo del tiempo 
(Barbat et al, 1994). La mas simple es considerar un pulso rectangular, de duration l , . Esta 

funcion no describe apropiadamente las variaciones en el valor cuadratico medio de los 
acelerogramas, que es la caracteristica de la no estacionariedad en el tiempo. 

Se probaron varias funciones, a fin de reproducir del mejor modo posible los registros 
muestras. Los registros disponibles en la base de datos tendian a una funcion de amplitud tipo 
gamma, por lo que fue esta la adoptada para el modelo. 

Se define la funcion gamma de la siguiente manera : 

<pit m ) = at m p e-* (17) 

en la que la variable t m —t — t 0 es la variable tiempo modificada por un parametro t 0 que tiene en 

cuenta el instante en que irrumpe el pulso, importante a tener en cuanta en el caso de sismos 
impulsivos, como lo eran la mayoria de la base de datos mendocina. En lo que sigue, se obviara 
escribirlo de esta manera, pero siempre debe tenerse en cuenta este hecho. 
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Tiempo 



Figura 2 Funcion de amplitud para registros de banda ancha (Casos A, B, C y D) 



A fin de estimar los parametros OC , (3 y J de las funciones de modulation (p(t ) en cada 
franja, se impone la condition de ergodicidad de los procesos no estacionarios (Bogdanoff y Kozin, 
1961) 

tf tj 

/L„= ft” <p 2 (t)dt = ft n a 2 ( t ) dt (is) 

0 0 



con lo que se obtienen las siguientes expresiones algebraicas, en las que los parametros (X , (3 y 
J estan en funcion de los momentos temporales obtenidos por regresion, 



a = 



P = 



7 = 



K Y P+ 1 

r(/? + i) 

% 

3, 0 A 2 — 
A 0 A 2 



-1 



(19) 



Los resultados de aplicar esta metodologia para un evento de banda ancha, de 32 km de 
distancia epicentral y una magnitud M = 5.4 grados se muestran como Caso A en la figura 2. 
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Figura 3 Registros artificiales de banda ancha (casos A, B, C y D) 



Otra manera de estimar la forma de la funcion de amplitud es mediante relaciones entre 
los momentos de orden 1 y 2 (que justamente describen el modo en que la energia se distribuye) 
con respecto al momento de orden 0 (que describe la energia contenida en rango frecuencial 
considerado) y cuyo valor es el determinado por la tecnica de regresion. 



Se definen entonces los parametros k x y k 2 como 



k \ 

k 2 



A 

K 

K 



( 20 ) 



En la Tabla 5 se muestran los valores medios de los parametros k x , k 2 y t Q para sismos 

cuyanos, obtenidos de los resultados del ajuste no lineal efectuado en cada franja con el modelo de 
simulation que usa el sismo muestra como semilla de generation. Como se aprecia el rango de 
valores en que tales parametros se ubican es amplio, como lo denotan los altos valores de sus 
desviaciones estandar (que es el valor que aparece entre parentesis). 

Los valores de k x y k 1 son fijados a voluntad, segun la apariencia de la funcion de 
modulation que se desea lograr en cada franja frecuencial. Por ello, los valores medios se 
muestran a titulo de ejemplo. A fin de asegurar que la modulation del proceso uniforme en cada 
franja sea mediante la funcion gamma, los valores de k x , k 2 y t 0 deben satisfacer los siguientes 
requisitos, a fin de que los parametros de la funcion no resulten negativos, en especial el valor de 
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^ <k x <^k ~ 2 



( 21 ) 



■ k x - ^3k x - k 2 < t 0 < -k x + ^3k x - k 2 



( 22 ) 



Tabla 5 Valores medios de k x , k 2 y t 0 (sismos mendocinos) 



Franja 


Banda estrecha 


Banda ancha 


K 


k 2 


^0 


k i 


k 2 


^0 


1 


2.51 

(0.59) 


9.04 

(3.24) 


1.61 

(0.78) 


6.84 

(3.02) 


71.82 

(59.61) 


2.49 

(2.24) 


2 


2.83 

(0.66) 


11.60 

(3.46) 


1.62 

(1.01) 


12.71 

(8.38) 


356.53 

(360.93) 


5.49 

(3.84) 


3 


2.88 

(0.59) 


11.53 

(3.78) 


1.87 

(1.01) 


10.04 

(5.25) 


195.09 

(183.03) 


4.40 

(2.70) 


4 








10.14 

(3.18) 


185.87 

(115.96) 


5.03 

(2.46) 


5 








11.84 

(4.26) 


287.58 

(224.81) 


5.04 

(2.50) 



La Tabla 6 muestra algunos ejemplos mas de pares de valores k x y k 2 que se agrupan 
en los denominados casos B, C y D y cuyas las funciones resultantes se aprecian en figura 2! las 
historias correspondientes y la evolucion temporal de la energia en cada ejemplo, se muestran en 
las figuras 3 y 4 respectivamente. 



Tabla 6 Ejemplos de valores de k x y k-, para sismos de banda ancha 



Franja 


Caso B 


Caso C 


Caso D 




K 


k 2 


k \ 


k 2 


K 


k 2 


1 


5.00 


43.17 


6.00 


59.51 


10.33 


140.24 


2 


4.00 


23.07 


15.00 


356.93 


7.95 


77.72 


3 


5.00 


36.52 


12.00 


195.09 


10.12 


135.00 


4 


5.00 


32.82 


11.00 


185.87 


8.79 


99.47 


5 


6.00 


52.70 


15.00 


287.58 


10.96 


168.38 



En todos los casos mostrados en la figura 2, se considero un valor nulo para el parametro 
t 0 . A fin de graficar su influencia en la forma de la funcion de amplitud buscada, se muestra en la 
figura 5 un detalle de las funciones de modulacion de la segunda y tercera franja correspondientes 
al caso D, en las que se vario el valor de t 0 . El intervalo de valores de este parametro adecuados 

para determinar una funcion gamma es reducido! para el caso denominado Dl, t 0 es de 1 segundo 

en ambas franjas! en el caso D2, los valores son de 2.5 y 2.80 segundos respectivamente y el caso 
D3 corresponde a valores de “1 segundo en cada una de ellas, manteniendo todas las demas 
condiciones constantes con respecto al caso D. La figura 6 muestra las historias resultantes. 

Por otro lado, en la figura 7 se observa el efecto de variar k x en la segunda y tercer franja 
frecuenciales. El caso denominado D4 corresponde a valores de k x cercanos a los limites inferiores 
del intervalo definido por (21) mientras que el caso D5 resulta de considerar valores muy 
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cercanos al llmite superior. La apariencia de los acelerogramas simulados se muestran en la 
figura 8. 




Figura 4 Intensidad de Arias para casos A, B, C y D. 




Tiempo 



Figura 5 Efecto de variar t 0 en la funcion envolvente de la segunda y tercera franja 



128 



Lia Haarala Orosco, Mika Haarala Orosco, Alex Barbat 



100 



42 0 



-100 




CASO D1 






100 



10 



15 20 



25 



30 



35 



4? 0 



-100 




CASO D2 






100 



4? 0 



-100 



10 



15 20 



25 



30 



35 



| CASO 03 

|l| I 




10 



15 20 

Tiempo (seg) 



25 



30 



35 



Figura 6 Registros artificiales de banda ancha (casos Dl, D2 y D3) 




Figura 7 Efecto de variar k x en la funcion envolvente de la segunda y tercera franja 
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Figura 8 Registros artificiales de banda ancha (casos D4 y D5) 



En el caso que se consideren valores de t 0 que no cumplan con la condition (22) pueden 

obtenerse funciones envolventes del tipo que se muestra en la figura 9, para las dos franjas 
mencionadas. Si el valor del parametro es mayor que el limite superior, la funcion se hace muy 
estrecha en la punta y deja de ser del tipo gamma! la historia simulada con este caso se muestra 
en figura 10. 




Figura 9 Funcion envolvente obtenida cuando t 0 no cumple la condition (22) 
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Analizando el caso de registros de banda angosta, en la Tabla 5 se muestran los valores 
medios y las desviaciones estandar de los parametros k l , k 2 y t 0 , mientras que en la tabla 7 se 

anotan los valores correspondientes a tres combinaciones de valores de las relaciones entre los 
momentos espectrales con las que se obtuvieron las funciones de modulation de amplitud que se 
detallan en la figura 11. Este ejemplo corresponde a un registro de magnitud 5.2 y una distancia 
epicentral de 15 km. Para la distribution de energia entre las franjas se adoptaron valores 
cercanos a los medios, esto es 70%, 20% y 10 % respectivamente. La figura 12 muestra los 
acelerogramas simulados segun las funciones de modulation detalladas en la figura 11. Para el 
caso D los resultados son similares a los obtenidos con la regresion 1, por lo que se muestra solo 
uno de ellos. 

Las funciones obtenidas de aplicar en las formulas (19) los momentos temporales 
estimados directamente por regresion son las denominadas “Regresion 1” y “Regresion 2” en la 
figura 11. 



Tabla 7 Ejemplos de valores de k { y k 2 para sismos de banda angosta 



Franja 


Caso B 


Caso C 


Caso D 


ki 


k 2 


K 


k 2 


ki 


k 2 


1 


2.55 


9.06 


3.00 


14.50 


2.55 


7.00 


2 


2.90 


12.08 


3.00 


16.00 


2.80 


8.00 


3 


2.92 


11.59 


3.00 


17.00 


2.91 


8.00 




Los resultados obtenidos segun una u otra formula de regresion difieren notoriamente en 
este caso, lo que se detalla en las figuras 13 y 14, en las que se designa como regresion 1 las 
curvas envolventes obtenidas por (12) y como regresion 2 las calculadas segun (13) y (14). En el 
curso de las experimentaciones numericas llevadas a cabo durante el desarrollo de este trabajo, se 
observo que las mayores diferencias entre ambos modelos de regresion se dan al considerar como 
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base de datos el conjunto denominado U1 y en especial para registros de banda angosta y de 
campo cercano! por otro lado, al realizar la regresion sobre el grupo U2 las funciones resultantes 
tienden a ganar estabilidad en su forma y a ser similares, tanto para registros de banda angosta 
como ancha. La figura 13 muestra un ejemplo de funciones envolventes para un sismo de banda 
angosta de magnitud 5, distancia epicentral de 10 km y profundidad de foco de 15 km mientras 
que la figura 14 hace lo propio para un sismo de magnitud 5.9, distancia epicentral de 30 km y 
profundidad de foco de 20 km. 




Figura 11 Funciones de modulation para registros de banda angosta 
(Casos A, B, C y D) 
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Figura 13 Diferencias entre las funciones de amplitud resultantes de 
considerar U1 o U2 (M=5, A = 10 km, H=15 km) 
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Figura 14 Diferencias entre las funciones de amplitud resultantes de 
considerar U1 o U2 (M=5.9, A = 30 km, H= 20 km) 

Se aclara que en cada ejemplo de este paragrafo se mantuvieron constantes todas las 
variables del problema salvo aquellas cuyo efecto se queria demostrar. En cada caso el proceso 
estacionario es el mismo, pues se mantuvieron constante los nucleos del generador de numeros 
aleatorios. 
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6.4 Diseno de la funcion de densidad espectral de potencia 



La generation de acelerogramas artificiales con este modelo de regresion puede realizarse 
utilizando como funcion espectral de potencia una funcion constante dada por la expresion (5) o la 
que incluye el espectro de fuente de Boore, expresada por (6). En general puede decirse que si se 
desean registros de banda angosta, como los sismos de tipo impulsivo por ejemplo, la funcion que 
incluye el espectro de fuente da mejores resultados, mientras que para los sismos de banda ancha, 
una distribution constante para la densidad espectral es suficiente. 



6.5 Control de la energia 



La energia total del registro se determina por regresion lineal, en funcion de la magnitud 
y la distancia epicentral. Por medio de la division porcentual de esa energia en las franjas 
consideradas segun el tipo de registro que se desea, se determina la energia de cada banda y se 
disena la funcion de amplitud correspondiente en base a lo expuesto anteriormente. Una vez que 
el proceso de generation es llevado a cabo banda a banda, la energia en cada una de ellas es 
controlada mediante la potencia (energia por unidad de tiempo) de la funcion de amplitud y la de 
la historia simulada correspondiente a esa banda frecuencial. Una vez realizada la correction si 
correspondiere, se suma la contribution de esta banda al registro total que se genera. Asi se logra 
que la energia del proceso generado es la que se obtuvo por regresion lineal. 

6.6 Duration total y efectiva del registro generado 

Con respecto a la estimation de la duration del registro que se busca generar, se pueden 
considerar distintos criterios. Asi por ejemplo, Crempien y Araneda (1993) hacen una regresion 
lineal sobre la duration de la fase fuerte del movimiento, con la misma formula de regresion 
usada para los momentos. Hay varios aspectos que influyen en la duration del registro en un 
determinado sitio, la que depende de factores como la frecuencia y la geomorfologia del lugar 
entre otros (Trifunac y Novikova, 1994) 

En este trabajo estos parametros tambien se determinan por regresion lineal, 
considerando las dos formulas propuestas. Los valores de las constantes obtenidos por medio de 
la regresion 1 se muestran en las tablas 3a y 3b para ambos grupos de datos. Por otro lado, en la 
tabla 4 se detallan los resultados de la regresion realizada segun las formulas del modelo 
alternativo de regresion. 

Considerando los datos con que se cuenta, se juzga que esta simple metodologia dio 
resultados aceptables. De todos modos, es posible dar valores distintos a los calculados inferidos 
por necesidades de la aplicacion del registro generado (por ejemplo, para probar casos de 
licuefaccion de suelos se puede exigir sismos de mayor duration total y efectiva), o por condiciones 
particulares del suelo y la topografia en el sitio de estudio. 

Una vez definida la duration total del evento a simular, a los fines del algoritmo de 
simulation se debe definir el paso de discretization de la variable tiempo para determinar el 
numero de puntos del acelerograma artificial. La funcion temporal de amplitud no depende del 
intervalo de discretization At pero si hay diferencias en el campo frecuencial. Cuanto mas grande 
sea este paso de tiempo, el contenido espectral es menor y si la duration total del acelerograma 
generado es pequena se perdera information. 

Como ejemplo, para simular acelerogramas con una duration de 5 a 6 segundos, un paso 
de discretization temporal de 0.005 segundos dio buenos resultados! en el caso de registros de 10 
segundos de duration total, un paso de 0.01 segundos dio resultados satisfactorios. Puede 
sugerirse que si la duration total es mayor a 20 segundos, un intervalo de tiempo de 0.02 
segundos es lo aconsejable. 
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6. 7 Aceleracion maxima del registro generado 

Finalmente, con respecto a la aceleracion maxima deseada, esta puede ser definida en 
base a la obtenida por alguna formula de regresion, por condiciones geomorfologicas, exigencias 
de algun reglamento, etc. En general por este metodo se obtienen aceleraciones cercanas a las 
exigidas, dependiendo del buen diseno de la funcion de amplitud principalmente! se admite una 
ligera correction del pico, siempre y cuando esto no implique deformaciones inadmisibles de la 
forma del registro. 

7. EJEMPLO DE APLICACION 



Con los conceptos antes expuestos se ha confeccionado un programa de computadora 
(PREDICE) con el que se ha realizado una amplia experimentation numerica aplicada 
principalmente a la zona en estudio, dando resultados satisfactorios tanto para sismos moderados 
como intensos. Se muestra un ejemplo que tiene como objetivo simular acelerogramas de un 
evento similar al ocurrido el 24 de Junio de 1984 en Granada (Espana). No es el principal 
proposito de este modelo generar acelerogramas concordantes con alguno tornado como muestra, 
pero se desea analizar los alcances de la aplicacion del modelo y ver de que manera es posible 
controlar sus parametros a fin de obtener historias con caracteristicas determinadas a priori. El 
espectro de Boore fue estimado con los siguientes valores ( Morales et al., 1996): 



Magnitude 5 

Distancia epicentrab 40.4 km 
Profundidad: 12 km 
Momento sismico: 3.0 10 16 N-m 
Caida de tension: 0.100 10 8 N/m 2 
Frecuencia maxima: 25 Hz 

Coeficiente de atenuacion dependiente de la frecuencia: 80 f 3 - 65 

La historia temporal de las aceleraciones registradas en Santa Fe corresponde a un 
evento de banda angosta por lo que se procedio a dividir el espectro en tres franjas frecuenciales, 
con un contenido energetico en cada franja de 85%, 10% y 5% respectivamente. En cuanto a las 
funciones de amplitud que modula el proceso estacionario en cada franja dieron buenos resultados 
las disenadas tanto en base a parametros determinados por regresion, como para relaciones entre 
momentos espectrales correspondientes al caso B de la tabla 7 y valores del parametro t 0 de 
0.06, 0.05 y 0.10 segundos para cada franja respectivamente. 




Figura 15 Registros SFLO y SIM1 
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En las figuras 15 a 19 se observan los resultados obtenidos para dos casos de seis 
simulaciones efectuadas. Se muestran las historias simuladas y sus respectivos espectros de 
amplitudes de Fourier, cotejandolas con el registro y espectro de Fourier originales. Las distintas 
simulaciones resultan de considerar distintos nucleos para el calculo de los procesos aleatorios 
correspondientes. La figura 19 muestra los espectros de respuesta elasticos para otras tres 
simulaciones y puede observarse una buena correspondencia con el del registro original. 




Figura 16 Espectros de Fourier de SFLO y SIM1 




Tiempo (seg) 



Figura 17 Comparacion entre SFLO y SIM5 
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Figura 18 Espectros de Fourier de SFLO y SIM5 




Figura 19 Espectros de respuesta elastica de SFLO y tres simulaciones 



8. CONCLUSIONES 



En cuento a las bases conceptuales sobre las que se define este modelo, se contempla la 
no estacionariedad temporal y frecuencial, mediante la aplicacion del concepto de espectro 
evolutivo donde la evolucion del contenido frecuencial en el tiempo se materializa mediante una 
discretizacion del proceso en el campo frecuencial, considerando luego en el campo del tiempo a 
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cada franja componente como un proceso estocastico uniformemente modulado, metodologia 
utilizada previamente por varios autores y aplicada a zonas que cuentan al menos con un registro. 

La determination de los parametros de la funcion de amplitud que modula cada proceso 
estacionario en las distintas franjas se los determina en funcion de los momentos temporales del 
registro a simular, obtenidos por medio de la tecnica de regresion multilineal en funcion de dos 
magnitudes macrosismicas de uso corriente en sismologia^ magnitud y distancia epicentral. 

En cuanto a la base de datos utilizada en la regresion consta de registros cuyas 
magnitudes se ubican en el intervalo [4,7-4] los que fueron divididos en dos grupos- 

• El primero, al que se denomina Ul, contiene eventos de magnitudes comprendidas entre 4 y 6 
grados relacionadas a la frecuente actividad sismica moderada de la zona cuyana. Amen de 
los registros del sismo de Mendoza'85 se incluyen otros considerados del mismo tipo que los 
del evento mencionado, esto es sismos moderados, superficiales y resultantes de 
deslizamientos de fallas debido mayormente a un proceso tensional de compresion. 

• Para incluir los dos acelerogramas del sismo de Caucete, de magnitud 7.4, se anadieron 
algunos acelerogramas de sismos superficiales intensos del Oeste Americano, con lo que se 
conforma el conjunto U2 de datos. Ademas, teniendo en consideration que la ciudad de 
Mendoza y sus alrededores esta asentada sobre varias fallas activas muy superficiales, se 
juzga factible que ocurran sismos de este tipo en la region bajo estudio, habida cuenta de que 
los informes de microzonificacion realizadas en esta area hablan de similitudes geotectonicas 
entre ambas regiones. 

• Se comprobo que tanto la aceleracion maxima como la energia de los registros de la base de 
datos utilizada siguen una distribution similar con respecto a la distancia epicentral y la 
magnitud del sismo. Por ello en primer lugar se hizo la regresion multilineal en base a la 
formula de atenuacion de aceleraciones maximas determinada por profesionales del Instituto 
National de Prevention Sismica de la Argentina, con la que se obtiene los valores de A 0 , y 

A-, . La misma expresion fue aplicada para determinar las duraciones efectiva t e y duration 
total 1 1 de los registros a simular, obteniendose buenos resultados. 



• Analizando las propiedades estadisticas de los resultados de la regresion se concluye que con 
esta formula la constante que multiplica a la distancia epicentral tiene en general poca 
significancia para explicar una relation con los parametros que se desean estimar! se infiere 
que en el caso de los momentos de orden 1 y 2 y en especial para el conjunto de datos Ul solo 
la magnitud tiene relevancia en la definition de estos parametros. 

• A fin de mejorar estos aspectos se propone una segunda formula de regresion o mas 
propiamente dos formulas! primeramente se realiza una regresion para determinar el 
momento de orden 0 (energia del registro) y las duraciones efectiva y total, los parametros 
que estadisticamente tanto la magnitud como la distancia epicentral explican 
significativamente con las formulas de atenuacion. Los momentos de orden 1 y 2 se 
determinan tambien por regresion, mediante una expresion que es funcion de la energia y la 
duration efectiva determinadas por regresion. 

• Con estas expresiones el modelo (denominado de “regresion 2”) presenta mejores niveles en 
las pruebas de confidencia y significancia estadisticas que el modelo basado en las formulas 
de atenuacion existentes para la zona bajo estudio (identificado como modelo de “regresion 
1 ”). 



• Se realizaron simulaciones utilizando ambos conjuntos de datos y con ambas formulas de 
regresion. Se detectaron diferencias entre las funciones de amplitud logradas con uno u otro 
metodo al utilizar el grupo Ul de datos, sobre todo para registros de campo cercano. Al 
ampliar la base de datos y hacer uso del grupo U2, se observa mayor estabilidad en la formas 
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de las funciones como asi tambien una reduction apreciable de las diferencias en los 
resultados segun ambos metodos, tanto para eventos de campo cercano como lejano. 

• Las regresiones se realizaron sobre los datos sin hacer especiales correcciones ni filtrados de 
ruidos o fallas de muestreo. Por ello, los momentos de orden 1 y 2 resultantes conllevan estos 
efectos que como se explico antes, dependen del intervalo considerado antes del primer 
impulso de la serial, esto es del parametro t 0 y el ruido de la serial (parametro S). 

• En este caso, el modelo no cuantifica automaticamente estos parametros por lo que las 
funciones de amplitud resultantes pueden ser defectuosas y no describir apropiadamente la 
variation temporal de la serial. Se determinaron formulas explicitas para corregir los 
momentos temporales de orden 1 y 2 a fin de incluir estos efectos de manera simple y clara en 
la determination de los mismos. 

• Esta correction se impone en caso de generar sismos impulsivos de banda angosta y no es tan 
critica en el caso de eventos de banda ancha. 

Una vez determinados los momentos temporales totales, se calculan los mismos para cada 
franja como fracciones de los respectivos momentos totales, considerandose los mismos 
porcentajes del contenido energetico de la franja con respecto a la energia total del evento. 

Cabe tambien en el modelo la posibilidad de manipular los momentos de orden 1 y 2 a fin 
de permitir la construction de historias probables de ocurrir en un lugar determinado, con 
caracteristicas especificas de acuerdo a la aplicacion que se hara de las historias simuladas, 
manteniendo el momento de orden 0 (que representa la energia de la franja) inalterable. 

• Se proponen a tal fin los dos parametros k x y k 2 , que simplemente relacionan los momentos 
espectrales de orden 1 y 2 con el momento de orden 0 del registro, respectivamente. 

• En este caso, el momento temporal determinado por regresion es el de orden 0 y en funcion de 
este se calculan los de orden superior, en funcion de k x y k-> . Se ha pensado esta alternativa 
para controlar de alguna manera el diseno de las funciones de amplitud en cada franja, 
variando el modo en que la energia se distribuye, pero no su intensidad. 

• A fin de asegurarse que la funcion de amplitud sea del tipo gamma los parametros k x , k 2 y 
t 0 deben satisfacer ciertas condiciones que se explicitan en el trabajo. 

• Se especifican los valores medios de los parametros k x , k 2 y t 0 obtenidos para sismos 

cuyanos, obtenidos de la experimentation numerica llevada a cabo con el modelo de ajuste no 
lineal! sus desviaciones estandar denotan una gran variabilidad de los mismos. Los valores 
explicitados en el trabajo son solo indicatives y pueden ser tornados como punto de partida 
para definir los que se ajusten efectivamente al problema a resolver en cada caso. 

• Se ha experimentado con varias combinaciones entre estos parametros, de lo que resultaron 
algunas reglas simples que se pueden sugerir a fin de disenar funciones de caracteristicas 
especificas determinada. Si se desean funciones gamma angostas se debe considerar valores 
bajos de k 2 y si se desean que las funciones gamma sean simetricas, se debe considerar un 
valor de k x cercano al limite superior del intervalo de valores que resulta de aplicar los 
requisitos de funcion gamma. 

Con los conceptos antes expuestos, se ha confeccionado un programa de computadora al 
que se denomino PREDICE, en el que se requieren cuantificar los siguientes parametros que 
definen las caracteristicas de la serial generada : 
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Numero de franjas frecuenciales para registros de la zona cuyana se sugiere considerar tres 
franjas para sismos de banda estrecha y cinco para sismos de banda ancha. En el caso de sismos 
cercanos (menos de 30 km de distancia epicentral) conviene optar por registros de banda angosta 
y para sismos lejanos es apropiado considerar registros de banda ancha. Como tambien las 
condiciones locales del suelo determinan el contenido espectral, se pueden tener en cuenta estas 
particularidades y por ende hacer otras consideraciones. 

Contenido energetico de las franjas frecuenciales la distribution energetica en cada franja 
frecuencial para sismos moderados, cercanos y producidos por deslizamiento de fallas 
superficiales pueden seB 68 %, 28% y 4% para sismos de banda estrecha (division en tres franjas) 
y del 4%, 23%, 51%, 19% y 3% para sismos de banda ancha (division en cinco franjas 
frecuenciales). Estos son los valores promedios de la base de datos usada pero pueden camhiarse 
a voluntad segun la necesidad de la simulation. 

Diseno de la funcion de amplitud en cada franja, que puede realizarse mediante la aplicacion de 
formulas cerradas para la determination de los parametros de la funcion gamma tradicional en 
funcion de los tres momentos espectrales determinados por regresion lineal o bien en funcion de 
los dos parametros propuestos k x y k 2 que simplemente relacionan los momentos espectrales de 
orden 1 y 2 con el momento de orden 0 del registro, respectivamente. 

Election de la funcion de densidad espectral de potencia, en cada franja uniformemente 
modulada- puede ser considerada constante o en funcion del espectro de Boore. 

La aceleracion maxima, duration efectiva y duration total, se determinan mediante formulas de 
regresion con los parametros determinados para la base de datos usada o bien por consideraciones 
especificas del lugar de estudio y fines de la simulation. En cuanto a la duration de los sismos 
generados, los parametros determinados por regresion lineal condujeron a resultados 
satisfactorios en todos los casos analizados! de todos modos, se pueden determinar a priori y por 
otras consideraciones, valores distintos a los determinados por regresion, que como se acoto 
repetidas veces, resumen caracteristicas de la base de datos utilizada en este trabajo. 

• En cada franja se controla que la energia del proceso uniformemente modulado generado 
tenga tenga minimas diferencias con respecto a la energia asignada a la banda frecuencial. 
Una vez obtenido el registro artificial se controla que la energia total del mismo sea similar a 
la obtenida por regresion. En general, se obtienen valores picos del orden de los determinados 
por la formula de atenuacion para la zona cuyana para suelo duro! el programa admite la 
posibilidad de una correction que no altere las caracteristicas dinamicas basicas de la historia 
generada. 

• Una vez generado el registro artificial se procede a su correction de linea base, tal como se 
realiza con los registros reales, mediante el clasico procedimiento de Berg y Housner. 

Se realizo una amplia experimentation numerica, aplicando el modelo principalmente a la 
zona bajo estudio. Dio resultados satisfactorios para predecir posibles movimientos en un 
determinado sitio, tanto moderados como intensos. Debido a la exigua base de datos con que se 
contaba al momento de realizar esta investigation, la mayor confiablidad corresponde a sismos 
moderados en el area de la capital mendocina y del Gran Mendoza, en la Republica Argentina. Se 
considera de todos modos que la metodologia propuesta es valida para cualquier sitio, con solo 
considerar los datos apropiados. 
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RESUMEN 



Con el apoyo financiero de la Direction Academica de Investigation (DAI) de la Pontificia 
Universidad Catolica del Peru (PUCP), se efectuo en el Laboratorio de Estructuras de la PUCP, 
un estudio experimental que tuvo por objeto impedir la falla por deslizamiento en los muros de 
albanileria armada construidos con bloques de concreto vibrado, cuando ellos se ven sujetos a 
terremotos. En primer lugar, se analizaron las posibles causas que pueden originar la falla por 
deslizamiento, para luego plantear soluciones economicas y sencillas de ejecutar en obra. 

Posteriormente, estas soluciones se aplicaron en la construction de cinco muros, los que 
luego fueron ensayados a carga lateral ciclica obteniendose resultados satisfactorios. 



1. INTRODUCCI6N 



En los edificios estructurados con muros de albanileria armada sometidos a terremotos, 
puede ocurrir, por una serie de razones, una falla por deslizamiento (cizalle o corte’friccion). Esta 
forma de falla es muy peligrosa, porque la section transversal del muro se reduce 
significativamente figura 1, debido a la action simultanea de las cargas sismicas coplanares y 
perpendiculares al piano del muro, y tambien, porque esta falla puede causar incluso la rotura del 
refuerzo vertical y el disloque del muro por encima del piano de falla, lo que genera una doble 
excentricidad de la carga vertical. De esta manera, el objetivo del proyecto es tratar de eliminar la 
falla por deslizamiento, o al menos retardar su aparicion hasta que los muros alcancen la 
distorsion angular inelastica maxima permitida para las construcciones de albanileria armada 
por nuestra norma sismica (0.005, Ref.l) cuando ocurren sismos severos, transformando este tipo 
de falla en otras menos peligrosas, como som las fallas por flexion o por fuerza cortante. 



2. POSIBLES CAUSAS DE LA FALLA POR DESLIZAMIENTO 



2.1. F alta de adherencia bloque-mortero. 

Los bloques tienen una suction relativamente elevada, lo que conduce a que el mortero 
colocado sobre el endurezca rapidamente, perdiendose la adherencia con el bloque de la hilada 
inmediata superior, originandose de este modoun piano potential de deslizamiento (figura. 2). 

Recibido: Julio de 2003 
Aprobado 1 Marzo de 2004 
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Cabe indicar que los bloques se asientan en su estado natural (secos) y el unico tratamiento que 
reciben es la limpieza de particulas sueltas. 




Figura 1 Falla por deslizamiento. 




Figura 2 



2.2. Traslape del refuerzo vertical. 



Tradicionalmente, para facilitar el proceso constructive de la albanileria armada, se 
acostumbra dejar espigas ancladas en la cimentacion, para despues de terminar con la 
construction de la albanileria, insertar las barras verticales traslapandolas con las espigas, 
creandose de esta manera un piano de debilidad en la zona donde terminan las espigas, por donde 
puede ocurrir la falla por deslizamiento figura 3. 

Para evitar el problema que ocasionan los traslapes, en otros paises se fabrica bloques H 
figura 4, que permiten el uso de refuerzo vertical continuo, con lo cual se descongestionan las 
celdas reduciendose la posibilidad de que se formen cangrejeras en el concreto liquido (“grout” o 
material con que se rellenan los muros). El emplear refuerzo vertical continuo y bloques 
convencionales (no en forma de H), retarda el proceso de asentado durante la primera jornada de 
trabajo, por el tiempo que demora insertar y bajar al bloque a traves del refuerzo vertical. 




Figura 3 



Estudio experimental de cuatro alternativas para impedir la falla por 
deslizamiento en los muros de albanileria armada 
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2.3. Falta de adherencia grout-cimentacion. 



Puede ser causado por una serie de razones a) tratamiento inadecuado de la superficie del 
cimiento (rayado poco profundo, presencia de lechada superior, falta de limpieza al iniciar la 
construction de la albanileria y tambien antes de vaciar el grout); b) segregation del grout figura5 
producida por: su alto contenido de agua, compactacion inadecuada, altura excesiva de vaciado; 
etc. 




Figura 5. Segregation del grout. 



2.4. Contraction de secado del grout. 



Lo ideal es que quede integrado como una sola unidad la union bloque'grout-mortero - 
acero. De separase el grout del bloque y del acero, el refuerzo perderia efectividad ante las 
acciones de flexion, cizalle, etc. El polvo de la arena y la cantidad de agua que requiere el grout, 
incrementan este fenomeno. 



2.5. Distribution uniforme del refuerzo vertical. 



Tradicionalmente, se distribuye el refuerzo vertical a lo largo del muro. Puesto que el 
refuerzo en los extremos del muro es pequeno, la base del muro rota con facilidad. Este hecho 
provoca el aplastamiento de los talones y el corrimiento de las fisuras horizontales por flexion, las 
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que al encontrarse por la action ticlica de la carga sismica, terminan transformandose en una 
falla por deslizamiento figura 6. 




Figura 6 Transformation de la falla por flexion en falla por deslizamiento. 



2.6. Uso de cintas de mortero longitudinales y verticales. 



El grout no rellena el espacio entre las cintas, lo que debilita la union bloquennortero, con 
la consecuente formation de fisuras. 




Figura 7 Cintas de mortero en las juntas. 



3. SOLUCIONES PARA EVITAR LA FALLA POR DESLIZAMIENTO 



En este proyecto se trata de encontrar soluciones al problema del deslizamiento, que sean 
economicas y faciles de ejecutar en obra. Sin embargo, en este acapite se describen todas las 
posible soluciones, ordenadas en funcion del costo y de la dificultad constructiva, a fin de tener un 
panorama amplio que permita eliminar las causas del problema, indicadas en el acapite 2, y 
seleccionar las soluciones mas adecuadas. 

1. Profundizar el rayado de la superficie superior de la cimentacion (minimo 6 mm), limpiarla de 
particulas sueltas y de lechada, para luego humedecerla antes de asentar los bloques y antes 
de vaciar el grout. El objetivo de esta operation es mejorar la adherencia y la resistencia a 
corte'friccion en las interfases mortero'timentacion y grout-cimentacion. El muro adoptado 
como patron de comparacion (Ml, traditional) tuvo un rayado conventional. 

2. Manteniendo la cuantia total de refuerzo vertical, concentrar la mayor area de acero en los 
extremos del muro, empleando en la zona central una cuantia minima igual a 0.001. El 
objetivo es disminuir el giro por flexion para controlar el corrimiento de las fisuras 
horizontales. El muro patron (Ml) tuvo una distribution uniforme del refuerzo vertical. 
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3. Invertir la cara de asentado de los bloques correspondientes a la primera hilada, creando 
aberturas en las 3 tapas transversales, de manera que el grout al pasar a traves de esas 
aberturas, rellene una mayor area en su zona de contacto con la cimentacion. 




4. Compactar al grout empleando una madera con section transversal de l”x 2” en reemplazo de 
la varilla lisa de 5/8” que se utiliza tradicionalmente. El objetivo es tratar que el grout se 
expanda lateralmente en mayor grado al apisonarlo con la madera. Esta solution es aplicable 
solo cuando se utiliza bloques de 19 cm. de espesor, con celdas grandes. 

5. Disminuir la cantidad de confitillo en el grout, de una proportion cemento- arena- confitillo 1- 
2 Vz- 1 Vi (traditional) a 1 : 3^ 1, y tamizar la arena a traves de la malla #200 para eliminar el 
polvo. El objetivo es atenuar el problema de segregation y contraction de secado. 

6. Similar a la solution 5, pero empleando cemento puzolanico y reduciendo el slump de 10” 
(traditional) a 9”. Se tiene conocimiento que el cemento puzolanico atenua la contraction de 
secado, mientras que al disminuir la cantidad de agua (al bajar el slump), se reducen los 
problemas de segregation y contraction de secado. 

7. Similar a la solution 5, pero agregando 1/10 de cal para retardar la contraction de secado. 

8. Tratamiento con agua para incrementar la adherencia bloqueunortero y atenuar la 
contraction de secado del grout. Con una brocha, humedecer las juntas 2 horas despues de 
haber asentado los bloques, una vez al dia, hasta el dia del vaciado del grout. Antes del 
vaciado, regar internamente las celdas, encofrar las ventanas de limpieza, vaciar y compactar 
al grout, luego curar al muro regandolo externamente una vez al dia durante una semana. 

9. Adicionar espigas en las celdas que carecen de refuerzo vertical, ancladas en la cimentacion. 
El objetivo es incrementar la resistencia a corte'friccion en las dos primeras hiladas del muro 
y en la interfase muro'timentacion, ademas, que las fisuras por flexion se originen en las 
partes mas altas del muro, de tal manera que la cimentacion no estorbe su trayectoria. 

10. Dejar orificios en la cimentacion, en las zonas ocupadas por las celdas. El objetivo es crear 
llaves de corte con el grout, para transferir la fuerza cortante hacia la cimentacion. Esta 
solution es menos efectiva que la numero 9, debido a que la falla de deslizamiento se puede 
trasladar a las hiladas superiores. Tambien, es dificil hacer los orificios y limpiarlos. 

11. Emplear refuerzo vertical continuo, con bloques recortados en forma de H para facilitar el 
asentado. El objetivo es tener mas espacio en las celdas para poder compactar al grout, 
evitando la posibilidad que se formen cangrejeras, asi como evitar los problemas que causan 
la discontinuidad del refuerzo vertical. Esta solution mejora al 100% las juntas verticales al 
quedar llenas de grout, pero debe anadirse otras soluciones para mejorar la union horizontal 
bloqueunortero, grout-cimentacion y atenuar la contraction y segregation del grout. 

12. Usar grout fino (cementcarena'cal V- 3- Vi) en la mitad inferior del muro y grout grueso 
(cemento'arena'confitillo'cal 1- 3- 1- 0.1) en la mitad superior. Con esta solution se evitaria el 
problema de segregation del grout en la base del muro, pero se incrementaria el problema de 
contraction de secado en la mitad inferior. Esta solution es de dificil aplicacion en obra, por 
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tener que variarse la mezcla en un mismo piso, sin embargo, no se descarta ya que el grout se 
vacia siempre en dos etapas, con un lapso de Vi hora entre ellas. 

13. Uso de aditivo en polvo Wacker (culminal mas vinnapas) en el mortero, en reemplazo de la 
cal. Este aditivo mejora notoriamente la adherencia bloque - mortero (37%, segun la referenda 
4), pero es costoso (segun la referenda 4, el costo de la albanileria se eleva en 18%). 

14. Uso de aditivo Intraplast (de Sika) en el grout para mejorar su plasticidad y disminuir su 
contraccion de fragua. En el mercado existen aditivos que permiten reducir la cantidad de 
agua en el grout, o superplastificantes (Rheobuild 1000 de MBT) para producir grout 
autocompactado (no necesitan vibradora ni chuzo). Pero, estas soluciones son costosas. 

15. Uso de vibradora de aguja de 14” para compactar al grout. Para esto se inserta apagada la 
vibradora a traves de las celdas de los bloques, se vacia el grout, se enciende la vibradora y se 
la retira lentamente, debiendose hacer esta operation celda por celda. Esta solution se 
descarta por ser poco practica. 

16. Low Lift Grouting. Consiste en construir la primera mitad del muro, llenar con grout al dia 
siguiente, al tercer dia construir la segunda mitad y al cuarto dia volver a llenar con grout. 
Esta solution atenuaria la posibilidad de segregation del grout, ya que el vaciado se hace 
desde una altura de 1.2 m, pero no se utiliza por la posibilidad de que se generen juntas frias 
entre las capas de grout y tambien porque el proceso constructive se retarda un dia. Cabe 
indicar que para limpiar los desperdicios de mortero que caen a traves de las celdas durante 
la construction de la segunda mitad del muro, habria que abrir ventanas en los bloques 
ubicados a media altura del muro, esto no es posible porque este tipo de muro es caravista. 



4. PROBETAS PRELIMINARES DE GROUT 



El objetivo que se busca en este acapite es definir el grout a utilizar en los muros 
mejorados, tratando de reducir los problemas de contraccion de secado y de segregation. Para esto 
se fabricaron 10 mezclas distintas (Gi en la figura 8), empleando como molde un pilar de bloques 
U (medios bloques), con una altura de 2.4 m (figura 10). Estos moldes permitieron la 
transferencia natural de agua desde el grout hacia los bloques. Por otro lado, se utilizo encofrados 
transparentes con el objeto de observar lo que ocurria con el grout durante su vaciado y 
compactacion. Al retirar el encofrado, pudo observarse si existian los problemas mencionados. 





Figura 8 Moldes con bloques U. 



Las mezclas (Gi) utilizadas fueron. 



Gl. Grout tradicional, cemento normal-arena normakconfitillo 1 1- 2 %'■ 1 !4, slump 10”. 
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• G2. Cemento normal-arena tamizada-confitilkv P 3: 1, slump 10”. 

• G3. Similar a G2, pero con 0.1 de cal aerea normalizada. 

• G4. Similar a G2, pero con aditivo Intraplast. 

• G5. Similar a G2, pero con cemento puzolanico IP. 

• G6. Similar a G5, pero con 9” de slump. 

• G7. Similar a G6, pero con 0.1 de cal. 

• G8. Similar a G6, pero sin confitillo (cemento-arena P 3) en la mitad inferior y con confitillo 
(cemento-arena-confitillo P 3^ l) en la mitad superior. Slump 9”, sin cal. 

• G9. Similar a G8, pero con !4 de cal en la mitad inferior y 0.1 de cal en la mitad superior. 

• G10. Similar a G5, pero con tratamiento con agua (regado interno de celdas antes del 

vaciado). 

Externamente se pudo observar que el grout tradicional (Gl) tuvo los mayores problemas 
de segregacion y contraccion de secado, mientras que estos problemas se redujeron en las mezclas 
G5, G6, G7 y G10 (con cemento puzolanico y arena tamizada), pero, las mezclas G6 y G7 (con 9” 
de slump), se notaron muy secas, lo que podria provocar cangrejeras internas en los muros, por lo 
que es conveniente utilizar 10” de slump. Para revisar internamente al grout, el pilar fue izado y 
volcado con un puente grua, para luego cortar su base con una amoladora a la altura de la 
primera hilada. De la base del pilar se extrajeron muestras de grout, con el objetivo de observar 
su estructura interna, cabe destacar que la mayor concentracion de piedras (segregacion en la 
base) se encontro en la probeta Gl (figura 10). Como conclusion de estos ensayos preliminares, se 
utilizo en los muros mejorados: G2 para M2, G5 para M3, y G10 para M4 y M5. 




Figura 9 Preparacion del grout y construction del pilar. 
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Figura 10 Extraction de muestra de grout de la base del pilar y muestras G5 y Gl. 



5. CONSTRUCTION DE LOS MUROS 



En total se construyeron cinco muros a escala natural (ver las figuras 11 a 15), 

empleandose la tecnica constructiva tradicional para el muro Ml y manteniendose constante los 

siguientes parametros (la mayor parte de ellos se emplean tradicionalmente en obras reales). 

• Geometria del muro las dimensiones fueron 2.4 x 2.4 x 0.14 m. 

• Aparejo se utilizo un aparejo de soga (14 cm. de espesor). 

• Calidad del bloque se utilizo bloque tipo Bill, serie 15.14, procedente de la fabrica FIRTH, con 
ranuras en las tapas transversales que permitieron el paso del refuerzo horizontal. 

• Tratamiento de los bloques las superficies de asentado fueron limpiadas con una brocha de las 
particulas sueltas antes de asentarlos en su estado natural (secos). 

• Calidad del mortero cemento-caharena 1- X A'- 4. Se uso arena gruesa tamizada entre las mallas 
ASTM #4 y #200. El cemento fue Portland tipo 1, no se uso cemento puzolanico porque 
proporcionaria un color marron a las juntas. 

• Calidad de la cal se utilizo cal aerea hidratada y normalizada de Agregados Calcareos. 

• Grosor de las juntas las juntas horizontales y verticales (brunadas) tuvieron 1 cm. de grosor. 

• Calidad del acero de refuerzo se utilizo acero corrugado de grado 60 de SIDERPERU. 

• Batido del grout se utilizo una mezcladora. El slump medido en el cono de Abrams fue 10”. 

• Construccion de la albanileria se construyo en 2 jornadas con una mano de obra calificada. 

• Tecnica de limpieza de las celdas los bloques de la primera hilada tuvieron ventanas de 
limpieza de 3”x4” en cada celda. Terminada cada jornada de trabajo, se limpio internamente 
las celdas con una varilla, sin tocar al muro. Antes del vaciado del grout, la base se rasqueteo 
con un tortol, las celdas se limpiaron con aire comprimido y se humedecio la cimentacion. 

• Tecnica de vaciado del grout fue vaciado en dos capas, cada una fue compactada con una 

varilla lisa de 5/8”, transcurrido 5 minutos, cada capa fue recompactada. El intervalo de 

vaciado entre capas fue 14 hora. El borde superior de cada capa quedo aproximadamente a 10 
cm. por debajo de la superficie superior correspondiente a la ultima hilada de cada jornada. 
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• Cuantia de refuerzo. El refuerzo horizontal tuvo una cuantia de 0.0025 (l f 3/8” @ 20 cm.) y la 
cuantia de refuerzo vertical total fue 0.0024. 

• Vigas. Las vigas de cimentacion y las vigas soleras tuvieron la misma seccion transversal, la 
misma calidad de concrete y el mismo refuerzo. 

Matriz de Variantes en la tabla 1, se indica las variantes que tuvieron los 5 muros, 
indicandose que el muro Ml, construido con la tecnica tradicional, corresponde al patron de 
comparacion, asimismo, las mejoras varian progresivamente, desde M2 hasta M5. 



Tabla 1. Matriz de Variantes. 



VARIANTE 


Ml 


M2 


M3 


M4 


M5 


Rayado de la 
cimentacion 


Normal 


Profundo 


Profundo 


Profundo 


Profundo 


Distribution del 
refuerzo 

vertical 


Distribuido 


Concentrado 


Concentrado 


Concentrado 


Concentrado 


Tamizado de la arena 
del grout 

por mallas #4 y #200 


no 


si 


si 


si 


si 


Tipo de junta vertical 


Cintas 


Llena 


Llena 


Llena 


Cintas 


Traslape del refuerzo 
vertical 


100% en la 

misma 

seccion 


alterno 
h = 70 cm. 
y 95 cm. 


alterno 
h = 70 cm. 
y 95 cm. 


alterno 
h = 70 cm. 
y 95 cm. 


continuo, 

sin 

traslapes 


Grout 

cemento-arena - 

confitillo 


l: 2 y 2 : 1 14 


l: 3: 1 


l: 3: 1 


l: 3: 1 


l: 3: 1 


Tipo de cemento en el 
grout 


Portland 


Portland 


puzolanico 


puzolanico 


puzolanico 


Tratamiento con 
agua en las 
juntas, celdas 
internas y muro 


no 


no 


no 


si 


si 


Espigas adicionales 


no 


no 


no 


si 


si 


Bloque de hilada 1 


normal 


normal 


invertido 


invertido 


H 





Figura 11. Caracteristicas comunes de los 5 muros. 
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Figura 12. Refuerzo vertical en Ml y M2. 
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Figural3 Refuerzo vertical, bloque de la primera hilada y espigas adicionales en M3, M4 y M5. 



Estudio experimental de cuatro alternativas para impedir la falla por 
deslizamiento en los muros de albanileria armada 



151 




Figural4. Construction de muros Ml (izquierda), M3 (centro) y M5 (derecha). 




Figura 15. Tratamiento con agua en M4 y M5. Ver el parrafo 3.8. 



6. ENSAYOS DE CONTROL 



6.1. Grout 



Empleando como moldes a los bloques forrados con papel filtro figura 16, para permitir la 
transferencia de agua del grout hacia los bloques, se fabricaron 9 probetas de 10x10x19 cm., tres 
provinieron del muro Ml (tradicional), las otras tres de M2, mientras que las tres restantes 
fueron hechas con el grout de M3, M4 y M5. Estas probetas no fueron curadas, y luego de 28 dias 
fueron ensayadas a compresion, obteniendose los siguientes resultados promedios. 
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• Grout l: 2 V* '■ 1 X A (cemento Portland tipo I- muro Ml): f' c = 240 kg/cm 2 

• Grout l: 3 :l (cemento Portland tipo I - muro M2): f'c = 193 kg/cm 2 

• Grout l: 3: 1 (cemento puzolanico tipo IP- muros M3, M4 y M5) : f'c = 197 Kg./cm. 2 

Para los tres casos, la resistencia supero al valor minimo especificado por la Ref.2 (140 Kg./cm. 2 

). 

6.2. Prismas de albanileria simple (no reforzada) 



Se construyeron un total de 6 pilas y 6 muretes, los cuales fueron ensayados a compresion 
axial y a compresion diagonal, respectivamente. Las pilas (14x39x60 cm.) estuvieron compuestas 
por 3 hiladas y los muretes por 4 hiladas (14x80x80 cm.). Tres pilas y tres muretes fueron 
rellenados con el grout tradicional de Ml (l: 2 14: 1 %) i mientras que el resto de prismas fueron 
rellenados con grout similar al empleado en M4 y M5 (l: 3: 1, con cemento puzolanico), ademas de 
recibir el tratamiento humedo especificado en 3.8 para M4 y M5. En la tabla 2 se indica los 
resultados, mientras que en la figura 17 se aprecia que los muretes rellenos con grout tradicional, 
tuvieron una falla escalonada a traves de las juntas, en tanto que los rellenados con grout 1:3:1, 
tuvieron algunas grietas que atravesaron los bloques, lo que implica una mejora en la adherencia 
mortero-bloque, producto del tratamiento que tuvieron esos prismas, sin embargo, el grout de esos 
prismas se noto que aun estaba humedo, pese a que tenia 28 dias de edad, lo que debe haber 
repercutido en los resultados, principalmente en la resistencia a compresion de las pilas. 



Tabla 2. Resistencia de los Prismas de Albanileria Simple 



Prisma 


Resistencia 


Grout V 2 % '■ VA Cemento I 


Grout l: 3: 1 Cemento IP 


Pilas 


Promedio fm 


104 kg/cm 2 


84 kg/cm 2 


Dispersion 


3.0 % 


6.0 % 


Caracteristica f'm 


101 kg/cm 2 


79 kg/cm 2 


Muretes 


Promedio vm 


13.1 kg/cm 2 


12.5 kg/cm 2 


Dispersion 


8.6 % 


3.5 % 


Caracteristica v'm 


12.0 kg/cm 2 


12.1 kg/cm 2 




Figura 17. Falla escalonada en muretes rellenos con gout l: 2 54 : 1 14 (cemento tipo I, izquierda) y falla 
mixta en muretes rellenos con grout l: 3 : 1 (cemento IP, derecha). 
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7. TECNICA DE ENSAYO DE CARGA LATERAL CICLICA 



Los cinco muros fueron ensayados a carga lateral ciclica figura 18, controlando el 
desplazamiento horizontal de la viga solera en 10 fases de distinta amplitud tabla 3 y figura 19. 
La carga lateral (V) se aplico mediante un actuador dinamico, mientras que los desplazamientos 
fueron medidos con LVDT (D en la figura 18). Las fisuras de cada fase se pintaron de distinto 
color. 



Tabla 3. Fases del Ensayo de Carga Lateral Ciclica. 



Fase 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


D1 

(mm) 


0.25 


1.00 


2.50 


5.00 


7.50 


10.00 


12.50 


15.00 


17.50 


20.00 


No. 

Ciclos 


1 


2 


3 


3 


3 


3 


4 


4 


4 


4 


Color 


- 


rosa 


lila 


negro 


Rojo 


azul 


marron 


verde 


naranja 


celeste 




Figura 18. Vista global de los dispositivos utilizados en el ensayo de carga lateral ciclica. 







154 



Angel San Bartolome, Wilson Silva, Eliana Melendez y Gino Castro 




Figura 19 Lazos histereticos V-Dl de los muros Ml y M4 (arriba) y envolvente V-Dl del 
primer ciclo para los 5 muros (derecha). 



8. COMPORTAMIENTO DE LOS MUROS 



A continuation se describe el comportamiento que tuvieron los muros durante el ensayo de 
carga lateral ciclica (figura 20), poniendo enfasis en el problema de deslizamiento. Cabe indicar 
que la norma sismica (referenda l) permite una distorsion inelastica maxima de 0.005 para los 
entrepisos de los edificios de albanileria armada, que en este caso se produce en la fase 7 del 
ensayo (distorsion = Dl/h = 12.5 mm/ 2500 mm = 0.005, figura 19‘derecha), por lo que si la falla 
por deslizamiento se produjera en una fase mayor, careceria de importancia! sin embargo, debido 
a que en el experimento no se contemplo una serie de factores reales, como la action dinamica con 
una gran cantidad de ciclos que imprime un sismo real, se considero prudente tomar como limite 
maximo un desplazamiento D1 = 15 mm (fase 8), asociado a una distorsion de 0.006. 
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Figura 20 Fase 10 D1 = 20 mm 



En el muro Ml se presentaron 2 pianos de deslizamiento, el primero estuvo ubicado en la 
base y se produjo en la fase 5, mientras que el segundo se ubico en la base de la quinta hilada, 
donde terminaban las espigas, esta falla se presento en la fase 7. En M2, el deslizamiento se 
produjo a traves de la base durante la fase 5. En M3, la falla por deslizamiento en la base se initio 
en la fase 7, en tanto que M4 no tuvo problemas de deslizamiento. En M5, la falla por 
deslizamiento se produjo durante la fase 9 en la base de la tercera hilada, donde terminaban las 
espigas. 



8.1. Observaciones Post-Ensayo 



Terminado el experimento, se procedio a inspeccionar internamente el estado de cada 
muro (figuras 21 y 22). Esta operation se hizo con el objetivo de detectar cangrejeras, segregation 
del grout, el grado de adherencia mortero-bloque-grout-acero, el estado del refuerzo, etc. 
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Figura 21 Muro Ml. Falta de adherencia bloque-grout y juntas verticales vacias. 




Figura 22 Pandeo del refuerzo vertical por deslizamiento en la base de M3 y elimination del 
pandeo al evitarse el deslizamiento en M4. 



9. CONCLUSIONES 

El objetivo principal de este proyecto, que fue controlar o eliminar la falla por 
deslizamiento transformandola en otra menos peligrosa como es la falla por flexion o por corte, 
pudo lograrse. Para ello es imprescindible el uso de espigas de 3/8”, colocandolas en aquellas 
celdas donde no exista el refuerzo vertical calculado por flexion, y que estas espigas tengan 
distintas longitudes por encima de la cimentacion (se propone 30 y 50 cm.), intercalandolas cada 
20 cm., que es el espaciamiento entre las celdas del bloque, y anclandolas 30 cm. en la 
cimentacion. 

Cabe destacar que ninguna mejora que se haga en este tipo de muro sera efectiva o se 
traducira en incrementos de resistencia, si es que no se evita la falla por deslizamiento! por 
ejemplo, de nada servira incrementar el refuerzo horizontal, o la adherencia bloque'grout-acero, o 
la elimination de cangrejeras internas (como en M5), si no se elimina la falla por deslizamiento. 



9.1. Proceso Constructive 

1. El rayado profundo de la superficie de la cimentacion (M2 a M5), mejoro la union mortero - 
cimentacion. Sin embargo, en la base del muro las fisuras horizontales se trasladaron a la 
union bloqueunortero, por la baja adherencia que presentaron esos materiales en todas las 
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hiladas. De requerirse mayor resistencia sismica, el problema de adherencia bloqueunortero 
puede solucionarse (referenda 4) agregando al mortero el aditivo en polvo Wacker. 

2. El rayado profundo de la superficie de la cimentacion y el empleo en la primera hilada de 
bloques recortados en sus tapas transversales para que circule el grout ocupando una mayor 
area, no incremento la resistencia a corte'friccion en la base del muro, debido a que la base 
del grout fue cubierta por una capa lisa de lechada de cemento proveniente del mismo grout. 




3. La disminucion del confitillo en el grout, manteniendo el slump en 10”, de una proportion 
tradicional cemento-arena-confitillo V- 214: 114 (usada en Ml), a una dosificacion 1- 3- 1, 
disminuyo el problema de segregation del confitillo en la base del grout. 




4. El tratamiento con agua hecho a los muros M4 y M5 (curado de juntas con brocha, asi como el 
regado interno antes de vaciar el grout y regado externo despues del vaciado) mejoro 
ligeramente la adherencia bloque-grout, al atenuar la contraccion de secado del grout. 

5. El empleo de cemento puzolanico IP en el grout del muro M3, redujo ligeramente el problema 
de contraccion de secado del grout, observandose una mejor adherencia bloque'grout que la 
existente en los muros Ml y M2, donde se utilizo cemento Portland tipo I. Esta mejora no se 
reflejo en incrementos de resistencia porque M3 fallo por deslizamiento. 

6. En este proyecto no se presentaron cangrejeras en el grout en la zona de traslapes de las 
barras verticales, debido a que esas barras tuvieron poco diametro 04”), por lo que no 
causaron mayor congestion en las celdas de los bloques. 

7. Todas las barras horizontales quedaron desplazadas hacia uno de los lados del muro, debido a 
que se instalaron sueltas y se movieron durante el vaciado y compactacion del grout. Esto 
podria causar con el tiempo problemas de corrosion y falta de adherencia con el grout. Para 
solucionar este problema se sugiere utilizar barras verticales continuas, con lo cual, las barras 
horizontales pueden amarrarse con la verticales, esta solution exigiria el empleo de bloques H 
para facilitar el asentado, o que las fabricas produzcan bloques con menor separation entre 
las ranuras existentes en sus tapas transversales (lugar donde se aloja el refuerzo horizontal). 

8. El espacio existente entre las cintas verticales de mortero, empleadas en el muro Ml 
(tradicional), no fue rellenado por el grout, lo que dio lugar a una zona de debilidad. En el 
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resto de muros, se utilizo junta vertical llena, pero esta mejora no se reflejo en la resistencia 
de los muros (excepto en M4), debido a que fallaron por deslizamiento. 

9. El empleo de arena seca colocada en la base del muro durante su construction, para evitar 
que los desperdicios de mortero queden pegados sobre la cimentacion formando juntas frias en 
la interfase grout'timentacion, deberia descartarse, debido a que la arena seca absorbe el 
agua del mortero debilitando en la base del muro las uniones bloquennortero y mortero - 
cimentacion. 

En reemplazo de la arena, es conveniente utilizar retazos de plasticos. 



9.2. Espigas Verticales 



En este proyecto se utilizaron dos tipos de espigas verticales ancladas en la cimentacion y 
son l) Las necesarias para traslapar las barras verticales requeridas por flexion (tipo A); y, 2) las 
adicionales para incrementar la resistencia a corte 'friction en la base del muro (tipo B, de 3/8” de 
diametro y 30 cm. de longitud), obteniendose las siguientes conclusiones. 

1. El empleo de espigas tradicionales tipo A de la misma longitud (barras traslapadas en la 
misma section transversal), causo en el muro Ml dos pianos de deslizamiento- en la base y 
donde terminan las espigas, lo que resulta muy peligroso, especialmente cuando el muro este 
sujeto a cargas sismicas perpendiculares a su piano, en vista que se comportaria como una 
barra triarticulada (inestable). Incluso, los bloques ubicados en los extremos de Ml llegaron a 
cizallarse horizontalmente en su parte intermedia, en la zona donde terminaban las espigas. 

2. El empleo de espigas tipo A con distintas longitudes (barras traslapadas en distintas 
secciones transversales) ubicadas alternadamente, evito en los muros M2 y M3 el segundo 
piano de deslizamiento! en estos muros la falla por deslizamiento se concentro en la base. 

Cabe destacar que en M3 la falla por deslizamiento se initio a partir de la fase 7 del 
ensayo (D1 = 12.5 mm), coincidente con la distorsion maxima permitida por la norma E'030 para 
sismos severos, mientras que en Ml el deslizamiento se initio en una etapa temprana del ensayo 
(fase 5, D1 = 7.5 mm), es decir, en muros reales como el Ml, la falla por deslizamiento podria 
presentarse ante la action de sismos moderados. 

3. En el muro M4, donde se utilizo espigas tipo A con distintas longitudes y se adiciono espigas 
tipo B en aquellas celdas que carecian de refuerzo vertical, el comportamiento sismico mejoro 
notablemente. En este muro no se produjo la falla por deslizamiento, sino mas bien una falla 
mixta por flexion y corte, elevandose la resistencia maxima en 21% con relation a la 
alcanzada por Ml. Cabe destacar que en M4, no se produjo el pandeo de las barras verticales 
extremas, a diferencia del resto de muros que tuvieron falla por deslizamiento. 

4. En el muro M5, donde se utilizo barras verticales continuas y espigas tipo B en aquellas 
celdas que carecian de refuerzo vertical, se evito la falla por deslizamiento en la base, pero 
esta se traslado a la base de la tercera hilada, donde terminaban las espigas, por lo que es 
necesario que estas espigas tengan distintas longitudes. Cabe destacar que a diferencia de 
M4, el muro M5 carecio de las espigas tipo A, por lo que en la zona de traslape de las barras 
verticales, M5 tuvo la mitad de refuerzo vertical que es el elemento que controla la falla por 
deslizamiento. 

Adicionalmente, el deslizamiento se initio despues que los talones se trituraron durante la 
fase 9, y la trituration se debio a que el grout aun estaba humedo producto del tratamiento con 
agua que tuvo M5, pese a que habian transcurrido 28 dias hasta el instante del ensayo. 
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RESUMEN 



Se presents un modelo matematico para identificar las matrices dinamicas de rigidez y 
flexibilidad de estructuras reticulares a partir de datos de mediciones de vibration forzada. En el 
planteamiento se considers que el comportamiento del sistema estructural es lineal, elastico y no 
amortiguado. 

La definition de las matrices dinamicas de rigidez y flexibilidad se consigue al identificar 
mediante mediciones de vibration; la matriz modal, las frecuencias naturales de vibration y la 
frecuencia de excitation. Es necesario proponer la matriz de masas del sistema por no haber sido 
posible deducirla experimentalmente en esta etapa de la investigation. 

Se muestra la aplicacion de la metodologia desarrollada a un modelo estructural sencillo y 
se comentan los resultados obtenidos. 

Palabras Clave 1 identification, mediciones, vibration, rigidez, flexibilidad 



ABSTRACT 



A mathematical model to identify the dynamic matrix of stiffness and flexibility for 
reticular structures from measurement data of forced vibration is presented. The model considers 
that the system is linear, elastic and undamped. 

The dynamic stiffness and the flexibility matrix are defined once the following data is 
obtained by means of vibration measurements; the modal matrix, the natural vibration 
frequencies, and the frequency of the driving force. It is necessary to propose the structure’s mass 
matrix since it is not possible to deduce it experimentally in this stage of the research. Some 
results for a simple structural model are included and discussed. 

Keywords 1 identification, measurements, vibration, stiffness and flexibility 



1. INTRODUCCION 



En general los sistemas estructurales estan sometidos a cargas que varian en el tiempo 
tal como las que se originan por viento, oleaje, sismos y maquinaria. El efecto de estas cargas 
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sobre la estructura origina una respuesta dinamica que es de suma importancia para el ingeniero 
disenador debido a que la respuesta dinamica puede exceder a la respuesta estatica en un factor 
significante. 

El gran avance de la computation como ayuda de calculo ha permitido el desarrollo de 
numerosos modelos matematicos para analizar el comportamiento dinamico de las estructuras 
reticulares. Sin embargo, la complejidad de estos sistemas hace necesario confirmar y mejorar 
estos modelos mediante tecnicas experimentales. La utilization de metodos para identification de 
sistemas permite que los resultados experimentales puedan usarse para definir con mejor 
precision sus caracteristicas dinamicas y calibrar sus modelos matematicos. 

Debido a la importancia de la infraestructura industrial para la economia de Mexico se ha 
hecho necesario el desarrollar algunas metodologias para la determination de las caracteristicas 
dinamicas de estructuras sometidas a carga de maquinaria. Entre estas metodologias destacan 
aquellas que se basan en el analisis de la serial de vibration en el dominio de la frecuencia. A 
partir de estas metodologias es posible evaluar el comportamiento dinamico de las estructuras 
industriales para contar con una estimation realista de los parameters dinamicos que 
caracterizan al sistema estructural y con base en ello determinar las fuerzas dinamicas 
generalizadas que actuan sobre la estructura en condiciones de operation. 

En este trabajo se desarrolla una metodologia para utilizar los resultados de mediciones 
de vibration efectuadas directamente en campo y a partir de estos datos hacer predicciones del 
comportamiento de estructuras sometidas a vibration forzada. 

2. IDENTIFICACION DEL SISTEMA ESTRUCTURAL 



2.1 Analisis de senales 

Mediante un analisis en el dominio de la frecuencia es posible determinar el sentido y la 
interdependencia estadistica de los movimientos de una estructura (Goldman, 1999). Para tal fin 
es necesario estudiar estos movimientos en varios puntos en las frecuencias de resonancia del 
sistema (Espinoza Barreras, 1999). El procedimiento experimental consiste en analizar dos 
senales simultaneamente, en donde uno de los acelerometros se coloca como referenda mientras 
que el otro se situa en puntos estrategicos sobre la estructura con el proposito de definir mejor el 
modo que se pretende medir. 

Como primer paso se evalua la Transformada Finita de Fourier de la serial (empleando el 
algoritmo conocido como Transformada Rapida de Fourier) mediante la siguiente expresion 
(Brighan, 1997): 



j27t " dt 



(i) 



En esta expresion x (f) es la serial registrada, t define la ventana de medicion, j es (*l) V4 y f 
es la frecuencia de analisis. La ecuacion (l) define la forma en que una serial se distribuye con la 
frecuencia. Si al resultado de esta expresion se le calcula su modulo y se eleva al cuadrado se 
obtiene el espectro de potencia de x (t). 

Por otra parte, durante el proceso de medicion se pueden registrar varias senales 
simultaneamente, por lo que se puede calcular su funcion de correlation cruzada, esto es: 

R xy ( T ) = y £ x (t)y(t + ^)dt 



( 2 ) 
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En la Relation 2, y (t+x) representa otra de las senales registradas siendo x un tiempo de 
retraso. Esta funcion permite estimar en que orden una serial se correlaciona con otra, en que 
medida las senales fueron originadas por la misma fuente y con que retraso, y detectar la 
presencia de ruido o senales extranas. Cuando la Ecuacion 2 se aplica sobre la misma serial se 
obtiene la funcion de autocorrelation^ 

R xxW= Ylo X ( t ) X ( t + X ) dt (3) 

Utilizando el concepto de Transformada Finita de Fourier, ecuacion (l), en las Ecuaciones 
2 y 3, se obtienem 



F xy ( f ) = l> xy (t)e- j2ltft dt (4) 

Fxx (f - ) - £ R xx (t)e _J2lcft dt (5) 



Las Expresiones 4 y 5 son, respectivamente, las transformadas de las funciones de 
correlation y autocorrelation. Una aplicacion de la funcion F xy (f) consiste en el calculo del angulo 
de fase entre las componentes imaginaria y real de una serial. Este angulo se define de la 
siguiente manera- 



e xy (f) = tair' 



f Qxy( f ) l 

lCx y (f)j 



( 6 ) 



En donde Q xy (f) y C xy (f) representan respectivamente las partes imaginaria y real de F xy (f). 
La interpretation fisica de este concepto se base en el hecho de que para cada frecuencia natural, 
la funcion de transferencia de fase correspondiente (relation de fase entre ambas senales) es 
cercana a 0 o 180 grados. 

Las transformadas de las funciones de correlation y autocorrelation resultan de utilidad 
en el calculo de las funciones de coherencia y transferencia 1 




Mff 

Fxx(f)F yy (f) 



(7) 



H xy (f ) 



F w (f) 

Fxx (F) 



( 8 ) 



En estas expresiones yxy(f) es la funcion de coherencia, H xy (f) es la funcion de transferencia 
y Fyy(f) denota la transformada de la funcion de autocorrelation de la serial y(t). La principal 
utilidad de la funcion de coherencia radica en que proporciona una medida del grado de linealidad 
entre dos senales. Los valores de esta funcion van de 0 a 1. Cuando se alcanza el valor de 1 se dice 
que existe una perfecta relation lineal entre ambas senales. 
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2.2 Determination de parametros estructurales 



La principal utilidad de los conceptos descritos, dentro de un contexto de analisis 
dinamico estructural, radica en la estimation de las frecuencias y modos naturales de vibration 
que corresponden a los sistemas analizados. Para ello se precede de la siguiente manera- 

a) De los espectros de potencia, ( | x(f) | ) 2 , se determinan las frecuencias en las cuales se generan 
los picos mas grandes. Un pico esta asociado a un valor maximo de ( | x(f) | ) 2 . Estos picos 
representan estimaciones preliminares de las frecuencias naturales de vibration. 

b) Se calculan las Transformadas Finitas de Fourier de las funciones de correlation entre dos 
senales registradas para una misma medicion. A estas transformadas se les obtiene su modulo y 
se eleva al cuadrado. Los resultados se dibujan en funcion de la frecuencia. En estas figuras se 
detectan los valores de frecuencia en que ocurren los picos o valores maximos de | Fxy(f) | 2 . 
Aquellos que coincidan con los detectados en el inciso anterior se podran considerar como 
asociados a las frecuencias naturales de vibration. 

c) Se obtiene la funcion de Transferencia de Fase entre las senales utilizadas en el inciso (b) y se 
dibujan en funcion de la frecuencia. Para aquellos valores detectados como frecuencias naturales 
de vibration se obtiene la magnitud de la Funcion de Transferencia de Fase. De acuerdo con ello 
se define el modo dominante de vibration (lateral, vertical, torsional) asociado a la frecuencia pico. 

d) Se dibujan las Funciones de Coherencia entre las senales analizadas con el objeto de estimar la 
magnitud de los posibles ruidos registrados durante el proceso de medicion de las senales. 

e) Conocidas las frecuencias de interes, se precede a calcular la forma relativa de los modos de 
vibration con ayuda de la expresion (Bendat and Piersol, 1961) 

<Pi (y i ) = - F y iyi ( f ,) (9) 



Donde F es la transformada de la funcion de autocorrelation para la frecuencia fi 

asociada al punto de medicion j. Esta ultima expresion equivale a dibujar el espectro de potencia, 
para el modo de interes, en diferentes puntos de medicion de la estructura. (yj) y despues unir los 
picos que correspondan a las frecuencias naturales detectadas (fi). 



2.3 Matriz dinamica de flexibilidad 



En este apartado se desarrolla un modelo para determinar la matriz dinamica de 
flexibilidad de una estructura reticular sometida a vibration forzada. Se acepta que el sistema 
estructural es lineal, elastico y no amortiguado. En consecuencia, por definition de modo de 
vibration (Clough and Penzien, 1975) 

1k]- Pi 2 [m]}(u 1 }={0} (10) 

En donde [K] es la matriz estatica de rigidez, pi es la frecuencia natural de vibration (pi = 
2jtfi), [M] es la matriz de masas, {Ui} es el vector modal asociado a pi. A partir de la Expresion 10 
puede demostrarse que (Ruiz Sibaja, 1991) 

{u i } T [k]{u j }= p j 2 {u ; } t [m]{u j } 



( 11 ) 
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Por otra parte para una estructura sometida a vibration forzada la ecuacion de 
movimiento es (Paz, 1997) 



{[k] - CO 2 [m]]{x o } = {F 0 } (12) 

En la Expresion 12 {xo} es el vector de amplitudes de desplazamiento, co es la frecuencia de 
excitation y {Fo} es el vector de amplitudes de la carga de excitation. La amplitud de la respuesta 
dinamica se obtiene como (Clough and Penzien, 1975): 

{ x o}=Z“i{U i } (!3) 

i=l 

En la cual oti es la amplitud modal y {Ui} el vector modal. A1 sustituir la Ecuacion 13 en la 
Ecuacion 12 y considerar la Ecuacion 11 se obtiene la expresion para la amplitud modal (Lopez 
Rios, 1990) 



a _ {u.I t {f u } 1 
1 {u j } T [m]{u i } (p j 2 -G) 2 ) 

A1 incorporar la Ecuacion 14 en la Ecuacion 13 resulta 



{x 0 }=Z 

i=l 




1 1 

-ctiLLMu,} 



Mujfo} 



Se establecen las siguientes definiciones 

B-P^-CO 2 



(14) 



(15) 



(16) 



ft i ={U i } T [MKu i } 



(17) 



Se obtiene, a partir de la Ecuacion 15, la expresion para la matriz dinamica de flexibilidad 

(18) 

i=i L> ; l>2 j 

En esta expresion es necesario conocer la frecuencia de excitation, la frecuencia natural de 
vibration y los vectores modales para determinar la matriz dinamica de flexibilidad de una 
estructura reticular, [Fdin]. 



2.3 Matriz dinamica de rigidez 

A partir de la Ecuacion 13 se puede demostrar que (Ruiz Sibaja, 1991) 

{U 1 } T [M]{x 0 }=a i {U 1 } T [M]{U i } (19) 

En esta expresion puede incorporarse la expresion para la amplitud modal, Ecuacion 14, 
de lo cual resulta 
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{U,} T [MKx 0 }=|^j (20) 

Que tambien puede escribirse 

B,{U 1 } T [M]{x 0 }={u j } T {Fo} (21) 

A1 premultiplicar por [M] y luego por {Ui} la Ecuacion 21 puede llevarse a la forma (Ruiz 
Sibaja, 1991) 



[mKu ; ) 1 B, {u, F [M]{x 0 } = {F 0 } (22) 

E incorporando la Ecuacion 12 en la Ecuacion 22 y considerando la participation de todos 
los modos de vibration resulta 

Z MUi K«i )"' B, {uJ T [M] = [k - C0 2 m] ( 23 ) 

i=l 

Finalmente 

[K din ] = t r B, {u, } T [M] ( 24 ) 

i=i 

Mediante la Ecuacion 24 se puede obtener la matriz dinamica de rigidez, [Kdin], de una 
estructura reticular. A1 igual que en la Ecuacion 18 es necesario conocer la frecuencia de 
excitation, la frecuencia natural de vibration, la matriz de masas y los vectores modales para 
determinar [Kdin]. 

3. CALIBRACI6N del modelo desarrollado 



3. 1 Modos de vibration medidos 



E=2,100,000kg'orr2 
W/d =7.85x10-3 kgfcrrfT j 
A= 5.876 or2 





I =0.6253 or4 



Wl v\2 v\8 v\4 v\5 




wl =4.042 hg 
\A2=4.176hg 
W3 =4.1 12 kg 
v\4=4.100kg 
vv5=4.113kg 







Axtadonesencm 



Figura 1 Modelo estructural experimental. 
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A fin de realizar una calibration preliminar del modelo desarrollado previamente, se 
llevaron a cabo pruebas de vibration forzada sobre una viga de acero doblemente articulada 
(figura l). Esta viga, ademas de su propio peso, carga cinco pesos concentrados tambien de acero. 
La vibration forzada de la viga se consiguio mediante un motor electrico de 117 V, el cual hace 
girar un vastago con un peso de 20 grs. concentrado en el extremo del mismo. En la prueba initial 
la frecuencia de excitation se establecio en 12.20 Hz. 

Se utilizo un equipo de medicion a base de dos acelerometros piezoelectricos y un 
analizador de espectros con dos canales de entrada. Los acelerometros utilizados tienen una gama 
muy extensa de frecuencias (de 1 a 1000 Hz) y muy buena linealidad en todas ellas (Endevco, 
1971). Estos acelerometros se conectaron a cables blindados para enviar la serial de vibration al 
analizador de espectros. El analizador que se uso cubre el intervalo de frecuencias de 0.02 Hz a 25 
Khz. (Hewlett Packard, 1979). A la serial registrada en cada uno de los acelerometros se le hizo un 
muestreo y se proceso mediante un filtro digital antes de almacenarla en la memoria del aparato. 
Estos datos almacenados en la memoria se analizaron mediante la Transformada Rapida de 
Fourier. 

El procedimiento de medicion consistio en analizar dos senales simultaneamente. Uno de 
los acelerometros se coloco fijo como referenda (canal A), mientras que el otro se situo en puntos 
estrategicos sobre la viga (canal B), con el fin de definir la respuesta dinamica a medir. 




De esta manera se seleccionaron cuatro posiciones diferentes de los acelerometros (figura 
2), para realizar series de mediciones sobre la viga de acero. En las figuras 3 y 4 se muestran, 
para la primera position de medicion, los espectros de transferencia de fase y de coherencia en el 
dominio de la frecuencia, en el intervalo de 0 a 100 Hz, despues de efectuar 32 promedios RMS. 
En todos los casos se tomo como referenda el acelerometro del canal A. 

Se observa que en los 6 Hz el valor de la funcion de transferencia de fase (Figura 3) es 
cercano a 0° (-4°), lo cual indica que las senales de los canales A y B estan en fase en esa 
frecuencia. Ademas el valor de la funcion de coherencia (Figura 4) cercano a 1.0 (0.94) senala que 
se trata de la misma respuesta dinamica (en este caso el primer modo de vibration). 
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Figura 3 Funcion de transferencia de fase para la primera posicion de medicion. 



En 24.10 Hz el valor de la funcion de transferencia de fase (figura 3) es de 180°, indicando 
asi que la serial registrada por el acelerometro del canal B esta desfasada en 180° con respecto a 
la senal del canal A. El valor de 0.93 de la funcion de coherencia (Figura 4) confirma que se trata 
del segundo modo de vibracion. 




FRECUENCIA (Hz_) 

Figura 4 Funcion de coherencia para la primera posicion de medicion. 



Identification de estructuras reticulares mediante 
datos de mediciones de vibration forzada 



169 



En la tabla 1 se resumen los valores obtenidos, a partir de los espectros en el dominio de 
la frecuencia, para los cinco primeros modos de vibration que pudieron medirse en la viga. De 
acuerdo con los datos de la tabla 1 se pueden formar los vectores caracteristicos que corresponden 
a los modos de vibration de la estructura, esto se consigue al considerar las amplitudes 
registradas en cada position y el comportamiento de las funciones de transferencia de fase y de 
coherencia. 



3.2 Matriz dinamica de flexibilidad medida 



En la figura 1 se resumen las caracteristicas geometricas y de material del modelo 
estructural experimental utilizado durante las pruebas de vibration forzada. De acuerdo con los 
datos de esta figura se calculo una matriz teorica de masas concentradas (ecuacion 29). Mediante 
la aplicacion de la ecuacion 18 se obtuvo la matriz dinamica de flexibilidad medida. En el calculo 
se considero la participation de los modos medidos de vibration 1,2, 3, 4 y 5 y de la matriz de 
masas teorica, el resultado se muestra en la ecuacion 25 



Tabla 1 Resumen de resultados de mediciones de vibration 



Frecuencia 

(Hz) 


Amplitud 

(mV) 


Transferencia 

Amplitud 


Transferencia 

Fase 


Coherencia 


Canal A 


Canal B 


Primera 

Position 




6.00 


15.40 


8.79 


0.571 


- 4 ° 


0.94 


12.20* 


29.86 


31.63 


1.059 


0° 


0.98 


24.10 


10.22 


9.37 


0.917 


-180° 


0.93 


54.60 


1.29 


20.16 


15.625 


182° 


0.92 


93.40 


8.44 


7.74 


0.917 


-9° 


0.90 


133.80 


8.52 


5.00 


0.587 


165° 


0.76 


Segunda 

Position 




6.00 


15.40 


15.49 


1.006 


-4° 


0.98 


12.20* 


29.86 


49.94 


1.672 


2° 


0.97 


24.10 


10.22 


10.31 


1.009 


-183° 


0.91 


54.60 


1.29 


1.29 


1.000 


-175° 


0.91 


93.40 


8.44 


7.82 


0.927 


-171° 


0.88 


133.80 


8.52 


8.25 


0.968 


O 

T— 1 


0.74 


Tercera 

Position 




6.00 


15.40 


17.66 


1.147 


-3° 


0.99 


12.20* 


29.86 


46.90 


1.571 


2° 


0.98 


24.10 


10.22 


0.63 


0.062 


183° 


0.91 


54.60 


1.29 


20.96 


16.250 


-6° 


0.91 


93.40 


8.44 


0.54 


0.064 


170° 


0.87 


133.80 


8.52 


9.76 


1.145 


167° 


0.76 


Cuarta 

Position 




6.00 


15.40 


8.82 


0.573 


4° 


0.99 


12.20* 


29.86 


22.98 


0.770 


1° 


0.98 


24.10 


10.22 


9.35 


0.915 


3° 


0.92 


54.60 


1.29 


18.14 


14.063 


-171° 


0.88 


93.40 


8.44 


7.72 


0.915 


•170° 


0.88 


133.80 


8.52 


4.82 


0.566 


166° 


0.76 


* frecuencia de excitation 
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[Fdinl: 



1 

-0.264 



SIMETRICA 



-9.475 



La norma Frobenius de esta matriz es 



[^din ]l, 2,3,4 



-9.306 


-11.058 


-7.039 


1 




-15.383 


-17.424 


-11.059 


2 




-16.731 


-15.698 


-9.658 


3 xlO' 3 


(25) 




-9.357 


-4.503 


4 








-0.444 


5 




ZC=5 


(N 

1 

o 

X 

o 
cn 
i r) 






(26) 



3.2 Matriz dinamica de rigidez medida 



Mediante la aplicacion de la expresion 24 se obtiene, considerando la participacion de los 
modos medidos de vibracion 1, 2, 3, 4 y 5 y de la matriz de masas teorica 



Kd, 



1 

746.728 



4 



-739.631 303.722 -100.454 

1095.434 -834.956 333.080 

1155.042 



SIMETRICA 



5 

40.133 ' 
-83.104 
-877.672 353.329 

1159.323 -783.481 
722.416 



(27) 



Cuya norma Frobenius es 



k,d,:,„ L= XV =3302.933 



(28) 



3.3 Matrices de flexibilidad y rigidez teoricas 



Se llevo a cabo el analisis modal del modelo estructural para encontrar las frecuencias 
naturales teoricas de vibracion. De esta manera se obtuvieron los modos naturales de vibracion 
cuyas frecuencias respectivas son pi = 6.172 Hz, p 2 = 24.698 Hz, p3 = 55.355 Hz, pi = 95.664 Hz y 
ps = 136.678 Hz. A1 considerar la matriz de masas teorica (Ecuacion 29), los modos de vibracion y 
sus frecuencias, se obtuvieron las matrices dinamicas de flexibilidad y rigidez teoricas, [F’din] y 
[K’dirJ. 



[M] = 



1 

4.120 

0 

0 

0 

0 



2 

0 

4.257 

0 

0 

0 



3 

0 

0 

4.192 

0 

0 



4 

0 

0 

0 

4.179 

0 



5 

0 

0 

0 

0 

4.193 



xlO 



-3 



(29) 
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A estas matrices se les calculo la norma Frobenius para comparar el resultado contra lo 
que se presenta en las Ecuaciones 26 y 28. En las Ecuaciones 29 y 30, se muestran los valores 
obtenidos para | | [F'din] I 1 2 , y I I [K'din] I I 2 , respectivamente. 



1 [FLli.2.3.4,5 || 2 = X f ij 2 =5.607 xlCT 2 




(30) 


II [K'J, ,«,« || 2 = ./!>,/ =3447.673 




(31) 


A1 comparar los valores dados en las ecuaciones 26 y 28 (| | 
contra lo que se presenta en las Ecuaciones 30 y 32 ( | | [F'diJ I I 2 , y I 


[FdiJ I 
1 [K'diJ 


1 2 , y II [KdiJ I 1 2 ), 
I | 2 ) se observa una 



buena aproximacion entre lo que se obtiene mediante mediciones de vibration forzada y lo 
pronosticado por la teoria. 



3.3 Comportamiento del modelo para diferentes frecuencias de excitation 



Para realizar una calibration mas amplia del modelo de analisis propuesto se estudio su 
respuesta ante diferentes frecuencias de excitation. Se aplicaron seis diferentes frecuencias de 
excitation: 01 = 76.655 rad/seg, 02 = 219.597 rad/seg, 03 = 439.195 rad/seg, 0)4 = 658.792 rad/seg, 
cos = 878.389 rad/seg y 0)6 = 1097.987 rad/seg. Para cada prueba de vibration forzada se obtuvieron 
las matrices [FdiJ y [Kdin], es decir, se determinaron seis matrices [FdiJ y seis matrices [KdiJ. Se 
calculo la norma Frobenius para cada una de estas matrices y los valores obtenidos se escalaron 
respecto de los mayores valores calculados para | | [Fain] I I 2 y I I [Kdin] I I 2 , respectivamente, de 
forma tal que se fijo una escala horizontal entre 0 y 1. 

En la figura 5 se muestra el comportamiento en funcion de la participation modal, de la 
norma | | [FdiJ I I 2 escalada, para diferentes frecuencias de excitation. 

A partir de la figura 5 pueden considerarse los siguientes intervalos de comportamiento de 

I I [FdiJ I I 2 



°<“i <^(Pi+P 5 ) 


(32) 


(Pi +p 5 ) < tOi < -(Pi +Ps) 


(33) 


(Pi + P 5 )< W i <^(Pi + Ps) 


(34) 



Para el primer intervalo se observa que | | [FdiJ I 1 2 tiende a mantenerse constante a 
partir de la participation del 60% de la matriz modal (los tres primeros modos de vibration). Esto 
significa que para este intervalo de excitation, la definition de [FdiJ no presenta variaciones 
apreciables a partir de la participation del tercer modo de vibration. En el segundo intervalo, el 
comportamiento de | | [FdiJ I 1 2 tiende a mantenerse constante a partir de la participation del 
80% de la matriz modal, es decir, que se requiere de los cuatro primeros modos de vibration para 
que la definition de [FdiJ no presente variaciones notorias. Conforme la frecuencia de excitation 
incursiona hacia valores superiores se hace necesaria la participation del 100% de la matriz 
modal, esto se advierte en el tercer intervalo, en el cual se necesitan de los cinco modos de 
vibration para que [FdiJ alcance una buena definition. 
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Figura 5 Comportamiento de | | [Fain] I 1 2 para diferentes frecuencias de excitacion. 
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Figura 6 Comportamiento de | | [Kdinl I I 2 para diferentes frecuencias de excitacion. 
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En la figura 6 se muestra el comportamiento en funcion de la participation modal, de la 
norma | | [Kdin] I I 2 escalada, para diferentes frecuencias de excitation. Pueden considerarse los 
siguientes intervalos de comportamiento de | | [Kdin] I I 2 : 



Para el primer intervalo, se advierte que es necesario contar con la participation del 100% 
de la matriz modal para alcanzar una buena definition de [Kdin]. En el segundo intervalo, el 
comportamiento de | | [Kdin] | 1 2 se mantiene casi constante a partir de la participation del 80% 
de la matriz modal, por lo tanto, solo se necesitan cuatro modos de vibration para que la 
definition de [Kdin] no presente variaciones apreciables. 



Se presento una metodologia para deducir las matrices dinamicas de flexibilidad y rigidez 
a partir de datos de mediciones de vibration forzada. En esta etapa de la investigation, es 
necesario proponer la matriz de masas del sistema por no haber sido posible deducirla 
experimentalmente. Aunque el estudio de un modelo de laboratorio con cinco masas concentradas 
no justifica una generalization, los resultados obtenidos con la metodologia propuesta muestran 
que es posible determinar, con una aproximacion adecuada, las matrices dinamicas de flexibilidad 
y rigidez, para estructuras sometidas a vibration forzada. 

En las pruebas realizadas lograron identificarse cinco modos de vibration en flexion del 
modelo estructural experimental. Las calibraciones del modelo matematico muestran que se 
obtiene una mejor definition de [Fdin] cuando la frecuencia de excitation se encuentra en la 
vecindad de las primeras frecuencias naturales de vibration. Por el contrario, se observo una 
mejor aproximacion hacia [Kdin] cuando la frecuencia de excitation incursiona en el intervalo de 
las frecuencias naturales de vibration superiores. Los resultados obtenidos muestran, en general, 
que puede obtenerse una buena aproximacion de [Fdin] cuando se conoce el 60% de la matriz 
modal. En cambio, la definition de [Kdin] requiere del conocimiento del 80% de la matriz modal. Se 
advirtio una gran similitud entre lo que se obtiene mediante mediciones de vibration forzada y lo 
pronosticado por la teoria, probablemente porque se trata de una estructura sencilla. 
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RESUMEN 



Se hacen comentarios sobre los llamados criterios de diseno sismico por desempeno, 
indicando que partes de estos criterios, a juicio del autor, constituyen aportaciones novedosas, y 
que partes son ideas antiguas con nuevos nombres. Se hace ver que ciertas porciones de los 
criterios de diseno por desempeno han estado incluidas desde hace muchos anos en las normas 
convencionales. Posteriormente, se senalan algunos inconvenientes practicos que hay en la 
aplicacion de los criterios de diseno por desempeno, muchos de los cuales no son exclusivos de este 
nuevo enfoque, sino que son, desde hace tiempo, limitaciones del diseno sismico a secas. 
Finalmente, se indican algunos de los puntos en que las normas mexicanas estan actualmente 
anticuadas y se esbozan algunos caminos para mejorarlas, no necesariamente en el marco del 
diseno por desempeno sino en el marco del diseno racional. 

1. INTRODUCTION 

Las normas de diseno por sismo en Mexico tienen mas de 50 anos de antiguedad. En 
varios momentos de su historia, las normas mexicanas han aportado ideas y practicas novedosas, 
algunas de las cuales han sido adoptadas despues en normas de otros paises (Fukuta, 1991): en 
1942, los factores de importancia! en 1957, la distribution de fuerzas sismicas con la altura del 
edificio, el metodo dinamico de analisis, la limitation de desplazamiento lateral en las estructuras 
y la obligatoriedad de uso de coeficientes mayores en la zona de lagoi en 1976, factores novedosos 
de reduction por ductilidad. 

En la actualidad, sin embargo, muchas de las normas mexicanas son relativamente 
antiguas, por lo que no reflejan todos los avances recientes en la materia. 

En el presente trabajo se hace una breve reflexion sobre el estado actual de la 
normatividad en nuestro pais y sobre que modificaciones serian pertinentes a la luz de los 
llamados criterios de diseno por desempeno. 

Pero antes de discutir sobre las normas mexicanas, a continuation se hacen algunos 
comentarios sobre los propios criterios de desempeno, con el fin de hacer ver algunos 
inconvenientes practicos para su aplicacion y con la intention de senalar que algunos de estos 
criterios han estado presentes en las normas mexicanas desde hace mucho tiempo. 

2. QUE NO ES NUEVO EN EL DISENO POR DESEMPENO 

2.1. Tener objetivos de diseno 



Diversos autores han senalado ( e.g :, Hadjian, 2002) que el diseno siempre ha tenido 
objetivos de desempeno. Las estructuras siempre se han diseriado para que “hagan algo”, es decir, 
para que tengan un cierto desempeno que alguien ha visualizado de antemano. Por esta razon, el 
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nombre “diseno por desempeno” es desafortunado y para muchos fines podria sustituirse 
simplemente por “diseno”. 

Los objetivos de diseno no siempre se han logrado! muchas estructuras han fallado a lo 
largo de la historia. Pero que no se hayan logrado siempre no se debe a que no hayan existido. Se 
debe a ignorancia, dolo, mala suerte, pero no a ausencia de objetivos de diseno. Esto nos conduce a 
una cuestion que el autor considera central- siempre han existido objetivos de desempeno! lo que 
ha fallado es el como. La diferencia principal esta entre el diseno bueno y el diseno malo. 



2.2 El concepto de “nivel de desempeno” 



Una lectura somera (o profunda) de los textos seminales del diseno por desempeno nos 
hace ver que “nivel de desempeno” es un sinonimo perfecto de “estado limite”. No hare mas 
comentarios, salvo insistir en lo innecesario e inconveniente que es llamar a la misma cosa de dos 
maneras distintas. Propongo que mantengamos el acostumbrado “estado limite” y evitemos 
confusiones, aun a costa de perder cierta sensation de modernidad y avance. 

2.3 Diseno multinivel o diseno para varios estados limite 



A continuation se reproducen algunos articulos del Reglamento para las Construcciones 
del DF vigente en la actualidad. Estos dos articulos, en esencia, han estado presentes en los 
reglamentos del DF desde 1966. Las cursivas son del autor: 

> Art. 172 (Titulo Sexto del Reglamento de Construcciones del DF). “Este titulo contiene los 
requisitos que deben cumplirse... para lograr un nivel de seguridad aceptable i^que ocurra 
poco frecuentemente /) contra fallas estructurales Qcolapso?) asi como un comportamiento 
estructural aceptable Qoperatividad?) en condiciones normales de operation Qfrecuente?) 

> Art. 182 “Toda estructura y cada una de sus partes deberan disenarse para cumplir con 
los requisitos basicos siguientes: 

> I Tener seguridad adecuada Qfallas poco frecuentes?) contra la aparicion de todo estado 
limite de falla posible Qcolapso?) ante las combinaciones de acciones mas desfavorables 
que puedan presentarse durante su vida esperada, y 

> II No rebasar ningun estado limite de servicio Qoperatividad?) ante combinaciones de 
acciones que correspondan a condiciones normales de operation Qcargas frecuentes?)' 

Como se observa, la necesidad de que un diseno satisfaga simultaneamente varios estados 
limite asociados a diferentes frecuencias de ocurrencia ha sido desde hace muchos anos una 
aspiration explicita de los criterios de diseno. Hadjian (2001) hace ver algo parecido, aunque mas 
dramatico: la aspiration explicita del documento Vision 2000 (Hamburger et a 1, 1995), 
considerado la piedra de toque del diseno por desempeno, es practicamente igual a la expresada 
por la SEAOC en 1967. 

Nuevamente, el problema ha sido el como, no la ausencia del deseo de disehar para 
satisfacer simultaneamente varios estados limite. 



2.4 Disenos mas conservadores para estructuras importantes 

Desde hace mucho tiempo las estructuras mas importantes han sido disenadas para 
satisfacer estados limites mas estrictos. Los llamados factores de importancia (los que se aplican a 
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estructuras del Grupo A, del Grupo B o del Grupo C) estan presentes en los reglamentos 
mexicanos por lo menos desde 1966. 

No esta claro que el uso de factores de importancia conduzca siempre a disenos adecuados! el uso 
de un factor de importancia independiente de la sismicidad de la zona parece tener poco sentido. 
Es, nuevamente, un problema de ignorancia sobre el fenomeno sismico y dificultad para hacer 
muchos calculos, pero no una falla esencial de los criterios de diseno. 



3. CONTRIBUCIONES REALES DEL LLAMADO DISENO POR 
DESEMPENO 

3. 1 El deseo de usar mejores indicadores de desempeno 



Ha sido muy importante para el progreso del diseno por sismo la insistencia, 
relativamente reciente, sobre usar mejores indicadores del desempeno estructural. En particular, 
la difusion del uso de los desplazamientos como medidas de la respuesta ha sido muy benefica. 
Esto esta relacionado con un deseo, correcto desde mi punto de vista, de saber en que condiciones 
va a quedar la estructura despues de sismos que ocurren con frecuencias dadas. 



3.2 El recordatorio de que satisfacer el estado limite de colapso no lleva a satisfacer 
todos los demas. 



Aunque el diseno multinivel —satisfacer varios estados limite simultaneamente - ha 
existido desde hace mucho tiempo, tambien ha existido la perception de que si una estructura 
puede lo mas, tambien puede lo menos. Es facil demostrar que esto, en muchos casos, no es cierto- 
una estructura puede tener un comportamiento aceptahle ante el sismo de colapso y un 
desempeno pohre ante sismos frecuentes de servicio! las estructuras flexibles de la ciudad de 
Mexico son un ejemplo notable de esta situation. Por ello, la insistencia en disenar para, 
realmente, satisfacer varios estados limite es muy pertinente. 



3.3 El senalamiento de las limitaciones del uso de factores de importancia 



Se ha hecho ver que solo con incrementar las cargas de diseno para estructuras mas 
importantes no necesariamente se consigue un comportamiento satisfactorio. En efecto, para 
estructuras importantes no basta con incrementar las cargas de diseno- es necesario especificar 
detalladamente los estados limite deseados asociados al desempeno deseado. El diseno por 
desempeno ofrece un excelente marco para resolver esta limitation de los factores de importancia. 



3.4 Tasas de excedencia versus tasas de falla 



En el estado actual del diseno por desempeno, se especifica la frecuencia de ocurrencia de 
las demandas asociadas a los estados limites considerados. Como se vera a continuation, esto no 
garantiza que se llegue a estructuras de confiabilidad uniforme, es decir, estructuras cuyas fallas 
acontezcan con la misma frecuencia promedio. Esto se debe a que, en general, las fallas pueden 
ocurrir aun con temblores de demanda nominal menor a la especificada en el codigo y, 
eventualmente, no ocurrir con temblores con demanda nominal superior. 
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Para ilustrar lo anterior proponemos el siguiente ejemplo. Supongamos que se tiene un 
tipo de estructuras que se disenaran para el estado limite a=c, donde a es la demanda y c la 
capacidad. Sea v A (a) la tasa de excedencia de la demanda, es decir, el numero medio de veces en 
que, por unidad de tiempo, la demanda excedera un valor particular, a. Supongamos que Va(ci) 
tiene la siguiente forma funcional (Esteva, 1970): 

V A (a) = Ka' (1) 

donde los parametros K y r (K> 0, r<0) dependen de la sismicidad regional. Se desea que este 
estado limite no sea excedido mas de, digamos, Vo veces por ano, lo cual se supone puede lograrse 
escogiendo como capacidad de diseno una cuyo valor numerico es igual al de la demanda asociada 
a una tasa de excedencia igual a v 0 (ver figura l). Llamaremos co a esa capacidad de diseno. 

Si la capacidad real de la estructura fuera determinista e igual a co, entonces, en efecto, la 
estructura fallaria cuando a>co, lo cual aconteceria exactamente v 0 veces por ano. En otras 
palabras, la tasa de falla de la estructura seria precisamente la tasa de excedencia de la demanda 
asociada al estado limite. 






demanda 




—□—Tasa de falla 





Figura 1 . Tasa de excedencia de la demanda y tasa de falla para un ejemplo hipotetico 



Sin embargo, por diversas razones, la capacidad estructural no es determinista, sino que 
existen incertidumbres en su calculo. Por sencillez, modelaremos esta incertidumbre de la 
siguiente manera: supondremos que la capacidad real, c, es una variable aleatoria con 
distribucion lognormal, con mediana igual a la capacidad nominal de diseno, co, multiplicada por 
un factor de sobrerresistencia, F, y una cierta desviacion estandar del logaritmo natural de c igual 
a CJ. La aplicacion del teorema de la probabilidad total nos lleva a que la tasa de falla cuando se 
adopta para diseno una capacidad nominal Co, Vc(co), esta da da por: 



Vc( c o ) ~ 



-J Pr( c<a\c Q ) 



dvja ) 

da 



da 



( 2 ) 



donde Pr(c<a / co) es la probabilidad de que la capacidad (real) sea menor que la demanda dado 
que se adopto para diseno una capacidad nominal co. Con la hipotesis simplificatorias 
mencionadas, esta probabilidad se determina de la siguiente manera 1 
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Pr(c < a\c 0 ) = j 

o 



1 

2noc 



— jfln — r 

2a ~ Fc ° dc 



(3) 



Sustituyendo las ecuaciones 3 y 1 en la 2, y resolviendo las integrales correspondientes, 
obtenemos que: 



v c (c g ) = K(Fc 0 ) 



l 2 2 1_ 

e 2 - v( Fc 0 )e 2 



= v 0 F r e 



-rV 

2 



(4) 



Se observa, como ya hablan senalado otros autores ( e.g :, Esteva y Ruiz, 1989) resolviendo 
casos mas realistas, que las tasas de falla pueden ser significativamente diferentes de las tasas de 
excedencia prefijadas, excepto en el caso determinista en que F=1 y <7=0, con lo cual Vc(co) = Vo- 

La ecuacion 4 ilustra varios hechos importantes. 

1) Las tasas de falla, las que verdaderamente importan, son diferentes de las tasas de 
excedencia de la demanda, que son las que actualmente se estipulan en las propuestas 
reglamentarias. 

2) La tasa de falla de una estructura en particular depende, aun en nuestro ejemplo 
simplificado, de su sismicidad regional. Dos estructuras iguales localizadas en diferentes 
partes del pais (r diferente, en nuestro ejemplo simplificado) tendran confiabilidades 
diferentes aun si ambas se disenan para una demanda asociada al mismo periodo de 
retorno. 

3) La tasa de falla depende de la incertidumbre que exista en el calculo de la capacidad 
estructural o, de manera mas general, de su vulnerabilidad. Dos estructuras construidas 
en el mismo sitio tendran tasas de falla diferentes si sus caracteristicas de vulnerabilidad 
(ejemplificadas con F y (7 en el ejemplo simplificado) son diferentes, aun cuando hayan 
sido disenadas ambas para la demanda asociada al mismo periodo de retorno. 

Esta situation no debera ser ignorada en la normatividad futura. La solution a secas es 
simple^ que los reglamentos estipulen parejas estado limite-tasa de falla y no parejas estado 
limite'tasa de excedencia de la demanda, como se hace en la actualidad. La solution practica, sin 
embargo, es sumamente compleja. 

En efecto, para determinar aceptablemente la tasa de falla de la estructura no bastaria 
con calcular su probabilidad de falla ante el sismo con tasa de excedencia Vo (la tasa de excedencia 
de la demanda de diseno) sino que, de acuerdo con la ecuacion 2, seria necesario determinar las 
probabilidades de falla condicionadas a la ocurrencia de varios niveles de demanda y, 
posteriormente, sumar ponderadamente dichas probabilidades, usando como factores de peso las 
frecuencias de ocurrencia de cada uno de los valores de demanda elegidos (En el fondo, la 
ecuacion 2 no es mas que una suma ponderada de probabilidades de falla). Esto es impracticable 
con las herramientas y los costos de hoy. 

Sin embargo, el resultado expresado en la ecuacion 4, aun tratandose de un caso 
simplificado, sugiere posibles caminos para acercarse al objetivo de diseno para confiabilidades 
predeterminadas y constantes: 

1) Utilizar factores de sobrerresistencia o factores de seguridad diferentes (el factor F en 
nuestro ejemplo) dependiendo del tipo de estructura y la sismicidad regional. 

2) Utilizar periodos de retorno diferentes para la demanda de diseno (diferentes Vo) 
dependiendo del tipo y localization de la estructura. 
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Esto, desafortunadamente, eliminaria la aparente sencillez y transparencia de los planteamientos 
conceptuales del llamado diseno por desempeno. El tema amerita reflexion. 



3.5 ^Tasas de falla uniformes? 



Hemos visto que, aun si se estipulan periodos de retorno uniformes a lo largo de un pais 
para las demandas de diseno, no necesariamente se consiguen estructuras con confiabilidad 
uniforme. Pero cabe la siguiente pregunta- nadie duda que es util disenar estructuras con 
confiabilidad conocida, pero £es deseable tener estructuras de diferentes tipos con confiabilidad 
uniforme a lo largo de regiones extensas?. Para responder a esta pregunta, que a mi juicio sigue 
abierta, se han esgrimido argumentos de muy diferente indole^ equidad, justicia, optimacion de 
costos, y otros. 

En opinion del autor, que una confiabilidad sea adecuada para una estructura en una 
parte del pais no implica que tambien lo sea para otra estructura en otro lado. 

Los analisis de optimacion de costos indicaron desde muy temprano en el desarrollo de la 
ingenieria sismica (Esteva, 1970) que la estrategia de confiabilidad uniforme conduce casi 
inevitablemente a disenos sub-optimos. La razon es clara- el costo de dar resistencia lateral a una 
estructura crece mas rapidamente que, digamos, el coeficiente adoptado para diseno. Esto 
significa que en las zonas de baja sismicidad, con niveles de diseno relativamente bajos, es mas 
barato “sobre'disenar” una estructura que en una zona de alta sismicidad, por lo cual uno se 
puede dar el lujo de ser mucho mas conservador. 

Ilustraremos esto siguiendo el planteamiento de Esteva (1970), hecho hace mas de 30 
anos. Supongamos que el costo de construccion, Cc, varia con la resistencia nominal ante cargas 
laterales, co, de la siguiente manera- 



C c (c 0 ) = 



K, 



{ K,+K 2 c 0 c 



sic 0 <c K 
si c 0 > c K 



(5) 



donde K 2 es una constante de normalization y Ki es el costo fijo de construccion que se tendria 
aun sin diseno ante cargas laterales. El exponente a es mayor que 1, indicando que el costo de 
suministrar resistencia lateral crece mas rapidamente que el valor de esta resistencia lateral. 

Por otra parte, puede demostrase que el costo esperado — actualizado a valor presente - de 
las perdidas por sismo, Cp, suponiendo que estas ocurren de acuerdo con un proceso de Poisson y 
que se utiliza para diseno una resistencia co, esta dada poL 

C P (c 0 )=D Vc(Co) (6) 

y 



donde D es el costo de la falla. El costo total de la decision de adoptar co como valor de diseno, Ct, 
es entonces la suma del costo de construccion mas el costo esperado de las perdidas por sismo. La 
mejor decision es la que minimiza el costo total, es decir, aquella para la cual 



dC T (c 0 ) _ q 
dc () 



(7) 



El coeficiente optimo de diseno, cop, resulta ser 
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KrD 

y K 1 ocF r 



Ira 



que resulta estar asociado a una tasa de falla, v 0 p , dada por 



( 8 ) 



^ op ~~ 



r K^ayK,^* 

rDF a 



1 r 2 a 2 a 

2 r+a 



(9) 



Como puede observarse, la tasa optima de falla depende de todos los parametros que 
hemos hecho interveniH sismicidad, vulnerabilidad estructural, costo de construction y costo de la 
falla. La tasa optima de diseno (el periodo de retorno de la falla) depende de la localization, el tipo 
estructural, el costo de construction y el costo de la falla. No hay, pues, tasas optimas de falla 
generales. 



Analicemos un caso simple : dos estructuras iguales (que tienen, por tanto, iguales F, <x, Kj, 
a y D) seran construidas en dos sitios diferentes, tales que en el sitio A la sismicidad es 
exactamente el doble que en el sitio B. (igual rpero diferente A) . Supongamos que en el sitio B —el 
de menor sismicidad- se ha calculado que 0.001/ano es la tasa optima de diseno (periodo de 



V opa V OPB 



mos). 


En 


ese caso 


la 




a 








r+a 




m 






= 0.001 










UJ 



( 10 ) 



Si <2=1.5 y r =-2.5 (valores tipicos de estos parametros), tendriamos que la tasa optima de 
diseno en el sitio A —el de mayor sismicidad- seria de 0.0028/ano, asociada a un periodo de retorno 
optimo para la falla de solo 353 anos. En otras palabras, para el sitio de mayor sismicidad no es 
optimo ser tan conservadores como para el sitio de menor sismicidad. 

Desde este punto de vista, no parece razonable, por ejemplo, que si se juzga adecuado 
disenar una estructura para 150 anos de periodo de retorno en Acapulco sea razonable tambien 
hacerlo para el mismo periodo de retorno en Aguascalientes o Ciudad Juarez. El analisis de costos 
— y en parte el sentido comun- indican que en Aguascalientes o Ciudad Juarez deberia aspirarse a 
mayores confiabilidades puesto que resulta barato. 

Quiza la estrategia mas adecuada sea imponer una confiabilidad minima, socialmente 
aceptable, para las zonas de alta sismicidad y aspirar a confiabilidades progresivamente mayores 
conformes la amenaza sismica vaya disminuyendo 1 . 



4. ESTADO ACTUAL DE LA NORMATIVIDAD EN MEXICO Y MEJORAS 
FUTURAS 



Como se sabe, la normatividad de la construction en Mexico es de competencia municipal. 
Por ello, en principio, cada uno de los cerca de 2500 municipios de Mexico podria tener sus propias 
normas de diseno y construction 1 . A esto hay que anadirle las normas que algunos sectores 
industriales desarrollan para sus propias obras (CFE, por ejemplo). 



1 Esto no ha ocurrido asi en la practica. Ha sido mas comun que se desarrollen normas de alcance estatal que 

posteriormente son adoptadas por todos los municipios del estado. 
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En el pasado se han hecho esfuerzos por desarrollar normas tipo que, adaptadas a la 
sismicidad y condiciones del suelo locales, constituyan la base de las normas estatales o 
municipales. Todos estos esfuerzos han fracasado. 

Existian hasta 2000 por lo menos 40 reglamentos diferentes en Mexico (Sanchez, 2000). 
En terminos generales, todas las normas mexicanas de diseno por sismo actualmente vigentes 
estan basadas en dos documentos principales 1 la Normas Tecnicas Complementarias para Diseno 
por Sismo del DF, y el Capitulo de Diseno por Sismo del Manual de Obras Civiles de la Comision 
Federal de Electricidad. Las normas del DF datan de 1985 (con actualizaciones en 1987 y 1992) 
mientras que las normas de CFE fueron emitidas en 1993. 

Para fines practicos, las normas de CFE fungen como una especie de reglamento nacional, 
puesto que, sin ser obligatorias, son usadas en muchas obras construidas en sitios para los que no 
se dispone de reglamento propio y, en algunos casos, aun en sitios donde si hay reglamento local. 
Puesto que no se preve un cambio drastico de esta situation a mediano plazo, hay que resaltar la 
gran importancia de que las normas de CFE se mantengan al dia. 

En los siguientes parrafos haremos mention de algunos inconvenientes que presentan la 
mayor parte de las normas mexicanas actualmente vigentes e iremos senalando algunas de las 
mejoras que deberian hacerse. 



5 NIVELES DE PELIGRO SISMICO ASOCIADOS A LOS ESPECTROS DE 
DISENO 



Salvo excepciones, los reglamentos mexicanos no especifican los niveles de peligro 
implicitos en los espectros de diseno. En otras palabras, rara vez se sabe a que periodo de retorno 
estan asociadas las ordenadas espectrales de diseno. El que los periodos de retorno no se 
especifiquen tiene el inconveniente de que el responsable del diseno o el dueno de la obra no 
saben, ni siquiera de manera aproximada, que tan frecuentemente se excederan los estados 
limite. Por esta razon, el disenador y el dueno tienen solo una vaga idea del desempeno futuro de 
la estructura. Pero la ausencia explicita de referencias a los periodos de retorno en los 
reglamentos quiza esconde un problema mayor: no tenemos criterios robustos y de aceptacion 
general para fijar esos periodos de retorno. 

Como se ha senalado en la introduction, el problema de asociar periodos de retorno a los 
estados limite no es trivial. Da la impresion de que no es razonable usar el mismo periodo de 
retorno para todo el pais ni para todos los tipos estructurales. Los acuerdos sobre este asunto que 
se han alcanzado en otros paises nos ayudan poco : que un periodo de retorno de 500 anos sea 
considerado adecuado para cierto estado limite en una region de los Estados Unidos y los tipos 
estructurales que ahi se construyen, no implica que tambien lo sea para todo Mexico, con sus 
propios tipos estructurales y un nivel de desarrollo economico muy diferente. 

A pesar de su complejidad, es necesario avanzar rapidamente en la solution de este 
problema. Si no, veremos en el futuro la proliferation de normas o estudios especificos de sitio con 
espectros asociados a muy variados niveles de riesgo tolerable. Esta situation traera confusion y, 
muy probablemente, uso sub-optimo del dinero. 

6. INCLUSION DE LOS EFECTOS DE SITIO 



Como se sabe, los efectos de amplification de las ondas sismicas producidos por la geologia 
local, tambien conocidos como efectos de sitio, son probablemente el principal factor de riesgo 
sismico para las construcciones. A pesar de ello, las normas mexicanas (con la exception de las del 
Distrito Federal) dan al asunto relativamente poca importancia. 
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A medida que la instrumentation sismica ha ido progresando, con lo que ha aumentado 
considerablemente el numero de acelerogramas registrados, ha ido quedando mas claro que 
clasificar los suelos en “firmes, intermedios y blandos” es decididamente inadecuado. 

En nuestro pais, el caso mejor estudiado es la Ciudad de Mexico. Quiza se trata de un caso 
extremo, pero las variaciones que hay en las caracteristicas del movimiento del suelo en unos 
cuantos cientos de metros es notable. Clasificar toda la zona del lago como “suelo blando” no 
ayuda a lograr mejores disenos para resistir temblores. 

Ademas de contar con muchos mas acelerogramas registrados, en la actualidad se dispone 
de mejores herramientas teoricas y metodos de campo relativamente baratos que permitirian 
mejorar nuestra capacidad para estimar movimientos fuertes en un mayor numero de sitios y con 
esto poder incluir en los reglamentos los efectos de sitio de mejor manera. 

En este sentido, podria avanzarse mucho llevando a cabo estudios de microzonificacion 
sismica de las ciudades mas expuestas en Mexico. En su mayoria, los estudios realizados hasta la 
fecha en nuestro pais son inadecuados, porque se limitan a clasificar los suelos segun los antiguos 
criterios (firme, intermedio y blando) y practicamente no dan indicaciones cuantitativas utiles 
para la estimation de movimiento fuerte. Algunos paises, como Colombia, han vuelto obligatoria 
la elaboration de estudios de microzonificacion sismica para ciudades con ciertas caracteristicas! 
quiza deberiamos seguir el ejemplo. 

A falta de buenos mapas de microzonificacion, las normas futuras podrian reglamentar el 
uso de espectros especificos de sitio y exigir su ejecucion para estructuras de cierta importancia. 



7 FORMA Y TAMANO DE LOS ESPECTROS DE DISENO 



En general, los espectros de diseno en las normas mexicanas no reflejan el verdadero 
tamano ni la forma real de los espectros de los temblores para los que se quiere disenar. Aun 
cuando las fuerzas y desplazamientos de diseno sean correctamente calculados, utilizando 
factores correctivos adecuados, el tener espectros de diseno que no se parezcan suficientemente a 
los de respuesta genera problemas de falta de claridad y perdida de perspectiva durante el 
proceso mismo de diseno. Algunos ejemplos- un espectro de diseno para suelo blando que tenga 
una meseta demasiado ancha, no transmite correctamente la idea de que hay que evitar que el 
periodo fundamental de la estructura este cerca del periodo predominante del suelo! o un espectro 
de aceleraciones que no tenga una caida en periodo largo compatible con las observaciones y la 
dinamica estructural, llevara a malas estimaciones de los desplazamientos. 

Un colofon de esta consideration es que las amplificaciones (con respecto al terreno firme) 
observadas por efectos de sitio en una ciudad deberian estar reflejadas en los espectros de diseno. 
En otras palabras, si se observa que, para cierta ordenada espectral las aceleraciones en sitios 
blandos son mucho mayores que en terreno firme, las ordenadas de diseno deberan ser tambien 
mucho mayores. 

Una de las razones que han hecho que en la practica mexicana los tamanos de los 
espectros de diseno no reflejen los tamanos de los espectros reales es que en las normas 
generalmente se han incluido reducciones a las fuerzas de diseno de las que el disenador no esta 
suficientemente advertido. Estas reducciones implicitas, que han tratado de tomar en cuenta 
efectos como la sobrerresistencia estructural, probablemente no llevan a disenos inseguros, pero si 
producen cierta falta de claridad que no es deseable. Por esta razon, en las normas futuras debera 
procurarse que todas las reducciones a los espectros elasticos de diseno se hagan explicitas, aun 
cuando los factores de reduction usados no puedan justificarse completamente con bases teoricas 
o empiricas. 
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8. DEFINICION MAS CLARA DE LOS ESTADOS LIMITE 



Todos los reglamentos mexicanos vigentes estipulan que las estructuras deben disenarse 
para satisfacer dos estados limite, generalmente llamados de colapso y de servicio. El estado 
llmite de colapso, con las salvedades que se han senalado anteriormente (periodo de retorno 
desconocido, formas espectrales inadecuadas, reducciones por sobrerresistencia implicitas), esta 
bien definido- la estructura debe tener propiedades de rigidez y resistencia tales que la demanda 
global de ductilidad sea inferior a la capacidad ductil de la estructura cuando ocurre un temblor 
caracterizado por el espectro de diseno. El estado limite de servicio, sin embargo, no esta 
claramente especificado. La revision requerida por los reglamentos mexicanos consiste en calcular 
los desplazamientos de la estructura ante las cargas correspondientes al estado limite de colapso 
y despues compararlos con unos desplazamientos permisibles, independientes del sistema 
estructura, que no son los que provocarian la falla de la estructura ni aquellos para los cuales se 
iniciaria el agrietamiento. Esta revision esta disenada para que los calculos sean simples y para 
que solo tenga que hacerse un analisis estructural; los desplazamientos permisibles estan 
calibrados para que el diseno lleve a rigideces adecuadas ante sismos de baja intensidad, cuyo 
tamano o periodo de retorno tampoco se especifica. 

El empleo de estos procedimientos de diseno lleva, generalmente, a estructuras 
razonablemente seguras. Sin embargo, es muy poco transparente y hace muy dificil que el 
disenador tenga una idea clara del desempeno de la estructura durante los sismos de diseno. 

Es deseable, entonces, que los estados limite se especifiquen con mayor claridad aun 
cuando esto suponga incrementos moderados en el trabajo de analisis estructural. En la nueva 
propuesta de Normas Tecnicas Complementarias para Diseno por Sismo del DF se ha avanzado 
en este aspecto. En dichas normas, se hace explicita la existencia de dos estados limite de 
desplazamiento, con distorsiones de entrepiso permisibles que reflejan mejor el desempeno 
estructural que se quiere obtener en los dos estados limites considerados. Es deseable que estas 
modificaciones se extiendan al resto de los reglamentos mexicanos. 

9. FACTORES DE REDUCCION POR RESISTENCIA 



Los criterios contemporaneos de diseno admiten que la estructura tenga incursiones en el 
rango no lineal ante el sismo que caracteriza el estado limite de colapso. Esto permite limitar las 
demandas de fuerza en los elementos estructurales (para, por tanto, utilizar resistencias de 
diseno menores), a costa de que se presenten demandas de ductilidad limitadas y cierto nivel de 
danos provocados por la fluencia de algunas secciones de la estructura. 

Para modelar el comportamiento no lineal, practicamente todos los reglamentos del 
mundo estan basados en analisis de un sistema de un grado de libertad con comportamiento 
elastoplastico. Es con este modelo que se determina la resistencia necesaria para limitar las 
demandas de ductilidad a un valor especificado, que llamaremos Q. Es usual expresar la 
resistencia necesaria para lograr una demanda de ductilidad dada, C(T, Q), como una fraction de 
la resistencia necesaria para tener una demanda de ductilidad unitaria (comportamiento 
elastico), C(T,lX notese que ambas resistencias dependen del periodo estructural, T. Llamaremos 
Q ’ a este cociente: 



Q\t,q) 



C(T, 1 ) 
C(T,Q) 



( 11 ) 



La forma de la funcion Q ha sido ampliamente estudiada en los ultimos anos (Krawinkler 
et al., 1992! Miranda, 1993; Miranda y Bertero, 1994). En particular, Ordaz y Perez Rocha (1998) 
observaron que, bajo circunstancias muy generales, Q’ depende del cociente entre el 
desplazamiento espectral y el desplazamiento maximo del suelo. Esto implica que, como 
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observaron Meli y Avila (1988) al analizar registros del sismo de Michoacan obtenidos en la zona 
del lago del Valle de Mexico, los factores Q dependen fuertemente de las caracteristicas del suelo 
en que se desplanta la estructura. En vista de esto, es claro que las reglas que actualmente se 
utilizan en los reglamentos mexicanos, basadas en las derivadas por Newmark y Hall, son 
inadecuadas para suelos blandos por lo que es necesario sustituirlas por otras mas realistas. 



10. FACTORES DE IMPORTANCE PARA ESTRUCTURAS DEL GRUPO A 



Los reglamentos mexicanos vigentes clasifican a las estructuras en grupos dependiendo de 
su importancia. La mayor parte de los requisites de diseno son los mismos para todos los grupos, 
excepto que se multiplican por un factor de 1.5 las fuerzas de diseno cuando se trata de 
estructuras importantes, generalmente llamadas del grupo A. 

Esto tiene al menos dos inconvenientes. l) El incremento en confiabilidad estructural que 
se logra al disenar para fuerzas mayores depende de la sismicidad regional. En otras palabras, 
dependiendo del regimen sismico del sitio, un incremento de 50% en las fuerzas de diseno puede 
significar un incremento de confiabilidad mucho mayor que para la misma estructura localizada 
en otra parte. En general, seria preferible estipular una confiabilidad objetivo para estructuras 
del grupo A y despues ver que fuerzas de diseno deberian usarse para conseguir esa confiabilidad. 
Esto, como se ha hecho ver anteriormente, conduciria a factores de importancia variables por tipo 
estructural, tipo de suelo y localization geografica de la estructura. 2) No esta claro que en todos 
los casos se consiga un desempeno satisfactorio de una estructura del grupo A mediante el solo 
incremento de las fuerzas de diseno. Seria deseable, ademas de especificar una confiabilidad 
superior, especificar claramente objetivos especiales de desempeno para estructuras importantes, 
tal como lo senala el documento Vision 2000. 



11. CONCLUSIONES 



La normatividad en materia de diseno por sismo en Mexico presenta dos problemas 
principales. Por una parte, la ausencia de uniformidad en los criterios adoptados en las distintas 
normas puede conducir paulatinamente a estructuras con confiabilidades muy diferentes 
dependiendo de las normas que se hayan empleado en su diseno. Para combatir este problema es 
necesario que los diferentes responsables de la elaboration de normas (municipios, gobiernos 
estatales, empresas paraestatales, grupos universitarios) acuerden criterios generales de diseno 
comunes, alrededor de los cuales podrian elaborarse las normas locales. Por otra parte, ciertos 
criterios basicos de diseno necesitan cambios, drasticos en algunos casos, para poner las normas 
al dia e incluir los avances que recientemente se han tenido en diferentes aspectos del diseno por 
sismo. Se han mencionado en este trabajo los cambios que, a juicio del autor, son prioritarios. En 
algunos casos se han senalado posibles soluciones! en otros, sin embargo, solo se ha planteado el 
problema. 
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RESUMEN 



Despues de un terremoto fuerte, los danos en el area afectada son tan generalizados que 
no es posible que todas las evaluaciones de dano de los edificios sean hechas por ingenieros 
expertos. Por esto profesionales voluntaries con poca o ninguna experiencia en evaluation del 
dano sismico deben evaluar gran parte de los edificios afectados. Los evaluadores inexpertos 
presentan la tendencia a agravar o subestimar el nivel real del dano. En estos casos la 
information es subjetiva e incompleta y se utilizan calificaciones linguisticas. Por lo anterior, se 
desarrollo un sistema experto para la evaluation post-sismica del dano en edificios, utilizando 
tecnicas como las redes neuronales artificiales y la teoria de los conjuntos difusos. Esta 
herramienta facilita la participation de inexpertos en el proceso de evaluation. Este modelo de 
Inteligencia Artificial tiene en cuenta los posibles niveles de danos en elementos estructurales y 
no estructurales, los efectos del sismo en el suelo, asi como las condiciones preexistentes que 
pueden hacer mas vulnerable el edificio, como la calidad de los materiales de construction, 
irregularidades en planta y altura y una mala configuration estructural. El sistema decide sobre 
la habitabilidad y reparabilidad del edificio aplicando una base de reglas difusas a la information 
obtenida del edificio. 

1 INTRODUCTION 



En el caso de un sismo fuerte, debido a la magnitud del desastre, es necesario que el 
proceso de evaluation del dano en los edificios sea hecho por profesionales con experiencia en esta 
tarea y que esten familiarizadas con la ingenieria estructural. Desgraciadamente el personal que 
cumple con estas caracteristicas es muy poco, por lo cual en casos de dano generalizado en un 
centro urbano, es necesario involucrar ingenieros y arquitectos inexpertos. Como consecuencia, es 
comun la sobreestimacion o subestimacion del dano, llevando a la ocupacion peligrosa o la 
demolition innecesaria de las edificaciones. Motivado en esta problematica, este trabajo propone 
un sistema experto para apoyar el proceso de evaluation de danos, usando herramientas como las 
redes neuronales artificiales y la teoria de conjuntos difusos. 

2 EVALUACION DEL DANO DESPUES DE UN TERREMOTO 

Como resultado de los sismos ocurridos en diferentes paises ubicados en zonas sismicas, 
ha sido necesario que estos desarrollen manuales para la evaluation del dano en edificios, con el 
proposito de decidir si los edificios pueden ser usados o no lo mas pronto posible. Despues de un 
sismo moderado o severo, es importante identificar las edificaciones que han sufrido serios danos 
y que pueden representar peligro para sus ocupantes o la comunidad en general, asi como 
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aquellas que pueden ser utilizadas como refugios temporales para personas que han sido 
evacuadas. Algunos paises han desarrollado manuales y procedimientos para evaluar el dano en 
edificios, algunos son: Japon, Macedonia (antes Yugoslavia), Italia, Grecia, Estados Unidos, 
Mexico y Colombia, entre otros. Las evaluaciones de dano no son solo utiles para la toma de 
decisiones sobre las estructuras, sino que tambien sirven para hacer mejoras a los codigos sismo - 
resistentes, identificando los tipos de falla en los sistemas estructurales y las razones por las que 
se presentaron. Ademas, con la informacion a cerca del dano en los edificios es posible conformar 
funciones de vulnerabilidad empirica. 



2.1 Problemas con la evaluacion del dano 



Cuando el dano en el area afectada es generalizado, los expertos locales en ingenieria 
estructural son insuficientes para hacer la totalidad de las evaluaciones, asi que profesionales 
voluntaries pero inexpertos deben hacer parte del trabajo. A pesar de la existencia de manuales 
de evaluacion con completas descripciones de los diferentes niveles de dano, la informacion que se 
maneja es altamente subjetiva y depende de la perception del evaluador en cada caso, por lo cual 
los evaluadores inexpertos tienden a calificar erroneamente los danos. Pueden cometerse errores 
como permitir la ocupacion de edificios que debian ser evacuados, o dar orden de demolition a 
estructuras que eran reparables. Estos errores llevan a serias complicaciones, especialmente en el 
caso de edificios esenciales, como los hospitales. Los niveles de dano son clasificados en los 
manuales con calificaciones linguisticas como leve, moderado, severo o fuerte, estos conceptos 
pueden tener diferentes significados dependiendo de la persona que los maneje, y no existen 
limites claramente definidos entre ellos. 

3 MODELO PARA LA EVALUACION POST-SISMICA DEL DANO 



El modelo propuesto utiliza la teoria de conjuntos difusos motivado en el caracter 
subjetivo e incompleto de la informacion. La evaluacion post-sismica del dano utiliza 
calificaciones linguisticas cualitativas que pueden ser manejadas apropiadamente con la teoria de 
los conjuntos difusos. Por otro lado, se utiliza una red neuronal artificial para calibrar el sistema 
utilizando el criterio de especialistas. Esto posibilita el uso de estas herramientas para la 
evaluacion del dano por inexpertos. 

Para el desarrollo del modelo se tuvieron en cuenta varios manuales de evaluacion 
existentes, tales como los desarrollados en Mexico (Rodriguez, 1995), (SMIS, 1998), Japon 
(CENAPRED, 1996), Estados Unidos (ATC, 1989), Italia (Goretti, 2001), Macedonia (antes 
Yugoslavia) (IZIIS, 1985), la metodologia usada despues del sismo de 1999 en Colombia (Campos, 
1999) y la metodologia desarrollada para la ciudad de Bogota tambien en Colombia (AIS, 2002). 
Ademas miembros de la Asociacion Colombiana de Ingenieria Sismica apoyaron tecnicamente 
este trabajo. El modelo ha sido implementado en Visual Basic 6.0. El programa es de uso 
amigable al usuario, y ofrece ayudas al evaluador, tales como fotografias y descripciones de los 
tipos de danos. 

Este modelo ha sido implementado para ser una herramienta de uso oficial para la 
Alcaldia de la de ciudad de Manizales, Colombia y se realizo en el marco del proyecto “ Evaluacion 
de edificaciones despues de un sismd’, en el cual se desarrollo un “Manual de campo para la 
inspection de edificios despues de un sismo” acompanado de un formulario de evaluacion. 



3.1 Estructura de la red neuronal artificial 



La red neuronal tiene tres capas. Las neuronas en la capa de entrada son agrupadas en 4 
tipos, elementos estructurales (EE), no estructurales (NE), condiciones del suelo (CS), y 
condiciones preexistentes (CP). Cada una contribuye con informacion a una neurona intermedia, 
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la cual corresponde al grupo de variables al que pertenece. El numero de neuronas de entrada o 
de variables en el modelo, no es constante. Este numero depende del sistema estructural que se 
evalua y de la importancia de los diferentes grupos de variables seleccionados para la evaluation. 
En la tabla 1 se ve como las variables de entrada cambian de acuerdo con el sistema estructural. 
Se asigna una calificacion dependiendo del dano observado, usando cinco posibles niveles de dano 
que son tratados como conjuntos difusos. Para elementos estructurales y no estructurales los 
niveles corresponden a las siguientes calificaciones linguisticas^ ninguno o muy leve (N), leve (L), 
moderado (M), fuerte (F) y severo (S) o destruction. En la figura 1 puede observarse algunos 
niveles de dano que pueden presentarse en diferentes elementos estructurales. La figura 2 ilustra 
las funciones de pertenencia para estas calificaciones linguisticas. Los conjuntos difusos son 
basados en indices de dano seleccionados, de los cuales se hablara mas adelante. 



Tabla 1 . Elementos estructurales de acuerdo con el sistema estructural 



Sistema Estructural 


Elementos Estructurales 


Porticos en concreto reforzado (con) muros de 
cortante 


Columnas/ muros, vigas, nudos y 
entrepisos 


Porticos de acero o madera 


Columnas, vigas, conexiones y entrepisos 


Mamposteria No reforzada / Reforzada/ 
Confinada 


Muros de soporte y entrepisos 


Muros de bahareque o tapia 


Muros de soporte y entrepisos 



El dano en los elementos no estructurales no pone en peligro la estabilidad del edificio, 
pero puede representar un peligro para la seguridad de los ocupantes. Los elementos no 
estructurales estan divididos en dos grupos- elementos comunes y elementos opcionales. La tabla 
2 presenta esta clasificacion. Algunos ejemplos de los danos que pueden observarse en elementos 
no estructurales pueden observarse en la figura 3. 



Tabla 2. Elementos no Estructurales 



Elementos Comunes 


Muros divisorios 


Muros de fachada 


Escaleras 


Elementos 

Opcionales 


Cielo rasos y 
luminarias 


Instalaciones 


Cubierta 


Tanques elevados 
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c) d) 

Figura 1. Danos en elementos estructurales, a) Dano severo en un nudo de concreto reforzado, b) Dano 
moderado en una viga de concreto reforzado, c) Dano fuerte en un miiro de mamposteria d) Dano fuerte en un 

muro de tapia. 




Muy Leve — — Leve - - - Moderado — - Fuerte — - Severo 



Figura 2. Funciones de pertenencia para las calificaciones linguisticas. 
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c) 



d) 



Figura 3. Danos en elementos no estructurales, a) Dano severo en un muro divisorio, b) Dano severo en una 
escalera, c) Dano moderado en un muro de fachada d) Dano severo en un tanque. 



Las variables de condiciones del suelo y preexistentes son valorados utilizando 
calificaciones linguisticas, de acuerdo con el estado observado en el momento de la evaluation. 
Las calificaciones linguisticas son: muy buena (MB), buena (B), regular (R), mala (M), y muy 
mala (MM). En las condiciones del suelo se incluyen efectos que pueden inducir una falla en la 
estabilidad de la edification que se estudia, tales como los asentamientos, subsidencias y 
licuacion del suelo y la ocurrencia de deslizamientos o fallas en el terreno! ejemplos de estas se 
pueden observar en la Figura 4. Las condiciones preexistentes se refieren a la calidad de los 
materiales de construction, irregularidades en planta y altura del edificio, y la configuration 
estructural (ejemplos en la figura 5). 

En la capa oculta o intermedia, se obtiene un indice para cada grupo de variables por 
medio de la agregacion y desfusificacion. Teniendo en cuenta los cuatro indices, es posible definir 
en la capa de salida el dano en el edificio utilizando una base de reglas difusas con las 
valoraciones del dano estructural y no estructural. La habitabilidad del edificio es obtenida 
teniendo en cuenta las condiciones del suelo, y finalmente utilizando las condiciones 
preexistentes se define el nivel de reparation necesario para el edificio. De este modo, este 
sistema puede dar recomendaciones sobre la habitabilidad y reparabilidad del edificio. Asi con 
este sistema puede llegarse a resultados como “habitable con reparaciones menores”, “restringido: 
utilizable despues de reparation”, o “peligroso : posible demolition”. 
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a) 



b) 



Figura 4. Condiciones del suelo, a) Asentamientos y licuacion del suelo, b) Deslizamientos de tierra y 

grietas en superficie. 




c) d) 

Figura 5. Condiciones preexistentes, a) Mala calidad de construction, b) Irregularidades en altura, piso 
blando; c) Irregularidad en planta, d) Mala configuration estructural, no hay continuidad vertical en los 

elementos. 



3. 1. 1 Capa de entrada 



Los conjuntos difusos para cada variable i (columnas, muros o vigas), en la capa de 
entrada, son obtenidos de las calificaciones linguisticas del dano dadas por los evaluadores Di en 
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cada nivel j y su extension wj. La extension del dano (porcentaje de elementos de un tipo en cada 
nivel de dano) varia de 0 a 100 y es normalizado, Ecuacion (l). La calificacion acumulada del 
dano en cada variable Di es obtenida con la union de los conjuntos difusos escalados, teniendo en 
cuenta las funciones de pertenencia del dano Poj(Dj) y su extension o peso (asignada por el 
evaluador). La union entre conjuntos difusos esta definida por la maxima dependencia, esto es 
mas claro en las ecuaciones (2) y (3). Por medio de la desfusificacion, usando el metodo del 
centroide de area (COA), se obtiene un indice de calificacion C, para cada variable de cada grupo 
de neuronas. En la Ecuacion (4) describe el metodo de desfusificacion utilizado. En la Figura 6 se 
ilustra el ingreso de datos al sistema. 



(l) 

N 

D i = ( D N (jD L ( ~ ,D M ( - jD F^ ,D s) ( 2 ) 

P A (22) = max(w Ar j . *|X Djr {D Ni w Si (3) 

C,. = |max (w v , *p Djv {D Ni ),...,w st *p Dj (d w ) L W e (4) 




Figura 6. Ingreso de la extension del dano de cada nivel para un elemento estructural. 



Cada variable tiene predefinidas las funciones de pertenencia basicas de los conjuntos 
difusos correspondientes a los cinco niveles de dano posibles. Las calificaciones linguisticas 
cambian en cada caso como se menciono antes. 



3.1.2 Capa intermedia 



En esta capa se tienen cuatro neuronas, que corresponden a cada uno de los grupos de 
variables^ elementos estructurales, elementos no estructurales, condiciones del suelo y 
condiciones preexistentes. En la figura 7 se presenta un esquema general del proceso de 
evaluation. 
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En este modelo de red neuronal las neuronas de entrada no tienen influencia sobre todas 
las neuronas intermedias, solo lo hacen sobre la neurona que corresponde al grupo de variables al 
que pertenecen. Los parametros de entrada de estas 4 neuronas son las calificaciones C, 
obtenidas para cada variable de cada grupo de neuronas y su peso W, o el grado de importancia 
dados por el evaluador de acuerdo con su criterio. Estos pesos son normalizados y calibrados por 
medio de una funcion de aprendizaje de la que se hablara mas adelante. Para el proceso de 
entrenamiento, los valores iniciales han sido definidos teniendo en cuenta el criterio de expertos. 
Con estas calificaciones y pesos de cada variable i se obtiene un indice global para cada grupo k, 
por medio de la union y desfusificacion de los conjuntos difusos escalados. Esto se hace con las 
funciones de pertenencia Hcu(Cid) y los pesos Wki como se indica en las ecuaciones (5) y (6). 



\^cee (c) — max (fV EEJ * \i- C]!Ei {C EEI ),■■■, W FFj * |i Cw . (C EEj )) (5) 

I ee = [max {fV EE ! * ( (C EE1 ).•••. W EFi * li C j (C EEi ))\ centroide (&) 



Las variables de condiciones del suelo y preexistentes pueden ser tomadas o no en cuenta 
dentro de la evaluacion, dependiendo de su importancia en los resultados, esto hace que en 
algunos casos la habitabilidad y reparabilidad puedan ser obtenidas unicamente a partir del dano 
estructural y no estructural. 




Figura 7. Estructura de la red neuronal propuesta 
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3. 1. 3 Capa de salida 



En esta capa se obtienen indices globales para los elementos estructurales, no 
estructurales, condiciones del suelo y preexistentes, corresponden a una calificacion linguistica 
final en cada caso. El nivel de dano es obtenido de acuerdo con la proximidad del valor del indice 
obtenido a la funcion de dano de referenda. Los niveles para los elementos estructurales y no 
estructurales son definidos de igual forma que antes, como- ninguno, leve, moderado, fuerte y 
severo. Tambien para las condiciones preexistentes y del suelo se definen asi : muy buenas, 
buenas, regulares, malas y muy malas. En esta capa tiene lugar el proceso de entrenamiento, los 
indices que identifican cada nivel cualitativo (centro de cluster) cambian de acuerdo con los 
indices calculados en cada evaluation que se hace con la red, teniendo en cuenta una tasa de 
aprendizaje. Esta ultima determina la velocidad de entrenamiento de la red. El procedimiento de 
aprendizaje es descrito en el siguiente apartado. Despues de hallar las calificaciones finales 
utilizando bases de reglas difusas se determinan el dano global en el edificio, su habitabilidad y 
reparabilidad. En secciones posteriores se explicara este procedimiento. 

3.2 Aprendizaje de la red neuronal 



La red neuronal es calibrada en la capa de salida, cuando se definen las funciones de dano 
en relation la matriz de indices de dano. Para dar inicio a este proceso se le dio un punto de 
partida a la red, se definio un valor initial para cada nivel de dano. Para la determination de 
estos indices iniciales han sido considerados los propuestos por el AT013 (ATC, 1985), Park, Ang 
y Wen (Park et al., 1984), las curvas de fragilidad usadas por Hazus-99 (FEMA, 1999) y los 
indices utilizados por Sanchez-Silva (Sanchez-Silva et al. 2001). Los valores de los indices 
iniciales corresponden a los centroides de area de cada funcion de pertenencia a cada nivel de 
dano. La tabla 3 presenta una comparacion entre los indices aqui propuestos, los de Park, Ang y 
Wen, y los de Sanchez-Silva y Garcia. 



Tabla 3. Tabla comparativa de los indices de dano 



Nivel de dano 


Park, Ang y Wen 


Sanchez Garcia 


Propuestos 


Muy leve 


<0.1 

0.07 


0.10 


0.07 


Leve 


0.10-0.25 

0.175 


0.20 


0.17 


Moderado 


0.25-0.40 

0.325 


0.35 


0.33 


Severo 


0.40-0.80 

0.6 


0.60 


0.55 


Destruction 


O 

GO CO 

o o 
A 


0.90 


0.76 



Los indices iniciales para la matriz de dano estan basados en los indices de Park, Ang y 
Wen, dado que estos han sido calibrados con diferentes estudios y danos de terremotos reales y 
son los que cuentan con mayor aceptacion. Sus autores consideran la existencia de colapso desde 
0.8, mientras que Stone y Taylor (Stone et al. 1993) lo proponen a partir de 0.77. Considerando 
esto se tomo 0.76 como valor para el indice del nivel de destruction o colapso. En la election de los 
indices de los niveles severo y moderado se ha intentado ser conservatives, dado la intensa 
discusion sobre estos para reducir su valor. 



La calibration se hace por nivel de dano y solo para los indices correspondientes a los grupos 
de variables considerados en cada evaluation que calibra el sistema. En la capa de salida el 
aprendizaje utiliza el algoritmo de Kohonen, como se ve en la ecuacion (7). 



196 



Martha L. Carreno, Omar D. Cardona, y Alex H. Barbat 



(7) 

donde I kj es el valor del indice de un grupo de variables k que presenta un nivel j, que se recalcula 
considerando una tasa de aprendizaje a y la diferencia entre el indice resultante de la presente 
evaluacion y el indice previo en cada nivel de dano j. En la ecuacion (8) se define la tasa de 
aprendizaje de la ecuacion (7). 



a (t)=0.1*Exp(- 0.1*t ) 



( 8 ) 



donde t es el numero de veces que el indice o peso que se calibra ha sido usado. Para el 
entrenamiento de la red se utilizaron las evaluaciones hechas despues del terremoto del Quindio 
en Colombia en 1999. En el entrenamiento no se tuvieron en cuenta sistemas estructurales como 
porticos en madera y acero, esto debido a la ausencia de estos sistemas estructurales en la zona 
afectada por el sismo. Se utilizaron pocos edificios con porticos en concrete reforzado con muros 
de cortante, esto debido al reducido numero de evaluaciones de estructuras de este tipo. 



3.3 Bases de reglas fuzzy 



La habitabilidad y reparabilidad del edificio que se evalua son estimadas basandose en el 
nivel de dano de los elementos estructurales y no estructurales, el estado del suelo y las 
condiciones preexistentes. La Figura 8 ilustra el uso de reglas difusas para estas estimaciones. El 
nivel de dano global del edificio es deducido a partir de los resultados de dano estructural y no 
estructural! y tiene cinco posibles calificaciones : ninguno, leve, moderado, fuerte y severo. El 
estado global del edificio y la habitabilidad son estimados teniendo en cuenta el dano global del 
edificio y las condiciones del suelo en dos diferentes bases de reglas fuzzy. La habitabilidad puede 
ser clasificada en cuatro niveles- habitable, uso restringido, prohibida la entrada y peligroso. La 
reparabilidad depende de otra base de reglas que contempla las condiciones preexistentes de la 
edification, aqui se contribuye a definir la factibilidad tecnica y economica de una reparacion. La 
reparabilidad tiene cuatro posibles niveles : ninguna o reparaciones menores, reparacion, 
refuerzo, y posible demolition, estos niveles son identificados usualmente con colores en los 
manuales de evaluacion, asi : verde, amarillo, naranja y rojo. Para cada edificio que se evalua con 
este sistema es posible ayudar en dos decisiones criticas que son basicas en la fase de respuesta a 
la emergencia despues de ocurrido un sismo fuerte. En la figura 8 se ilustra la presentation de los 
resultados dada por el programa. 



4 CONCLUSIONES 



Despues de una completa revision de los diferentes manuales para evaluacion post-sismica de 
dano, se propone un sistema experto innovador que apoye el proceso de evaluacion. Para el 
desarrollo de esta herramienta se han tenido en cuenta las ventajas y desventajas de cada uno 
de los metodos. 

El sistema experto fue desarrollado utilizando tecnicas de inteligencia artificial como las redes 
neuronales artificiales y la teoria de conjuntos difusos. Este tipo de herramientas es apropiado 
en esta practica, debido al caracter subjetivo de la evaluacion de los danos y que la 
information suele ser incompleta. 

Evaluaciones hechas por ingenieros expertos despues del terremoto del Quindio, en Colombia, 
en 1999, han sido utilizadas para el entrenamiento del sistema. 
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Indice Estructural 
Dano Estructural 



Indice No Estructural 
Dano No Estructural 



N. No 

^Estructural 

Estructural^. 


Ningun 

0 


Leve 


Moderad 

o 


Fuerte 


Severn 


Ninguno 


N 


L 


L 


M 


M 


Leve 


L 


L 


M 


M 


M 


Moderado 


M 


M 


M 


F 


H 


Fuerte 


F 


F 


F 


S 


S 


Severn 


S 


S 


S 


S 


S 



Base de reglas para Dano en el Edificio 
Dano Estructural vs. Dano No Estructural 






Indice de Condiciones Preexistentes 
Condiciones Preexistentes 



Indice de Estado Global 
Estado del Edificio 



T. 



y 



Preexistentes 

Estado^^^ 

Edificio 


Muy 

buena 

s 


Buena 

s 


Regulare 

s 


Mala 

s 


Muy 

Mala 

s 


Ninguno 


M 


M 


M 


R 


S 


Leve 


M 


M 


R 


S 


S 


Moderado 


R 


R 


S 


S 


S 


Fuerte 


S 


S 


S 


D 


D 


Severn 


D 


D 


D 


D 


D 



Base de reglas para Reparabilidad del Edificio 
Estado del edificio vs. Condiciones Preexistentes 











Nv Suelo 
Edificio 


Muy 

Buena 

s 


Buena 

s 


Regulare 

s 


Mala 


Muy 

Mala 

s 




Ninguno 


N 


L 


M 


F 


S 


Leve 


L 


L 


M 


F 


S 




Moderado 


M 


M 


F 


S 


S 
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Base de reglas para Estado del Edificio 
Dano en el Edificio vs. Condiciones del Suelo 
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Base de reglas para Habitabilidad del Edificio 
Dano en el Edificio vs. Condiciones del Suelo 





CALIF1CACIONES: 




Dano Estructural y No Estructural: 


Habitabilidad: 


N: Ninguno 


H: Habitable 


L: Leve 


R: Restringido 


M: Moderado 


P: Prohibido 


F: Fuerte 


D: Peligroso 


S: Severo 




Condiciones del Suelo y Preexistentes: 


Reparabilidad: 


MB: Muybuenas 


M: Menor o ninguna 


B: Buenas 


R: Reparation 


R: Regulares 


S: Refiierzo 


M: Malas 


D: Posible demolition 


MM: Muy Malas 






Figura 8. Procedimiento para determinar la Habitabilidad y Reparabilidad 



Se confirmo la necesidad de tener un metodo preestablecido para la evaluation del dano 
en edificios en caso de un sismo fuerte en ciudades ubicadas en zonas de sismicidad 
considerable. 

El uso de herramientas de inteligencia artificial en la ingenieria civil ha tenido muy poca 
difusion, es recomendable promover su uso para dar soluciones adecuadas y rapidas a 
problemas en este campo del conocimiento. 
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a) 




b) 



RESULTADOS 

Resultados J Estado } Habitabilidad } Reparabilidad~ 



Habitabilidad del edificio No Habitable 


Habitabilidad 


La ocupacion de la edification debe prohibirse porque 
presenta peligro para la seguridad de los ocupantes debido al 
muy grave estado del edificio y las graves condiciones del 
suelo. 


Medidas de seguridad 




Es necesario remover o anclar elementos que estan en peligro 
de caerse y tomar medidas preventivas debido a los graves 
problemas geotecnicos que se presentan. 


« Devolve! | Continue 



Guardar Informacion Aceptar 




c) 



d) 



Figura 9. Presentation de los resultados de una evaluation de danos realizada con el sistema “EDE- 

Manizales” 
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RESUMEN 



Se presenta un estudio de peligrosidad sismica para la ciudad de Salta mediante el cual se 
determina la aceleracion maxima del suelo en Roca para cuatro eventos sismicos denominados : 
frecuente, ocasional, raro y muy raro que tienen periodos de retorno de 43, 72, 475 y 970 anos. En 
base a estos resultados y considerando como base la forma del espectro elastico de la Norma 
INPRES CIRSOC 103 vigente al 2004 se determinan cuatro espectros para los sismos indicados, 
con los cuales se va a verificar el desempeno sismico de un edificio de Hormigon Armado 
construido en la ciudad de Salta en el 2000. 

Se determina la curva de capacidad sismica resistente del edificio, en forma espacial, 
empleando el programa CEINCI3 y los resultados se comparan con los obtenidos por el programa 
ETABS. Posteriormente se encuentra el espectro de capacidad del edificio y se hallan los 
desplazamientos laterales maximos ante cada uno de los sismos anotados en el parrafo anterior 
empleando el metodo del espectro de capacidad. Con los desplazamientos laterales se halla en la 
curva de capacidad sismica el desempeno que va a tener la estructura, todo esto en forma grafica. 

Finalmente se indica el desempeno que tendria la estructura en analisis modificando la 
calidad del hormigon y el detallado del refuerzo transversal tanto en vigas como en columnas. 



ABSTRACT 



A seismic hazard study for the city of Salta deternining the highest soil acceleration in 
rock is presented. Four seismic events are determined: frequent, occasional, rare and very rare 
with recurrence intervals of 43, 72, 475 and 970 years. Using these results and considering as 
basic the elastic spectyrum form of the INPRES CIRSOC 103 Norm in force at 2004, four 
spectrum for the indicated seisms are determined, with which the seismic performance of a 
reinforced concrete structure built in Salta in 2000 is verificated. 

The seismic resistant capacity of the building is determined, in 3'd form, using the 
CEINCI3 program and the results are compared with those obtained with the ETABS program. 
After that the capacity spectrum of the building is founded and the maximum lateral 
displacement for each of the seisms anotated in the previous paragraph are also solved using the 
capacity spectrum method. With these lateral displacements the performance in the seismic 
capacity curve is calculated, all of these in graphic form. 
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Finally the structural performance of the building modifying the concrete quality and the 
detailment of the transversal reinforcement in beams and in columns is indicated. 

1. INTRODUCTION 

Salta es una de las provincias con mayor superficie de la Argentina, tiene 155488 Km 2 y 
tiene una poblacion aproximada de un millon de habitantes de tal manera que su densidad 
poblacional es muy baja. La provincia se encuentra al noroeste y limita al norte con la provincia 
de Jujuy y la Republica de Bolivia, al este con las provincias del Chaco y Formosa y la republica 
de Paraguay, al sur con las provincias de Catamarca, Tucuman y Santiago del Estero y al oeste 
con la Republica de Chile. La capital de la Provincia es la ciudad de Salta que tiene una gran 
riqueza colonial. 

Historicamente la Provincia de Salta se ha visto afectada por los movimientos teluricos 
como se puede apreciar en la figura 1 en que se indican los principales eventos registrados desde 
la epoca de la colonia hasta inicios del siglo XX, Quispe et al (2003). Basta recordar el sismo del 
13 de septiembre de 1692 de magnitud 7.3 que destruyo la ciudad de Esteco una de las 
poblaciones mas importantes de la Provincia. 
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MGURA 1 : Mapa de la provincia de Salta- Sismos histdricos 



Por otra parte el gran avance que ha tenido el Analisis Sismico por Desempeno, a nivel 
mundial, Seaoc (1995), Chopra y Goel (1999), Priestley (2003), ha originado el desarrollo de 
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importantes programas de computation como el CEINCI2 que sirve para el analisis en dos 
dimensiones y el CEINCI3 para el analisis en 3 dimensiones, Aguiar (2002 y 2003). 

A1 contar con los programas mencionados u otros, es factible predecir en forma virtual el 
probable desempeno, en terminos estructurales y economicos, que tendra un edificio ante 
diferentes sismos. A1 conocer en un ordenador la forma como se va a comportar se pueden tomar 
las acciones necesarias para mitigar los danos que puedan registrarse en caso de un sismo 
determinado y de esa manera se estaria precautelando las vidas y los bienes materiales. 



2. SISMICIDAD DE LA PROVINCIA DE SALTA 



En la figura 2 se aprecia los epicentros de los sismos registrados en la Provincia de Salta 
de acuerdo al catalogo cedido por el Instituto National de Prevention Sismica INPRES (2004) y 
que han sido complementados con los registros historicos por Quispe et al (2003). Notese que 
existe un nido sismico en la zona nor oeste. 




En el siglo XX el sismo mas destructive que afecto a la Provincia de Salta fue el registrado 
el 25 de agosto de 1948 que tuvo una magnitud de 7 y una intensidad de IX y a mas de Salta fue 
afectada la provincia de Jujuy. 

3. ZONAS SISMOGENETICAS Y MAGNITUD MAXIMA 



Para determinar la Peligrosidad sismica de la ciudad de Salta se considero el area 
geografica comprendida entre los paralelos -22 y -26.3 grados de latitud sur y los meridianos -62 y 
-68.2 grados de longitud oeste, que es la que se ha venido presentando en las figuras 1 y 2. En 
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este rectangulo se analizo la ubicacion de los epicentros, las fallas geologicas, la informacion 
tectonica, se hicieron varios cortes y se ubico los hipocentros de los sismos. Con toda la 
informacion historica disponible y la instrumental se han definido cinco zonas sismogeneticas las 
mismas que se indican en la figura 3 y brevemente se describen a continuation. Los numeros que 
constan en los recuadros de dicha figura corresponden al numero de sismos registrados. 



3.1 Zona Sismogenetica 1 



Ubicada en el NE de la Provincia de Salta, con una sismicidad superficial y profunda. Del 
analisis de la tectonica de la region, coincide con el frente de deformation de las Sierras 
Subandinas, que comprende una serie de fallas activas, de tipo inversa, con rumbo aproximado 
norte'sur, en escamas. 

En cuanto la magnitud maxima que puede esperarse, solo puede darse una aproximacion 
considerando las longitudes totales de fallas extraidas de los mapas geologicos disponibles, 
aplicando abacos segun Slemmons (1982). Se obtiene para una falla de 200 Km. un sismo de 
magnitud Mmax=6.8 aproximadamente. Haciendo el analisis de los datos segun la ley de 
Gutemberg-Richter, se determina que para sismos superficiales cabe esperar una magnitud 
maxima Ms=7.1, mientras que si se consideran todos los eventos (superficiales y profundos) se 
obtiene una Ms=6.9. 

En el catalogo, la maxima magnitud correspondiente a la sismicidad superficial es de 
Ms=6.4, por lo que se considera apropiado considerar una magnitud maxima de 7, una distancia 
promedio de 270 Km. y una profundidad de 30 Km. 




3.2 Zona Sismogenetica 2 



Ubicada al SE de la Ciudad de Salta. Se caracteriza por una sismicidad intensa 
superficial, que ha afectado a la Ciudad de Salta grandemente. Los mayores sismos historicos han 
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sido originados en esta region sismogenetica. Del analisis de la tectonica de placas, la region esta 
surcada por varias fallas activas que son continuation del frente de ruptura norte, pero se 
considera una fuente separada ya que la caracteristica eminentemente superficial de sus eventos 
intensos hace que sean determinantes en la definition de la peligrosidad sismica de la zona en 
estudio. 



La magnitud maxima esperada (por longitud total de fallas) es de 6.7. En el catalogo se le 
han asignado a dos sismos historicos valores de 7 y segun ley de Gutemberg-Richter, se obtiene 
una magnitud maxima de Ms 7.44 para sismos superficiales que varia muy poco si se consideran 
tambien los eventos profundos, por lo que a los fines de este estudio y con los datos disponibles, se 
considera apropiado asignarle el valor de 7.5., una distancia media de 100km y una profundidad 
entre 20 y 40 Km. 



3.3 Zona Sismogenetica 3 



Localizada al NO de la ciudad, a una distancia promedio de 220 Km. Se caracteriza por 
una sismicidad profunda, cuyos efectos en general no son percibidos en la Ciudad de Salta. Los 
movimientos originados en esta region corresponden a la zona de Benioff por lo que se deben a la 
subduccion de la placa del Pacifico en la Sudamericana, de alii el caracter profundo de los 
mismos. Considerando la ley de Gutemberg-Richter, se tiene una magnitud maxima de Ms=6.6 y 
de caracter eminentemente profunda, 



3.4 Zona Sismogenetica 4 



Se halla al SO de la ciudad. Corresponde a una sismicidad superficial, debida al 
movimiento de fallas locales, que muestran en general una sismicidad moderada, y que afecta 
principalmente a las localidades ubicadas en los Valles Calchaquies. Esta ubicada a una 
distancia media a la Ciudad de Salta de 160 Km. con una profundidad de 45 Km. y una magnitud 
maxima probable de 6.7. El valor maximo reportado en el catalogo es de 6.5 para la sismicidad 
superficial de la zona. Segun la ley de Gutemberg-Richter, se debe esperar un sismo maximo de 
magnitud Ms=6.3 considerando los datos de sismos superficiales y de 6.9 considerando todos ellos. 
Por ello y dado la escasez de datos para hacer mayores inferencias, se le asigna un valor de 6.9 . 



3.5 Zona Sismogenetica 5 



Se ubican en esta zona escasos eventos cercanos a la Ciudad de Salta, son eventos 
superficiales con una profundidad focal media de 49.21 Km. y una magnitud maxima reportada 
en el catalogo de 5.8. 

Se considera que las zonas 2 y 5 son determinantes en la estimation de la amenaza 
sismica para la Ciudad de Salta, debido a su actividad (zona 2) y cercania (zona 5). En menor 
medida se ubican la zona 1, zona 4 y por ultimo la zona 3. 



4. SISMICIDAD DE CADA ZONA SISMOGENETICA 



La information sismica instrumental en la mayor parte de las regiones del mundo 
empieza a tenerse a partir de 1960 en que se comenzaron a instalar las redes sismicas. Antes de 
esa fecha en los catalogos sismicos no hay information de los sismos de baja magnitud, 
unicamente existe information de los sismos de gran magnitud. Por este motivo en todo estudio de 
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peligrosidad sismica se realiza en primer lugar un estudio que se denomina de completitud de la 
informacion sismica, que consiste en determinar el tiempo en el cual la tasa anual de ocurrencia 
de los sismos permanece constante para ello se clasifica la informacion sismica registrada en cada 
area fuente por magnitudes y por intervalos de tiempo en cada una de las zonas sismogeneticas. 

Por motivo de espacio no se indica el estudio realizado de completitud de la informacion 
sismica. Unicamente se manifiesta que en la zona de estudio y para todas las zonas 
sismogeneticas se asegura completitud en el periodo comprendido entre 1963 y 2003 es decir 41 
anos. Para este periodo de tiempo los resultados del estudio estadistico son los indicados en la 
Tabla 1. 

Los resultados que muestra la tabla 1 corresponden a la sismicidad superficial, entendida 
como tal la que se registra en los primeros 70 Km. de la corteza. La columna 1 designa a las cinco 
zonas sismogeneticas consideradas. La columna 2 muestra las medias de los valores de magnitud 
del catalogo sismico, y la columna 3 las desviaciones estandar. H denota la profundidad media de 
los focos, a y b son los parametros de la relacion de Gutemberg-Richter, determinadas por 
regresion. Las dos ultimas columnas muestran las magnitudes esperadas para periodos de 
retorno de 43, 72, 475 y 970 anos, respectivamente, determinados considerando la ocurrencia de 
sismos como un proceso de Poisson (Msti) y de valores de extremos Gumbel tipo I (Mst2). 



Ta bla 1 Resultados para un periodo de tiempo de 41 anos, y sismicidad superficial 



Zona 


Media 

Ms 


Std 


H 

Km. 


A 


b 


Msti 


MsT2 


1 


4.59 


0.38 


34.00 


4.00 


-1.07 


5.25 


5.53 














5.46 


5.68 














6.22 


6.24 














6.51 


6.46 


2 


4.70 


0.45 


32.90 


3.15 


-0.78 


6.11 


5.81 














6.40 


5.99 














7.44 


6.66 














7.84 


6.91 


3 


4.30 


0.36 


31.10 


5.78 


-1.53 


4.85 


5.19 














5.00 


5.34 














5.54 


5.87 














5.74 


6.07 


4 


4.19 


0.54 


60.00 


6.31 


-1.62 


4.92 


5.53 














5.05 


5.74 














5.56 


6.54 














5.75 


6.84 


5 


4.36 


0.62 


51.00 


3.98 


-1.11 


5.07 


5.89 














5.27 


6.14 














6.01 


7.06 














6.29 


7.41 



En la tabla 2 se muestran los periodos de retorno para sismos de magnitud 6, 6.5, 7 y 7.5 
correspondientes a cada una de las zonas sismogeneticas consideradas, para sismicidad 
superficial y un periodo de datos de 41 anos y considerando un proceso de Poisson para el 
tratamiento estadistico. 



Tabla 2 Periodos de retorno para sismos intensos (Poisson) 



Zona 


Ms=6 


Ms=6.5 


Ms=7 


Ms=7.5 


1 


277 


953 


3.280 


11.300 


2 


35 


87 


214 


526 


3 


2.440 


14.300 


82.000 


477.000 


4 


2.440 


15.700 


100.000 


650.000 


5 


463 


1.660 


5.930 


21.200 
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5. RESULTADOS DEL ESTUDIO DE PELIGROSIDAD SISMICA 



Se aplicaron varias formulas de atenuacion de la aceleracion maxima del suelo en funcion 
de la magnitud y distancia epicentral, y se determinaron los valores probables de aceleraciones 
maximas en roca para la Ciudad de Salta. Tales formulas no son locales, pues no se tiene aun 
una expresion basada en registros locales por lo que se utilizo las que fueron determinadas para 
la zona oeste-central de pais. 

Los valores de magnitud utilizadas fueron en cada caso el valor mayor de los obtenidos 
utilizando las dos distribuciones mencionadas anteriormente. 

Los resultados obtenidos, se comparan con los determinados en base a consideraciones 
realizadas para la determination de la zonificacion sismica por parte de los sismologos del 
Instituto National de Prevention Sismica y que dieron origen a los espectros de diseno vigentes en 
la Norma Sismorresistente actual. A tal fin se reprodujeron los resultados publicados y luego se 
hicieron los calculos para los periodos de retorno considerados por Vision 2000. 

Teniendo en cuenta la envolvente de las zonas sismogeneticas consideradas, tal como fue 
realizado por el proceso de zonificacion sismica realizada por el INPRES, se tienen los valores de 
aceleraciones maximas para roca en la Ciudad de Salta, como porcentaje de la aceleracion de la 
gravedad indicados en la Tabla 3. La aceleracion maxima en roca, en que se anclan los espectros 
de respuesta elasticos, para los tres tipos de suelos considerados se considera igual a 0.25 g (la 
norma INPRES CIRSOC considera aproximadamente un sismo raro). 



Tabla 3 Aceleraciones maximas en roca 



Tiempo de retorno 
en anos 


Este estudio 


INPRES 


43 


0.12 


0.065 


72 


0.14 


0.084 


475 


0.29 


0.222 


970 


0.37 


0.320 



Se vuelve a resaltar que siempre se eligieron los valores mayores, tanto de magnitud, 
como de las aceleraciones obtenidas en base a los distintos modelos de atenuacion considerados. 

6. ESPECTROS DE ANALISIS 



Se determinan los espectros de diseno para los sismos frecuente, ocasional, raro y muy 
raro en base al espectro de diseno de la norma INPRES'CIRSOC 103 en vigencia, mediante el 
simple metodo de multiplicar este ultimo por parametros iguales a la relation entre los valores de 
aceleraciones maximas en roca obtenidas. Resulta asi que el espectro del sismo frecuente se 
obtiene multiplicando el espectro de diseno existente por 0.48! el ocasional, por 0.56, el raro por 
1.16 y el muy raro por 1.48. 

La estructura que se va a analizar se halla ubicada en un suelo tipo SI, razon por la cual 
unicamente se indican las formas espectrales para ese tipo de suelo en la figura 4. El espectro 
notado con 1 corresponde al sismo frecuente, el 2 al ocasional, el 4 al sismo raro y el 5 al sismo 
muy raro. El espectro 3 es el que esta inmerso en la Norma INPRES-CIRSOC 103 que no se va a 
utilizar en el presente estudio. 



208 



Roberto Aguiar Falconi, Lia Haarala Orosco y Roberto Caro Bassani 




Figura 4 Espectros de analisis considerados en el estudio para un perfil de suelo SI 

7. ESTRUCTURA DE ANALISIS 



Se va a analizar el desempeno que tendria un Edificio de 3 pisos situado en la ciudad de 
Salta que fue construido en el 2000! con un hormigon que tiene una resistencia maxima a la 
compresion de 17 MPa y un acero de 420 MPa. No se pretende alarmar a los usuarios del edificio 
sobre el desempeno que va a tener el mismo ante los sismos con los cuales se van a analizar por lo 
que se evitara dar mayores detalles de la estructura. 

A groso modo se indica que la estructura esta conformada por un sistema aporticado en 
base a vigas y columnas rectangulares. La losa de piso es maciza de 15 cm. de espesor. El refuerzo 
transversal de las vigas esta constituido por un estribo simple de 6 mm. de diametro espaciado 
cada 12 cm. en los apoyos y 16 cm. en el centro de luz! las columnas tienen estribos simples de 8 
mm. de diametro y con la misma separation. 

A1 sentido mas largo de la planta indicada en la figura 5 se denomina direction X y al 
sentido perpendicular, direction Y. En sentido X se tienen 4 porticos con luces de 7.50 m, para los 
vanos exteriores y 2.50 m. para el vano central. En sentido Y se tienen 6 porticos espaciados cada 
6.0 m. 



Una columna tipica tiene una section transversal de 40/45 cm. y una armadura 
longitudinal de 16 (j) de 12 mm. en la planta baja, en los pisos superiores esta armadura va 
disminuyendo. Las vigas son de 30/55 en la direction X con armadura longitudinal de 2 (f) de 12 
mm. en la parte superior y 3 0 de 12 mm. en la parte inferior! ademas de ello existen algunos 
adicionales de refuerzo en los extremos de 12 mm. de diametro. En el sentido transversal las 
vigas que cubren la luz de 7.50 m. son de 30/70 y la viga del centro de 30/55 con armadura un poco 
mayor a la de los porticos en sentido X. 

El hierro utilizado es Acindar que tiene en realidad un limite de fluencia de 480 MPa 
pero que se lo comercializa como de 420 MPa. El esfuerzo maximo a la rotura es de 770 MPa! la 
deformation del acero al initio de la zona de endurecimiento es 0.018 y el modulo en la zona de 
endurecimiento es 2116.788 MPa. Estos datos fueron utilizados para definir el modelo trilineal 
del acero con el cual se realizo el analisis por desempeno. Para el hormigon se trabajo con el 
modelo de hormigon confinado de Park y Priestley (1982). 
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Se considero una carga muerta de 1.11 T/m 2 y una carga viva de 0.25 T/m 2 para los pisos 
uno y dos. Para el tercer piso la carga muerta fue de 0.90 T/m 2 y la carga viva de 0.10 T/m 2 . Por 
otra parte cada planta tiene un area de 525 m 2 . La altura de los entrepisos es de 3.40 m. 





Figura 5 Distribution en planta baja y corte del edificio analizado. 



8. MODELO DE PLASTICIDAD 



La estructura analizada tiene deficiencia en la armadura transversal, en consecuencia la 
falla que se va a presentar en los elementos es de corte antes que de flexion, razon por la cual se 
selecciono el modelo de Thom et al (1983) que concentra la inelasticidad en dos resortes uno de 
ellos que trabaja a flexion y otro que trabaja a corte. Estos dos resortes que trabajan en paralelo 
se ubican en los extremos del elemento como se aprecia en la parte inferior de la figura 6. En la 
parte superior de dicha figura se indica el modelo de Giberson (1969) que solo trabaja a flexion 
pero este modelo no es adecuado utilizarlo ya que supone que no se va a presentar la falla por 
corte. En cambio con el modelo de Thom et al (1983) se da la posibilidad de que el programa 
analice y vea en base a la capacidad resistente si se presenta o no la falla de corte. 
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Se ha denominado (EI) a , ( EI) a , ( EI) h a la rigidez a flexion en el nudo inicial, centro de 
luz y nudo final. Por otra parte ( K s ) , (K s ) b son la rigidez a corte en el nudo inicial y final. 



La matriz de flexibilidad f asociada a las rotaciones del elemento, indicado en la figura 6. 
Para el modelo de Giberson esta matriz resulta- 
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Figura 6 Modelos de plasticidad extendida de Giberson (1969) y de Thom et al (1983). 



donde L es la luz libre del elemento. Los terminos ( EI) a , ( EI) a , ( EI) b se obtienen del diagrama 

momento curvatura. Se destaca que al invertir la matriz de flexibilidad se encuentra la matriz de 
rigidez y a partir de esta matriz utilizando la matriz de transformation de coordenadas se 
determina la matriz de rigidez en coordenadas locales y globales del elemento considerando tres 
grados de libertad por nudo (Aguiar 2004) . Para el modelo de Thom et al. (1983) la matriz de 
flexibilidad es la siguiente: 
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siendo K s la rigidez al corte de una seccion del elemento. Existen varios modelos para 

determinar la rigidez al corte que esta sintetizado en Satyarno (2000), uno de ellos es el indicado 
en la figura 7 que relaciona la fuerza de corte con la deformation. 

Park y Paulay (1975) presentan en forma muy extensa la deduction de las ecuaciones para 
determinar la rigidez en el rango elastico K se y en el rango plastico K sp , por unidad de 

distorsion al corte, estas son las siguientes^ 
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Figura 7 Modelo numerico de calculo adoptado para definir la rigidez al corte 



Las variables todavia no definidas son: E c , E s modulos de elasticidad del hormigon y del 
acero, /? factor de forma que toma en cuenta la distribution no uniforme de los esfuerzos de corte, 
para secciones rectangulares /? = 1.2 , 9 sp es el angulo de inclination de la falla por corte, y es el 
angulo del estribo con respecto al eje del elemento. Se puede considerar que 9 sp = 45° para vigas 
y 9 sp = 30° para columnas. 

En Aguiar (2003) se presenta con detalle todas las ecuaciones que conducen al calculo de 
la rigidez de corte, lo importante es destacar que en las relaciones momento curvatura se debe 
incluir la relation corte deformation y de este acoplamiento se determina la verdadera rigidez de 
los elementos. 

En estructuras en la que se sospecha que la falla se va a presentar por corte debido a que 
tienen armadura transversal muy baja es importante encontrar su capacidad sismica resistente 
utilizando un modelo que incorpore el efecto de corte al efecto de flexion y para ello se debe 
seleccionar con detenimiento el modelo de plasticidad adecuado que para el presente caso se ha 
considerado el de Thom et al (1983). 



9. CAPACIDAD SISMICA EN DOS DIMENSIONES 



Utilizando el programa CEINCI2, Aguiar (2002) se obtuvo la curva de capacidad sismica 
resistente de cada uno de los porticos de la estructura, curva que relaciona el cortante basal con el 
desplazamiento lateral maximo en el tope del edificio, para el efecto se trabajo con el modelo de 
plasticidad de Thom et al (1983). En el modelo se considera que los elementos horizontales son 
axialmente rigidos y todos los demas totalmente flexibles. 
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La curva de capacidad sismica resistente se obtiene aplicando la Tecnica del Pushover que 
consiste en aplicar cargas estaticas incrementales en una sola direccion hasta llevar la estructura 
al colapso. En la figura 8 se indica a modo de ejemplo la curva de capacidad sismica del portico 
exterior en sentido X, en ella se aprecia el modelo bilineal que se obtiene de la curva con el criterio 
de iguales areas. En el modelo bilineal se determina una rigidez elastica k e y una rigidez 

plastica k p que sirven para el modelo con tres grados de libertad con el cual se determina la 

curva de capacidad sismica espacial. El limite de separacion de estos dos valores viene dado por el 
punto de fluencia del portico. 



Curva de capacidad resistente 



Portico: P 3 Al [zortapapeies 




Desplazamiento (m) 



Figura 8 Curva de Capacidad sismica del portico exterior en sentido X. 

10. CAPACIDAD SISMICA EN TRES DIMENSIONES 



El modelo numerico de calculo con el cual se determina la curva de capacidad sismica 
espacial en tres dimensiones considera que cada uno de los porticos es un elemento de una 
estructura que une los elementos por medio de una losa que se considera totalmente rigida en su 
piano. Cada portico esta definido por las rigideces k e , k p y los desplazamientos a nivel de 

fluencia d y corrimiento de fallo d tu . En base a estos datos se aplica la Tecnica del Pushover 
aplicando fuerzas laterales incrementales en el Centro de Masa del edificio, Aguiar (2003). 

En la figura 9 se indica la curva de capacidad sismica que se obtiene para la direccion X, 
con el programa CEINCI3 y en la figura 10 se indica la misma curva que se halla con el programa 
ETABS. En las dos figuras se aprecia que la estructura ingresa al rango no lineal con un cortante 
basal que esta alrededor de las 110 T. y que el desplazamiento de fluencia esta entre 2.5 y 3.0 cm. 

Por otra parte se aprecia en las figuras 9 y 10 que el cortante maximo asociado al fallo es 
del orden de las 130 T. y el desplazamiento de fallo que reporta el CEINCI3 es de 5 cm. en cambio 
el ETABS reporta un desplazamiento que esta por los 8 cm. Realmente existe una muy buena 
aproximacion entre las dos curvas. 

Tambien se aprecia en las figuras 9 y 10 la poca capacidad que tiene la estructura de 
incursionar en el rango no lineal es decir la estructura tiene muy poca ductilidad y esto 
fundamentalmente se debe a la deficiencia del refuerzo transversal. 
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Curva de capacidad resistente 




Desplazamiento (m) 

Figura 9 Curva de capacidad slsmica en sentido X del edificio analizado, obtenida con CEINCI3 
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Figura 10 Curva de capacidad slsmica en sentido X del edificio analizado, obtenida con ETABS 



ll.ESPECTRO DE CAPACIDAD 



Una vez que se tienen las curvas de capacidad sismica resistente para los sentidos X, e Y 
se obtiene el espectro de capacidad de la estructura que relaciona el desplazamiento espectral Sd 
con la aceleracion espectral S a para los dos sentidos de analisis. El espectro de capacidad de la 
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estructura es una definition que esta asociada a un sistema de un grado de libertad. Las 
ecuaciones con las cuales se determinan el espectro de capacidad son las siguientes : 



s - D ‘< 

* dj ~ 



Sa , = 



Yi <t>n 

F, 



1 a, Mj 



(5) 

( 6 ) 



Se ha anadido el subindice j a las variables del desplazamiento lateral D t y del cortante basal 
F para indicar que para cada punto de la curva de capacidad resistente se debe hallar el 
respectivo desplazamiento espectral S d y la aceleracion espectral S a . Las variables no definidas 
son: el factor de participation del modo 1, Ct l el factor de participation del modo 1 en el 

cortante, M T la masa total del edificio y (j) tX el maximo valor del primer modo de vibration para 

el punto j, normalmente se obtienen los modos de vibration de tal manera que el mayor valor sea 
la unidad. 




Figura 11 Espectro de capacidad y modelo bilineal de la estructura analizada en sentido X. 



En la figura 11 se indica el espectro de capacidad de la estructura analizada en sentido X 
y el correspondiente modelo bilineal definido por un desplazamiento de fluencia D ly = 1.65 cm. 
Para el sentido Y el espectro de capacidad es similar. 
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12. DESEMPENO ESPERADO 



Se utilizo el Metodo del Espectro de Capacidad, para determinar el punto de desempeno 
de la estructura ante cada uno de los sismos de analisis y se encontro que para la direction X, el 
edificio no es capaz de resistir el sismo raro y el muy raro. 

Para el sismo frecuente la estructura trabaja en el rango inelastico lo propio sucede con el 
sismo ocasional. En la figura 12 se indica el desempeno de la estructura en sentido X ante estos 
eventos. 



Desempeno sfsmico 




Sd (m) 



j Al gortapapeies i| Frecuente Qcasisr.o' Raro Muy raro 

Figura 12 Desempeno de la estructura en sentido X. 



Cerrar 



Para el analisis sismico en sentido Y, la estructura los desplazamientos maximos 
esperados son 2.16 cm., 2.77 cm., y 4.38 cm., para los sismos frecuente, ocasional y raro. El muy 
raro no resiste la estructura y las demandas de ductilidad son 2.14; 2.77 y 4.38. 



13. REFORZAMIENTO DE LA ESTRUCTURA 



Interesa ver cual es el desempeno que va a tener la estructura si se cambia la calidad del 
hormigon a 21 MPa y considerando los siguientes casos : 

• Los estribos que se utilizan en vigas y columnas estan espaciados a 12 y 16 cm. como se 
tenia inicialmente. Los estribos de columnas son de 8 mm de diametro como se tenia 
inicialmente. Unicamente se ha cambiado el diametro del estribo de vigas de 6 mm. a 8 
mm. Los resultados que se obtienen para este caso se indican en la figura 13. Notese que 
ahora la estructura si es capaz de resistir los cuatro eventos sismicos en la direction X. 

• Se mantienen los estribos indicados en el caso anterior pero adicionalmente se coloca una 
bincha (estribo suplementario) en las columnas de tal manera que tienen mas 
confinamiento y ductilidad. Los resultados del desempeno se indican en la figura 14, 
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nuevamente para la direccion X. Como era de esperarse la estructura se comporta de 
mejor forma. 



Desempeno sismico 




Sd (m) 




A^ortapapeie^^^cuent^D^siona^a^^u^ra^ 



Cerrar 



Figura 13 Desempeno sismico de la estructura en sentido X con estribos de 8 mm. en vigas. 




Figura 14 Desempeno sismico de la estructura en sentido X con estribo suplementario en 

columnas. 
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14. CONCLUSIONES Y RECOMENDACIONES 



Se ha presentado un estudio de peligrosidad sismica para la ciudad de Salta con el cual se 
han obtenido cuatro espectros para los sismos- frecuente, ocasional, raro y muy raro que tienen 
periodos de retorno de 43, 72, 475 y 970 anos. Por otra parte se analizo el desempeno que tendria 
una edificacion de tres pisos de hormigon armado construida en el 2000 ante los eventos 
mencionados. Del estudio realizado se desprenden las siguientes conclusiones : 

• Las aceleraciones maximas en roca para la ciudad de Salta para los sismos: frecuente, 
ocasional, raro y muy raro son 0.12 g, 0.14 g, 0.29 g, y 0.37 g, respectivamente. 

• La peligrosidad sismica de la Provincia de Salta es considerable razon por la cual se deben 
tomar mayores precauciones en el diseno sismico de las estructuras habida cuenta que la 
edificacion analizada no fue capaz de soportar el sismo raro y muy raro. Por otra parte se 
observo que la estructura no tiene suficiente ductilidad. 

• Se recomienda que todas las edificaciones de la Provincia de Salta sean construidas con 
un hormigon que tenga una resistencia a la compresion minima de 21 MPa. De igual 
manera que se tengan mayores precauciones en el detallado del refuerzo transversal. El 
diametro minimo de estribos en vigas y columnas debe ser por lo menos 8 mm. 

• A1 mejorar el detallado del refuerzo transversal se mejora notablemente el 
comportamiento sismico de la estructura, se incrementa la resistencia y la ductilidad del 
edificio. 
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